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Jahrestagung der Gesellschaft für angewandte Mathematik 
und Mechanik 


Die diesjährige Tagung der Gesellschaft für angewandte Mathematik und Mechanik fand 
in der Zeit vom 21. bis 25. April in Aachen statt. Der Besuch war gut, wenn auch nicht so stark 
wie der der vorigen Tagungen. Eine Reihe von Mathematikern aus dem Auslande nahmen auch 
dieses Mal wieder an der Tagung teil. Die Zahl der Vorträge (52) machte es erforderlich, daß zwei 
Vortragsreihen nebeneinander laufen mußten, eine für angewandte Mathematik und eine zweite 
für Mechanik, Hlastizität usw. 

Die Hauptversammlung der Gesellschaft fand am Mittwoch dem 23. April um 16 Uhr statt. 
Nach Ehrung der Toten des verflossenen Jahres gaben der Vorsitzende Prof. Waliher und der 
Geschäftsführer Prof. Görtler den Jahresbericht. Letzterer berichtete über die Mitgliederbewegung 
(z. 2.505, davon etwa 25%, im Auslande), und über den Kassenstand. Auf Antrag der 
Kassenprüfer wurde Entlastung erteilt. Um die Hauptversammlung künftig von den vielen Einzel- 
vorträgen zu entlasten, wurden Fachausschüsse gebildet. Diese sollen im Laufe des Jahres Kollo- 
quien abhalten, über die dann auf der Haupttagung nur berichtet werden soll. So war es im ver- 
gangnen Jahre mit dem Fachausschuß für die Entwicklung moderner Rechenautomaten in Deutsch- 
land, der eine Tagung in Göttingen abgehalten hatte, über die Prof. Walther berichtete. Dieser 
Ausschuß wird natürlich weiter arbeiten. Weiter wurde beschlossen den Fachausschuß für Quali- 
tätskontrolle umzubenennen in Ausschuß für statistische Methoden. Seinen Vorsitz legte Prof. 
Graf wegen anderweitig starker Belastung nieder. Künftig werden diesen Ausschuß Dr. Stange 
und Prof. Günther Schulz betreuen. Neugebildet wurden: 

I. Ein Fachausschuß für Rheologie unter dem Vorsitz von Prof. Schultz-Grunow, 

2. Ein Fachausschuß für Regelungstechnik unter dem Vorsitz von Dr. Oppelt, 

3. ein Fachausschuß für moderne Anwendungen der Funktionentheorie unter Leitung von 

Prof. Buchholzund Prof. Ullrich. 

Für die aus dem Vorstandsrat ausscheidenden Herren E. Schmidt und G. Schulz 
wurden die Herren Sauer und Schultz-Grunow gewählt. Der Geschäftsführer Herr 
Görtler wurde wiedergewählt. Das Amt eines Kassenprüfers übernahmen die Herren Prof. 
Tollmien und @. Schulz. Der Mitgliederbeitrag bleibt wie bisher 5 DM. Nächste 
Tagungsorte werden München 1954, Berlin 1955 und Stuttgart 1956 sein. Von dem vor- 
handenen Kassenbestand sollen 1000 DM zur Nachwuchsförderung verwendet werden, z. B. 
als Reisekostenzuschuß für jüngere Mitglieder zum Besuch der Tagungen. Ferner soll bei der 
Geschäftsstelle eine Vermittlungsstelle eingerichtet werden, auf der sich stellensuchende Mitglieder 
melden können. Sie wird so bei Anfragen aus der Industrie usw. passende Bewerber empfehlen 
können (s. Nachrichten in Heft 7, Seite 252 dieses Bandes). 

Natürlich war auch auf dieser Tagung für Geselligkeit und Gelegenheit zur Aussprache 
gesorgt, so veranstaltete die Stadt Aachen einen Empfang im Kurhotel Quellenhof, in der Mensa 
academica fand ein Unterhaltungsabend des Aachener Studentenkabarett „Pünktchen‘‘ statt. 
Ferner beteiligte sich ein großer Teil der Tagungsteilnehmer an den eigentlich für Damen vor- 
gesehenen Besichtigungen des Domes und des Rathauses, sowie an der Fahrt nach Monschau. Am 
25. April fand zum Abschluß eine gemeinsame Autobusfahrt nach Brüssel usw. statt, die unter der 
kunsthistorischen Leitung von Prof. Dr. Mennicken, Aachen stand. 

Ich glaube im Namen aller Tagungsteilnehmer zu sprechen, wenn ich Herrn Schultz- 
Grunow für die sorgfältige Vorbereitung dieser ausgezeichnet verlaufenen Tagung den herzlichen 
Dank aller Teilnehmer ausspreche. Willers 
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A. Praktische Analysis 


Ein Beitrag zur numerischen Integration bei nicht Sleichabständigen 
Abszissen und zur Berechnung von Kurvenintegralen 
Von Heinz Eltermann in Braunschweig 
Wenn von einer steligen Funktion f(x) bekannt ist, daß sie an gewissen Stellen &,, 
&,°"", %, in einem Intervall (a; b) die Werte f(z,) =f, ® =1, 2, ---) annimmt, und das 
bestimmte Integral über f(x) von a bis b berechnet werden soll, ist es in der praktischen 
Analysis üblich, die folgende Methode anzuwenden. 


Man bestimmt eine stetige Ersatzfunktion f(x), die sich stückweis aus geschlossen inte- 
grierbaren Funktionen, etwa Polynomen, aufbaut, so daß sie mit f (x) die vorgegebenen Funk- 


tionswerte gemeinsam hat, und berechnet das Integral über diese Ersatzfunktion. — Die 
üblichen Ansätze liefern f(x) als lineare Form der Ordinaten f,. Man kann also fin der Form: 
5, n 
(®) = yfn(a) - ne RE - (1) 


schreiben, wobei die Funktionen 9,(x) die ‚‚interpolativen“ Eigenschaften: 9,(#,) =1, 
9,(&,) = 0 für v»+ u haben. 
Für das Integral ergibt sich dann der Näherungswert: 


& 292 us b 
J= | H2)da = 7,9. mit m [6 (du ne a 
a 0 a 
J ist also auch eine lineare Form der Ordinaten f,, deren Koeffizienten g, man als Gewichte 
bezeichnet. 

Diese Gewichte hängen nun von der speziellen Wahl der Funktionen g,(x) ab, also 
nebenbei auch von den Abszissen &9, + * *, %,. Legt man sich aber von vornherein auf eine hin- 
reichend eng begrenzte Funktionenklasse für die @, fest — baut man etwa f (&) stückweis aus 
Parabeln p-ten Grades auf — dann bleibt im wesentlichen nur die Abhängigkeit von den 
Abszissen übrig. 

Für die Gewichte ergeben sich also Formeln, welche die Abszissen &,, * * *, „in mehr oder 
weniger komplizierter Form enthalten. So ergeben sich zum Beispiel für den Fall, daß man 
J(«) stückweis linear ersetzt und « = x, und b = x, setzt, die Formeln: 


1 1 1 
"Tg hı, ga, >, (hr), J=7Z (h, + h,), usw. nt (2) 


wobei stets h, = x, — %,_ı einzusetzen ist. 
Für den Fall, daß man f (x) stückweis durch Parabeln zweiten Grades ersetzt, ergeben 
sich in entsprechender Weise die Formeln: 


1 h 1 KT 
nat 2). art). 
1 2! 
1 h } 
= @ ne) (2 = Fir, Eh (2 a usw. 
2 3 


Bei gleichabständigen Abszissen entstehen hieraus die einfachen Formelnder Simpson- 
schen Regel: 


| 
er 


%=ah, Hr h, g=7 h, g=zh, USW. oe.) 


Im folgenden soll nun eine Methode aufgezeigt werden, die in allen Fällen lineare 
Formeln für die Gewichte liefert. 


Man faßt die Integrationsvariable & als Funktion eines Parameters tauf. Dann wird der 
Integrand gleichzeitig eine Funktion von t, und man hat jetzt das Integral: 


b 2 


7= | 1942 = | Sie): 7 dtzurechnen a a 


Den Abszissen x, darf man nun in recht willkürlicher Weise Parameterwerte i, zuordnen, darf 
also auch speziell t, =» setzen. Für die Funktionen f(t) = f(x(t)) und &= « (t) kann man die 
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linearen Ansätze: 
n 


n n 
Ben en. | ds ıl j 5 al 
Kt) en; Zu 9%» (t) &(t) PR Ly x Yu (t) dt Se > Lu £ Yu (t) U N (7) 
machen, wobei die Gleichungen: 9,(,) = pl) =1; lt) =yplt) =0 für vu 
erfüllt sein müssen. 
Jetzt berechnet man einen Näherungswert für J in der Form: 


- R, nn n 
J=/ 0 0 (t) de an RR = hr EEE RR 
mit 7 
Ai 
Iyu u 9 (t) y(t) dt u Be ee are oa. (9) 
und F 
n 
9 n Gy u Yun . . . . . . . Fr; . . . . . . (10). 


Die Gewichte sind also stets lineare Funktionen der Abszissen. Hierzu einige Beispiele: 

Die Funktion f(t) möge aus Parabeln p-ten Grades, und die Funktion x(t) möge aus 
Parabeln g-ten Grades aufgebaut sein. Es sit=1—=0, 1, =»v,1=t,—=n. Dann errechnet 
man folgende Formeln: . 


17 Beispielen liHingSil. 
N 1 
I=zh, = 5 (ht), usw. 
Dieses sind die schon bekannten Sehnentrapezformeln (3). 


ZB spLeh..n liste 2, 


7 1 1 7 | 1 
9% = 79 hh— 75 h,, = D) (tb), &= iD (h, h;) — iD (tu), 


1 
Ga >, (ha + hy) BE LLSIW.; 


saberspiel! pe2, = 


Li 
B) 1 2 B) 1 
% = 15 hı — 15 RR N — 3 (ei), Go —- 75 (A, + h;) 12 (A,+ Rh), 


2 
(ie (h5+h,), usw. 


AaeisDiel: on 2 2 .gr=2. 


1 1 2 1 1 
%"=5 h, =, h,, = 3 (h,+h), = 2 (ha + h;) — 6 (h, + hu), usw. 


Abschließend soll noch gesagt werden, daß man mit diesen Formeln natürlich auch Kur- 
venintegrale rechnen kann. Hierbei empfindet man die Linearität der Formeln besonders 
angenehm, da sie sozusagen invariant gegenüber Koordinatentransformationen sind. Leider 
erlaubt der zur Verfügung stehende Raum nicht, hierauf näher einzugehen. 


Verbesserte Konvergenz- und Fehlerabschätzungen für die 
Störungsrechnung 
Von Friedrich Wilhelm Schäfke in Mainz 
Die wichtigsten praktischen Anwendungen der Störungsberechnung beziehen sich auf 
Matrizenprobleme, auf reguläre Eigenwertprobleme bei gewöhnlichen Differentialgleichungen 


und auf Fredholmsche Integralgleichungen. Es genügt dabei, sich auf linearen Stör- 
parameter zu beschränken und ein Problem der Form 


Fy+iy+n@y=0 
zu betrachten. [A Eigenwert, u Störparameter, y Eigenlösung (Vektor oder Funktion), F, @ 
lineare Operatoren (Matrizen, Differentialoperatoren oder Fredholmsche Integralope- 
19* 
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ratoren)]. Bei diesen Problemen !) gibt es eine ganze analylische Funktion A zweier Variabler, 
derart, daß A(A, u)=0 für ein Eigenwertpaar A, u charakteristisch ist (charakteristisches 
Polynom bei Matrizen, charakteristische Determinante bei den genannten Differentialglei- 


1 
chungsproblemen) oder derart, daß A Er „)—0 für ein Eigenwertpaar 4 #0, u notwendig 


und hinreichend ist (Fredholmsche Determinante bei Integralgleichungen). 

Diese Tatsachen bilden den Ausgangspunkt der folgenden kleinen Theorie. Sie ist zwar 
in Richtung auf Zulassung auch kontinuierlicher Spektren nicht so allgemein wie die Störungs- 
theorie der Spektralzerlegung für selbstadjungierte Operatoren des Hilbertschen Raumes 
(Rellich, Nagy). Sie hat jedoch dafür folgende Vorzüge. Sie baut auf wenigen für die 
genannten praktisch wichtigen Probleme leicht nachprüfbaren Voraussetzungen auf. Sie um- 
faßt auch nicht selbstadjungierte Probleme; nicht einmal das ungestörte Problem (u —0) 
braucht selbstadjungiert oder hermitesch zu sein. Die Herleitungen erfordern sehr geringen 
Aufwand. Die Resultate sind wesentlich schärfer als die bisher bekannten (vgl. das Beispiel 
am Schluß). Sie beziehen sich hier allein auf die Störung von Eigenwerten, die einfache Null- 

1 
stellen von A(A, 0) bzw. A I ‚ 0) sind. 

Bezüglich weiterer Ergebnisse (mehrfache Eigenwerte, obere Schranken für die Konver- 
genzradien) sei auf meine Arbeit Math. Nachr. 6, 109—124 (1951) verwiesen; dort wird das 
ungestörte Problem hermitesch und der Störoperator beschränkt vorausgesetzt, jedoch sind 
die Ergebnisse zum großen Teil unter den hier zugrunde gelegten weiteren Voraussetzungen 
völlig analog zu erhalten. 

1. Vorbemerkungen, Voraussetzungen. ®R sei ein komplexer linearer 
Raum, f, ,°**, 9, 2, seine-Elemente,; a, ß,- - ,A, u," "reelle oder komplexe Zahlen. In 
N sei ein hermitesches positiv-definites Skalarprodukt, (f, g), erklärt. 

NR wird nicht als perfekt in bezug auf die durch die Norm 


Kal 


gegebene Metrik, also nicht als Hilbertscher Raum vorausgesetzt. 


Bei Matrizenproblemen ist NR ein endlich-dimensionaler komplexer Vektorraum und das Skalarprodukt das 
innere Produkt f-9; bei Diiferential- und Integralgleichungen ist i. A. N die Gesamtheit der stetigen (komplex- 
wertigen) Funktionen eines Intervallsa < x < b und 


b im 
LI =S Sa) ge) de. 


Seien weiter I und ll* lineare Unterräume von #. 


Für Matrizen- und Integralgleichungsprobleme fallen U und U* mit N zusammen. Für die Differential- 
gleichungsprobleme bezeichnen U bzw. I* Gesamtheiten von Funktionen aus R, die gewisse Differenzierbarkeits- 
eigenschaften haben und gewissen Randbedingungen genügen. 


In U sollen die linearen Operatoren F und 6, in Ü* der lineare Operator F* erklärt sein. 
F und F* seien zueinander adjungiert: Für vEU, vEll* gelte stets 
(Fan) lu Bu): ei 1 er En ze, 
Bei den Anwendungen: Adjungierte Matrizen, adj un gierte Integraloperatoren (K* (£, 7) = Kn, &)) und (kom- 
plex-) adjungierte Differentialoperatoren. 
Für uEN sei stets 
elle el). ee 


Dabei ist y (&,n) eine für. 0 <&,n< oo stetige, für O<E&,n< oo positive Funktion, die 
doppelt monoton, 
re SIENA en) - cr. 


und positiv homogen vom ersten Grade, 


veäcm)=iylen - - . (4) 


Bei Matrizen- und Integralgleichungsproblemen kann man stets mit beschränkter Störung auskommen : 
y(&,n)=y"E£. Die Gewinnung von Abschätzungsfunktionen y(£, n), insbesondere für Differentialgleichungspro- 
bleme mit nicht-beschränkter Störung, untersucht eingehend J. Schröder?). Die Ergebnisse werden dort in 
einer für die Anwendungen unmittelbar brauchbaren Form zusammengestellt. 


1ST% 


!) Kleine Erweiterungen, z. B. auf Integrodifferentialgleichungen, sind möglich. 
2) Diss. Hannover 1952. 
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Wir betrachten nun das 
Eigenwertproblem: Es werden Eigenwertpaare A, u gesucht, für die 
Fy+iAytuay=0, yeu 
nichttrivial lösbar ist. 
Dazu seien folgende Voraussetzungen erfüllt: 
A. Es gibt eine ganze Funktion A zweier Variabler, derart, daß A(A, u)—=0 für ein 
Eigenwertpaar 2, u charakteristisch ist. 
B. Für u—=0 gibt es, der Vielfachheit entsprechend gezählt, höchstens abzählbar viele 


Eigenwerte A, und zugehörige Eigenlösungen %,. Die Gleichungen (r durchläuft stets sämtliche 
Indizes von /,) 


Fti.v=0, TURN IN IR AN HEr (5) 
sind nichttrivial lösbar. %,, %n bilden ein normiertes Biorthogonalsystem: 
(Yn 5 Ym) = or Se DE se SUR nee ee eve (8 (6). 


Man hat eine verallgemeinerte Vollständigkeitsrelation O<ß<a<x) 


sy 4. ...0... (0 


für jedes faus R°). 
Statt A. kann zugrunde gelegt werden: *) 
: 1 
A’. Es gibt eine ganze Funklion A zweier Variabler, derart, daß A r .n) —=( für ein 
Eigenwertpaar A+0, u charakteristisch ist. Für w—=0 ist A=0 entweder Eigenwert oder 
Häufungspunkt von Eigenwerten. 
2. Hauptergebnisse. Es wird die Störung eines Eigenwertes 7, untersucht, der 


entweder (Voraussetzung A) einfache Nullstelle von 4(}, 0) oder (Voraussetzung A’) von 0 


1 
verschieden und einfache Nullstelle von A er 0) ST 


Wir führen dann folgende Bezeichnungen ein: 
Dr IA A en ats [8 
” Mm 


sei der Abstand des Eigenwerts A, zuu#=0 vom nächstbenachbarten oder evtl. (A’) von 0. Es 
sei 


7 
mark Zu eu (9 
DA En Bep Fr (9) 
für 
1 
A > A| nn. 
Es ist sicher "R 
5 [Am] 0 <2( ea) 10 
2 = me2l1+ 7) Henne (10) 
Schließlich bilden wir mit der Abschätzungsfunktion y von (2), (3), (4) und a, ß von (7) 
a ı 
em = 7 \a,’ Pr Bene... 0. (1). 


Dann gilt: 


Satz: A(A, u)=0 bzw. A(7 


Am (u) = Am = Am, u+ Am u? <h u 
aufgelöst. Sie konvergiert zumindest für 


lul < om 


: )- 0 wird um A=/#,„, u=0 durch eine Potenzreihe 


und es ist dort 


d 
[Am (16) = Am = = ; 


3) Für hermitesche Probleme gilt yt = „„a=ß=1. 
4) Vgl. die Einleitung. 
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: a n d, ER» 
Zu u mit |u| < 0m gibt es außer A, (1) keinen weiteren Eigenwert mit | — 4, < „Kunde 


Koeffizienten von A„(u) gilt 


2% u 


Ir SE On > er F 


Man hat die Fehlerabschätzung 
I 1/2 
ul mer 
uam 
Im = er Ira; 
Ds wel Hıl1ıSssauze 


Hilfssatz1: Sei für ein zEU, fER und ein A, das nicht Eigenwert oder Häu- 
fungspunkt?) von Eigenwerten zu u —0ist, 


BERTE ı i 
Am (u) ‚Fa Am >% zum k [7 


Fz+3/2=f. 
Dann gilt 
Ill < 5 sup A—Aul- Irall< 5 sup MN: 
Beweis: Aus 
Fe+iz=f, ze, W, 
Fa ENT EU? 
folgt mit (1) 
FAME Ye] We rn a 
und 
(Fz, y)=— Az m)+(} m)=— nt; Y) nee. 
Aus (*) ergibt sich mit (7) 
N <arsup AZ Id 


a? 
= ga SUP A— Al (EJ): 
n 
also die erste Behauptung. 
Ebenso liefert (£) die zweite. 


Hilfssatz2: Ist ), u Eigenwerlpaar, A nicht Eigenwert oder Häufungspunkt von 
Kigenwerten zu u N so gilt 


2 ni er; m | In 
lu B (sup A— 4, sup I) 4 


Beweis; Aus | 
Fy+iy=— uGy, yEU, Yy+0 ® 
und (2) folgt mit Hilfssatz 1 
N 
A— A)’ 


IeyllsrÄlsil, |Fyll) 
also wegen (4) und @y=#0 (A ist nicht Eigenwert zu u + 0) die Behauptung. 


n 


<y(& — uG Ey ” || — u@yl|su 
Aral u@y Fl „6yl| sup 


4.Beweis des Satzes. Zunächst zwei Vorbemerkungen: 


a) Für |u| <om ist| A — Am| = 7 unmöglich. Man hätte sonst nach Hilfssatz 2 


) 

18 
een 
also wegen (3) und der Definition von D,„ erst recht 

& 2 
2,(2,2)>1. 
Om B Y (2 m 
Das widerspricht aber der Definition von o,. 


& E 
Om ß Y dy ‚„ Sup 


N 


5) Das kann nur bei A’ und A = ( eintreten, 
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b) Aus a) folgt: Für |#,| < 0, sind die Funktionen von A 
A (A, U) bzw. A ie A 1) 


nicht identisch 0. Um jedes Eigenwertpaar /, 1 mit |w| < 0m und A(A,, w) —=0 bzw. 


Alt, 1) —0 wird mithin nach den Sätzen über implizite Funktionen A(A, u)—=0 bzw. 
0 


A e > n) —0) durch endlich viele bei u, höchstens endlich verzweigte Funklionselemente A(u) 


mit A(u,) = /, aufgelöst. 
Nun zum eigentlichen Beweis. 


x 1 
Da A, einfache Nullstelle von A(A, 0) bzw. A Be 0) ist, wird um A =, u. —=0 durch eine 


Potenzreihe A,, (u) aufgelöst. Daß diese für |u| < 0, Konvergiert, zeigen wir indirekt. 
Andernfalls müßte auf dem Rande des Konvergenzkreises eine Singularität zw, mit 
Ittu| < 0m liegen. Diese muß eine Verzweigungsstelle endlicher Ordnung sein: denn man hat 


wegen a) |Am(#) — Am| < ze Daher bleibt A,(w#) bei u— u, beschränkt. Es gibt daher eine 


Folge u; > u, für die der Grenzwert A„(1,) > A, existiert. Um },, u, kann nun wegen b) aufge- 


löst werden, A,„(z#) muß einer der zugehörigen Funktionszweige sein. — Da es sich um eine echte 
I r d 
Verzweigung handelt, kann man auf einem anderen Zweige — wegen a) immer mit | — „|< en 


— nach u —=0 zurück analytisch fortsetzen (Beweis wie oben). Dort würde man aber wieder 
)(0) =/„ haben müssen nach Definition von d,; das widerspricht jedoch der schon festge- 
stellten Tatsache der Auflösung um 4 =/,„, u = allein durch eine Potenzreihe. 

Damit sind die vier ersten Behauptungen des Satzes bewiesen. Die beiden letzten sind 
einfache Folgerungen. Die Koeffizientenabschätzung ist nichts anderes als der Gutzmer- 
sche Koeffizientensatz. Die Fehlerabschätzung ergibt sich wieder daraus durch Anwendung 
der Schwarzschen Ungleichung: 


(3 mellar)'< 8 mon 3 (RI. 
k=1+1 k=1+1 k=1+1 \Om 
5.EinBeispiel. Vergleiche. Wir betrachten mit A?als Störparameter die 
Mathieusche Differentialgleichung 
y'(l)+ (A — 2%? cos 2t) y(t)=0 
mit den Randbedingungen der Periodizität 
v0) =ym)., YO) =y a). 
Das Problem ist selbstadjungiert; man hat a—=ß = 1 und kann y(&, n) =2£ wählen. 
m = Pist einfache Nullstelle von cos me 1. Es wird d„ =4; D, wird nicht benötigt. 
Man erhält für die Potenzreihe 
Am (4) — 04 (h?) 


die untere Schranke 0, =1 des Konvergenzradius®). 

Hier ergeben die Abschätzungen von F. Rellich?), B. v. Sz. Nagy?) und 
J.Schröder?°)der Reihe nach die Schranken 

0,125; 0,55 0,738, . 
Der exakte Konvergenzradiusistnach C. .Bouwkamp!). 
1,468 768 52 ° - -. 

Auch bei nicht beschränkter Störung gilt für die Schranken der der Konvergenzradien von 
Rellich(g), Nagy (8), Schröder (o,) und die hier gegebenen (9,): 


80 =20, —@s; G@<Q<m- 

Ähnlich fällt der Vergleich der Abschätzungen der Fehler aus; nur tritt hier das Verhältnis der 
Schranken für den Konvergenzradius offenbar in der (l + 1)ten Potenz auf, auch liefert die 
Verwendung des Gutzmerschen Koeffizientensatzes noch eine zusätzliche Verbesserung. 
Für kleine Zsind jedoch die Fehlerschranken von J. Schröder?) teilweise schärfer. 

6) Vgl. auch meine in der Einleitung genannte Arbeit, S. 121/122. 

?) Math. Ann. 117, 356—382 ($ 2), (1940). 

8) Comment. math. Helvetici 19, 347—366 (1946). 


9) Vgl. Anm. 2). . 
10) Proc. Acad. Wet. Amsterdam 51, 891—893 (1948). 
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Zur Lösung von Potentialaufgaben mit Hilfe des Differenzenverfahrens 
Von Johann Schröder in Hannover 


Die Rechenarbeit bei Anwendung des gewöhnlichen Differenzenverfahrens auf Rand- 
wertaufgaben der (ebenen oder räumlichen) Potentialtheorie läßt sich in vielen Fällen wesent- 
lich herabdrücken. Man kann häufig in einfacher Weise ein reduziertes Gleichungssystem für 
etwa die Hälfte der unbekannten (Näherungs-) Funktionswerte aufstellen und erhält nach 
dessen Lösung die restlichen Funktionswerte aus dem entkoppelten System der zugehörigen 
gewöhnlichen Differenzengleichungen. — Bei Eigenwertaufgaben ergeben sich die (Näherungs-) 
Eigenwerte mit Hilfe der Determinante des reduzierten Systems. — Die üblichen Iterations- 
verfahren führen bei Anwendung auf das reduzierte System oft mit erheblich weniger Schrit- 
ten zum Zie] als bei Anwendung auf das System der gewöhnlichen Differenzengleichungen. 

Eine ausführliche Darstellung folgt in einem der nächsten Hefte. 


Zur Interpolation von tabellierten Funktionen durch 
Exponentialsummen und zur Berechnung von Eigenwerten aus den 
Schwarzschen Konstanten 
Von E. Stiefel in Zürich 

Anläßlich der Untersuchung von Algorithmen der Relaxationsrechnung stellte sich 
heraus, daß die im Titel genannten Aufgaben neben einigen anderen !) einheitlich durch ein 
Rechenverfahren behandelt werden können, das weitgehend die Züge des gewöhnlichen 
Differenzenschemas trägt und sich daher auch für den Einsatz von Rechenauto- 
maten eignet. Die folgende Form des Verfahrens wurde vonH.Rutishauser entwickelt. 

Gegeben sei eine mit dem Schritt A =1 tabellierte Funktion f (x) und fu, fi, fa, - - - seien 
ihre sukzessiven Werte. Es wird dann das folgende Quotienten-Differenzen- 


schema gebildet, dessen allgemeines Bildungsgesetz noch durch die Formeln (1) fest- 
gehalten ist. 


s | 0.4 |p=a4| @ | a \a=-R4rAp=Zr| 
Jo 
Ro 
fi r S 
Q, | % ri 
ip Aı N A, A, Si 
Q2 9% 2 0 
fs hh fi i 1i fi 
Qs 1% 
fi A; f3 
Q, 
IB 
a ok — ferı 1% — ok ok ge gen ge—1 k+1__yk Ak 1 
a Henn Ad ieh - years 0). 
ı 


(Der obere Index 1 wurde in der Tabelle weggelassen.) 


In der Q-Kolonne stehen also die Quotienten von zwei aufeinanderfolgenden 
Funktionswerten, in der A-Kolonne werden die Differenzen dieser Quotienten gebildet. 
Anschließend erfolgt die Bildung einer neuen Funktion f? durch 


f=f+ı' 4, 
mit der analog weitergefahren wird, wobei aber die Differenzen A? gemäß (1) zu A? korrigiert 
werden müssen, indem die auf derselben Zeile stehende erste Differenz hinzugeschlagen wird. 
Die Eignung des gewöhnlichen Differenzenschemas zur Interpolation gegebener Funk- 


tionen beruht bekanntlich auf der Tatsache, daß die (n + 1)-te Differenz eines Polynoms n-ten 
Grades identisch verschwindet. Analog gilt für das Quotienten-Differenzschema 


Satz. Ist f(x) eine n-gliedrige Exponentialsumme 
fl) = aM + aA 4 Eau /n 00 Do 
soverschwindetdie n-tekorrigierteDifferenzidentisch. 


!) Orthogonalisierung von Polynomen, analytische Fortsetzung durch Kettenbrüche, Bestimmung der 
Wurzeln von Polynomen. 
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Dabei soll f(x) eine reelle Funktion sein, jedoch dürfen paarweise konjugiert-komplexe 
4; auftreten. Zum Beispiel verschwindet die zweite korrigierte Differenz von sin ka, wo k 
irgend eine reelle Zahl ist. 

Satz 1 gewährleistet die Extrapolation einer durch endlich viele Ordinaten gegebenen 
empirischen Funktion durch eine Exponentialsumme, indem man in bekannter Weise von der 
als konstant angenommenen letzten Differenz nach links hin rechnend die Funktion aufbaut. 
Dabei fällt auch die Entscheidung, ob eine Exponentialsumme den Verlauf der Funktion 
einigermaßen wiedergeben kann. 

Wünscht man zurPeriodenanalyse also zur Bestimmung der }; aus den gegebe- 
nen Funktionswerten vorzustoßen, so benützt man die folgende Tatsache, die wir zunächst in 
demjenigen Spezialfall aussprechen, wo die A; reell sind und verschiedene Beträge aufweisen. 


. . wi . 
Satz2. Die Wurzeln A, seien nach abnehmenden Beträgen an- 
geordnet: 
Al>lAal> >|. 
Dann strebt die k-te Quotientenkolonne gegen A, also 
rk 
i>o 
Sind zwei konkugiert-komplexe }; vorhanden, deren gemeinsamer Betrag nicht mit dem 
Betrag einer anderen Wurzel übereinstimmt, so äußert sich dies in Schwankungen in zwei 
benachbarten @Q-Kolonnen, aus denen man approximativ das quadratische Polynom 
herleiten kann, das die beiden zugrunde gelegten konjugiert-komplexen /,; als Wurzeln besitzt. 
Es ist aber auch möglich, das Polynom n-ten Grades F(}), welches A,, As, *, A, als 
Nullstellen besitzt, auf rationale Weise aus dem Quotient-Differenzenschema herzuleiten. Zu 


diesem Zweck bilde man folgende zweidimensionale Anordnung von Polynomen pf (2) vom 
Grade k: 
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mit dem Bildungsgesetz: 
Er M=Ap7ı ee pi = 10:5 HAITI IE (3). 


Dann gilt 


Saeza se flzreine n-gliedrige Exponentialsumme Q),so stim- 
men die Polynome p#(A) in der letzten Kolonne überein 
wneındesdeıchdem PolvnomF(ü)swelches A4A,‘‘*, A) als 
Nunstellen hat. 


Außerdem ist 
lim pF(A) 


Eereivnomkuten Grades, das genau die Nullstellen 
erhabesitzt. 

Umgekehrt besteht die hier nicht näher auseinandergesetzte Möglichkeit, aus einem 
vorgelegten Polynom F(}) durch Kettenbruchentwicklung ein Quotienten-Differenzschema 
zu konstruieren, so daß Satz 2 zur Auflösung algebraischer Gleichungen verwendet werden 
kann. 

Wir wenden diese Ergebnisse nun auf das Problem an, die Eigenwerte einer Matrix A 
zu finden, die der Einfachheit halber als symmetrisch und positiv-definit angenommen sei. Ist 
A n-reihig, so seien 

A>%> 7 >An 
die als verschieden vorausgesetzten Eigenwerte und e@,, €, * , €, die orthogonalen und nor- 
mierten Eigenvektoren. Nach Wahl eines Ausgangsvektors 


en Er Ve Pe EEE ur (. )) 
bilden wir die iterierten Vektoren 
Aa he taNao rt min een.) 
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und die Schwarzschen Konstanten: 
= (4a, 0) = ML RE TER in ak A 
wobei das Komma das Skalarprodukt andeutet. 

Fassen wir i als unabhängige Variable auf, so ist (6) eine n-gliedrige Exponentialsumme 
und wir nehmen daher die Schwarzschen Konstanten als erste Kolonne in einem Quotien- 
ten-Differenzschema. Die bisherigen Resultate ergeben 
Satz la, Dieon-te korrigierbe Differenz verschwiınder wo,surs 

Extrapolation der Schwarzschen Konstanten er moyiıcae 
Satz2a. Der k-te Eigenwertist dern Grenzwert derık=ten Ouo- 
tLenben konn r 

Da die erste Q-Kolonne einfach Ra yleig h- Quotienten enthält ergibt sich so eine Aus- 
dehnung des Rayleighschen Prinzips auf die weiteren Eigenwerte. 

Satz3a. Das nach Satz 3 konstruierte Polynom F(Ü)ist.das cha- 
rakteristische Polynom der Matrix 4, 

An Stelle von Polynomen mit der Wurzel A kann man auch rationale Funktionen be- 
nutzen, welche die A; als Pole haben. In der Praxis berechnet man sie als Kettenbrüche. Die 
ersten Glieder solcher Kettenbruchentwicklungen ergeben manche der Konvergenzverbesse- 
rungsverfahren, die bisher zur Berechnung von Eigenwerten durch Iteration vorgeschlagen 
wurden. 

Es liegt auf der Hand, daß das auseinandergesetzte Verfahren auch mutatis mutandis 
auf selbstadjungierte Eigenwertprobleme bei Differentialgleichungen anwendbar ist. 


Bemerkungen zur Matrizentheorie 
Von U. Wegner in Heidelberg 


Das Vorgetragene ist mehr von didaktischem Interesse als es Neues zur Matrizentheorie 
bringt. Die zugrunde liegenden Tatsachen gehen auf Cayley,Silvester,H.Poin- 
care, Waring,EulerundFrobenius zurück. Das Prinzip war, die bei der mathe- 
matischen Lösung von technischen Problemen auftretenden Desiderata, die in den Bereich der 
Matrizentheorie fallen, durch einfache, einheitliche und gleichförmige Operationen innerhalb 
des Matrizenkalküls zu erfüllen und damit zugleich die Legitimität der Verwendung des 
Kalküls zu erhärten. 

A sei eine n-reihige, quadratische Matrix mit komplexen Elementen, 


oa) tr — ai + gr — It... Ei 
ihr charakteristisches Polynom und 
h(e)=1; h()=ee— co; h,()=ER° — 48 +6,;° 
h,(s)e tr — ar io? — ...2,=0(%) 


die zugehörigen Hornerschen Polynome, deren erzeugende Funktion 


o(z) — p(®) = h,(&) #1 +h(a)a02 +: + Mn-ı(%) WET RES (1) 


2—ı% 
ist, und für die 
her1(2) = uhr(e) — (— 1)*latı 
ist. h,(z#) könnte über n fortgesetzt werden durch 
Ant,(7) = art — 0, art IL. En. 
Wegen der Gültigkeit der Ringoperationen im Bereich der A gilt 
h,(A)= H,„= 4* — A140, 4? — ... +@E, 
H,nı= AH,„— (- 1) &HıE 
und auch 


f 
h,.=u— it oa "E60, 
also: : 


12-2075 A a 


% 
Yy 
h,=0; v=1,2,3 4 | 
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liefert rekursiv die &,, &, &,,°'*. Esist dann 

1 | Antr—1 
va) 

wenn die Integration über einen Kreis in der komplexen }-Ebene streckt wird, der alle Wurzeln 


Ay Ay Ag" '",A, von p(A)=0 enthält. Man setze in (I)J2z>AE, x > A(AE und A sind ver- 
tauschbar) dann ist: 


-V=3,; Re ashun. Aaklar » = (II). 


26% 
Ö 


PM) E— g(A)=(E Hr) el ENEREECTT, 


Aus 
po(A)E=F(A;))-AE— A), 
wobei F (A; }) die adjungierte Matrix zu} E — A ist, folgt durch Subtraktion 


py(A)=0 Vene y=Bealation)e a 2 ana (1), 
n—1 
A De N ee bie sales nie 2). 
—=0 
Aus (III) folgt bei Beachtung von 9(A)= 0, falls g a)=_ a, 4” eine ganze Funktion in Aist: 
ee a [ .4(A) BED Ei [ Ani 20) = ) 
g 2m ” AE—A = "2ni R a) Arm | 
Syn z 1 De r H. | ze) 
u: Oalrr Br ( )- ” DA 2) 


d. h. man erhält die explizite Darstellung jeder ganzen Funktion von Ai indenHornerschen 
Polynomen 4, mit Hornerschen Zahlen &%. (Wichtig für Systeme linearer Diffential- 
gleichungen, wobei g (A) = et ist). Setzt man in (III) A =4,, so folgt 


(0)= (E— 4.3 8,0) Se a SO 


n—1 
Jede Spalte der Matrix S' HA" *ist Lösungsvektor £ von: 
0 
(,E — A)t=0. 


Die Hornerschen Polynome H, können sehr leicht berechnet werden ohne Kenntnis 
der c;. Letztere bekommt man bei der rekursiven Berechnung der HZ, zugleich mit. Man setze 
in (I) «=, und addiere über £. Es ergibt sich: 


ont (nn — 12? + (n — 2) 52m — dm 


Din De — near gan [2 Be sn) re 
an—3 + Dr h(2,) 4 na) + ed [Dr a) - De,) 4 


Der Koeffizientenvergleich gibt: 


Spur (AH,-)= 5 Alu A TE en ut (5). 


Somit ist: 
H,„=AH,ı,— (-1)1cE=H,„=AH,n,— - _ Spur (4 ER) SE EPAE MER (5) 
und H,=E. 
Pers NET ‚26 
d.h. insbesondere 


Beer d- ge n de Spir (AH) SDR. 1.7) 


1) U. Wegner, Bemerkungen über das Li essche Verfahren zur Erzeugung endlicher linearer Transfor- 
mationen (Berichte der matb. -phys. Klasse der sächs. Akad. der Wissenschaften zu Leipzig, 1933. 
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Da H, = (0) ist, so folgt auch: 
1181 
H„—1A= = Spur (AH„-ı)'E 
d.h. falls Spur (AH HH) 05h 


(Iyrrie 


—1 — ne - z 
z Spur (A H„—ı) 


® VE . . . . . . . . . . . . (8). 
Setzt manin (I) z>E, 2-4 sofolgt: 


gl)(E-A= IH, 


v-=0 
und weiter: 
—1 n—1 
(1)(E +A)(E— A) Se, 
yo v— 
oder 
N 
\T x 
Be B+ 2 M.4+2 AH,_ı 
ren = 2 er Weatunteshe, Tas 
1 = = - Spur (AH,) 
da. 


Eat at Fair Et 


ist. Nun ist nach (5) 5 ne — N, die Summe der Matrizen, die die linke Seite von (5) aus- 


*— = 


machen. A H,_ı= N, die Summe der Matrizen, die an erster Stelle auf der rechten 
—ı 


Seite von (5) stehen und 
s— — De — Spur (AH, 1) 
x»=1 
die Summe der 2. Glieder der rechten Seite von (5), so daß 


Deu pr ee) 
=7 7 Tre) 


E+A : » 
ist. Die charakteristischen Wurzeln 4? von B= En haben aber die folgende Eigenschaft: 
Dann und nur dann ist “ar 


Realteil (44) 03: , falls #421 <15ist.2) 


AP araz < 1 kann aber sehr einfach iterativ festgestellt werden und damit kann die Stabili- 
tät der Wurzeln von A durch sehr einfache, und gleichförmige Rechenschritte aus (5) und (10 
“(mit der Feststellung |A# |ya, < 1) entschieden werden. Die .Darstellungen (2), (3), (4), (6), 
(7), (8), (9), (10) sind für die praktische Verwendung bei Benutzung von elektronischen Ma- 
schinen äußerst wichtig. Mittels der Berechnung von |4|y,, und (4) können die Wurzeln 
von A nacheinander berechnet werden. 


?)U.Wegner ,Remarque sur les valeurs propres des matrices, Comptes rendus des seances de l’Acadömie 
des Sciences t. 288, Paris 1949. 


Ein ausführliches Referat des Vortrages: 


R. Ludwig, Die Möglichkeiten von Iterationsverfahren zweiter und dritter 
Ordnung für Gleichungssysteme wird in einem der nächsten Hefte dieser Zeitschrift 
erscheinen. 
Referate der Vorträge: 
W. Richter, Koordinatentransformationen mit Hilfe eines Fluchtlinien- 
nomogrammes und Anwendungen auf die graphische Lösung von Differentialgleichungen, 
A.Walther, Die deutsche Entwicklung von Rechenautomaten im letzten Jahr 
sind nicht eingegangen. 
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B. Reine und angewandte Mathematik 


Ein Einschließungssatz für die charakteristischen Zahlen 


allgemeiner Matrizen-Eigenwertaufgaben 
Von Helmut Bartsch in Hamburg 
Für die charakteristischen Zahlen der allgemeinen Matrizen-Eigenwertaufgabe 
As—HBo 
gilt der folgende Einschließungssatz: 

A und Bseien n-reihige quadratische, hermitesche Matrizen und B überdies positiv-defi- 
nit. Mit einem Vektor & bildet man die Vektoren y=A zund z—= B x, sowie die Quotienten 
entsprechender Komponenten dieser Vektoren: 

Y; 

=, AR 

27 
(y»7=12,°**, n, sind die Komponenten von y usw.). Der Vektor z ist bis auf die folgenden 
Einschränkungen beliebig wählbar: Wenn z2;=0 ist, muß auch y;=0 sein und der ent- 
sprechende Quotient g,; wird ausgelassen. # darf nicht der Nullvektor sein. Liegen die Quo- 
tienten g;in der komplexen Zahlenebene in einem Kreis mit dem Radius v, so enthält der kon- 
zentrische Kreis mit dem Radius R=r Pmaz mindestens eine charakteristische Zahl der 

min 

Eigenwertaufgabe Av =x Bx, d.h. eine Lösung der Gleichung det.(4— x B)=0. Dabei 
sind Byaz und Bin die größte und kleinste charakteristische Zahl der Matrix B, d.h. die größte 
und kleinste Nullstelle der Gleichung det. (B— ßE) =0. (E = Einheitsmatrix). 

Der Beweis dieses Satzes wird voraussichtlich im ‚‚Archiv der Mathematik“ Bd. 4, Heft 2 
erscheinen. Man kann einfache Beispiele angeben, bei denen der Kreis mit dem Radius R der 
kleinste zum erstgenannten konzentrische Kreis ist, der eine charakteristische Zahl enthält. 


Anwendungen zahlentheoretischer Abschätzungen bei 
numerischen Rechnungen 


Von O. Emersleben in Berlin-Zehlendorf 


Einfachere, alt bekannte Eigenschaften der ganzen Zahlen 
werden laufend benutzt, z. B. die Eindeutigkeit der Zerlegbarkeit einer ganzen Zahl in Prim- 
faktoren bei Addition gewöhnlicher Brüche, wenn diese auf einen gemeinsamen Nenner ge- 
bracht werden. Weit weniger bekannt istde Anwendbarkeit neuerer Ergeb- 
nisse der Zahlentheorie, insbesondere außerhalb der Mathematik selbst, so in 
der Physik der Kristallgitter, wie bei Berechnungen des Potentials, vorzugsweise der 
Gitterenergie, von Kristallgitteranordnungen. 

Eine hierbei des öfteren auftretende zahlentheoretische Funktion ist r,(n), die Anzahl 


der Gitterpunkte auf einer p-dimensionalen Kugel vom Radius Yn um einen Gitterpunkt als 
Mittelpunkt oder — arithmetisch ausgedrückt: die Anzahl der Darstellungsmöglichkeiten der 
natürlichen Zahl n als Summe von p Quadraten ganzer Zahlen. 

Hierbei wird man zunächst zur Anwendung zahlentheoretischer Identi- 
täten geführt. So folgt aus der bekannten Tatsache, daß die Größe von r,(n) lediglich von 
den ungeraden Teilern der natürlichen Zahl » abhängt, die Identität 


HERAN ACT 2 RE 57 
auf der einige mathematische, auch in Physik und Technik anwendbare Ergebnisse beruhen: 
In einer schachbrettartig alternierenden quadratischen Netzebene erhält man mittels (1) für 
die Summe der reziproken Abstandsquadrate eines Gitterpunktes von allen übrigen, nach 
wachsender Entfernung summiert, den geschlossenen Ausdruck ?) 


ee) 2 192 = — 2.177586 .. aan ‚ur. (2), 
n=1 2 


1) 0. Emersleben, Die elektrostatische Gitterenergie eines neutralen ebenen, insbesondere alter- 
nierenden quadratischen Gitters. Z. Physik 127 (1950), 588—609, insbesondere S. 609, Ziff. 14 des Anhangs. 

2)0.Emersleben, Einige Identitäten für Epsteinsche Zetafunktionen 2. Ordnung. Math. Ann. 121 
(1949), 103—106, Gl. (1) u. (3). Vgl. auch O0. Emersleben, Über einige bei technischen Berechnungen vor- 
kommende unendliche Reihen regelmäßiger (gitterförmiger) Punktanordnungen. Die Technik 4 (1949), 479 —480 
und (Ergänzung u. Ber.) 5 (1950), 606. 
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während es früher einige Mühe machte, hierfür einen Wert mit geringerer Stellenzahl zu er- 
halten >) ®). 

Dieser numerische Wert interessiert außer bei Fragen der Kristallphysik !) auch bei Be- 
handlung der Torsion zylindrischer Stäbe gewisser Querschnitte ®) sowie von Flüssigkeits- 
strömungen ?). 

Auch bei p > 4 kennt man anwendbare allgemeine Identitäten für r,(n) ®). 

Auch Abschätzungen von Größenordnungen über das Wachstum zahlentheoretischer 
Funktionen können bei numerischen Rechnungen vielfach mit Erfolg Verwendung finden, so 
Abschätzungen der Funktion r,(n). 

Als erstes Anwendungsgebiet hiervon sind Konvergenzuntersuchungen 
zu nennen. Von einer in der physikalisch-chemischen Literatur auftretenden Reihe, die dort 
mitunter sogar als Definition der wichtigen Madelungkonstante des Steinsalzgitters 
dient, läßt sich mittels zahlentheoretischer Abschätzungen beweisen, daß diese Reihe nicht 
einmal konvergiert. 

Dabei handelt es sich um die Reihe: 


< r.(n) 6 12 8 6 
Aysca = TE ale ren 3.ORE 
urn 2 Yin‘ Hl erTee Sr abe 


die jedoch nicht konvergiert, da ihr allgemeines Glied nicht gegen O strebt, weil nicht: 
r,(n) = o(m}P) °). 


Die Reihe 
< r.(n 0 j 
> (— 2% 22 — oz (s), SEO! 1. 4 Dia EIER ERFER (7), 
n=1 


2 > B i . 
konvergiert daher nicht für s—= 1. Sie konvergiert für 0o> 52 (Wert für s=3 siehe nach- 


folgend Gl. (11).) Für 1<o <Z ist die Frage der Konvergenz noch nicht entschieden; je- 
doch ergibt sich der Wert der Madelung- Konstanten als analytische Fortsetzung der 
Reihe (7) über die Konvergenzhalbebene hinaus, für s=1; also Ayacı = — oz (1). 
Auch bei dieser funktionentheoretischen Behandlung kann vorteilhaft von zahlentheoretischen 
Eigenschaften Gebrauch gemacht werden ?). 

Weitere Anwendungen zahlentheoretischer Abschätzungen ergeben sich bei Fehler- 
abschätzungen, insbesondere bei Restgliedabschätzungen unendlicher Reihen. 

Die erstrebenswerte Angabe brauchbarer Fehlerschranken bei numerischen Berech- 
nungen bietet in vielen Fällen nicht geringere Schwierigkeiten als die Bestimmung eines 
Näherungswertes selbst. Jedoch ist eine einwandfreie Kenntnis der Fehlergrenzen mitunter 
eine dringende Notwendigkeit, z. B. bei den bereits erwähnten numerischen Berechnungen der 
Potentiale von Kristallgittern, von denen gewisse asymptotische Werte neuerdings bei einer 
verteilhaften Rechenmethode mit sehr großer Genauigkeit gebraucht werden. 

Bei ihrer Berechnung als Werte Epsteinscher Zetafunktionen, bei denen Potenzen 
‚bestimmter (ein Entfernungsquadrat darstellender) quadratischer Formen auftreten, hängen 
die Restgliedabschätzungen unendlicher Reihen von r,(n) ab. 

Diese Abschätzungen der Anzahl von Gitterpunkten kennt man als interessante Frage- 
stellungen in der additiven Zahlentheorie. Jedoch findet man nicht alles fertig vor, was man 
hierfür braucht. So hat sich die Zahlentheorie vielfach nur für de Größenordnung 
des Wachstums gewisser zahlentheoretischer Funktionen interessiert. 

Die Abschätzung r,;(n) = O(n!/?+®) für jedes e>0 reicht für Restabschätzungen jedoch 


nicht aus. Füre= „ist 6 eine Konstante, wie man sie hierbei benötigt, sogar die bestmög- 
liche, da 


0<r,;(n)<6n für ‚alle: 1.020 olireui a6 rue 
und für n == 1 sowie n — 2 in (8) rechts das Gleichheitszeichen gilt. 


») O. Emersleben, Gitterpotentiale und Zetafunktionen. Diss. Göttingen 1922; Anhang. 

4) O. Emersleben, Das Darcysche Filtergesetz. Physik. Z. 26 (1925), 601—610, insbes. S. 605. 

5) O0. Emersleben, Über zwei Epsteinsche Zetafunktionen 4. und 8. Ordnung. Ber. Math. 
Tagung Berlin 1953 (im Druck). 

6) O. Emersleben, Über die Konvergenz der Reihen Epsteinscher Zetafunktionen. Math. 
Nachr. 4 (1950/51), 468—480, insbes. $ 2. j 

?) Bei Durchführung numerischer Rechnungen dieser Art ist es nützlich, einige zahlentheoretische Eigen- 
schaften von r,(n) zu kennen, z. B. daß r,(n) = 0 dann und nur dann, wenn n =41(8m +7) mitZund m > 0 
ganz. 
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Als Beispiel soll die Restgliedabschätzung der Reihe 


ai rn Elhtar% Ken Breite 3 und da) 


durchgeführt werden, die für p = 3 interessiert, zumal man für p = 1, 2, 4 und 8 geschlossene 
Ausdrücke kennt °). 


Unter Anwendung des Eulerschen Integrals 2. Gattung und von Thetatransformationen 
erhält man Darstellungsmöglichkeiten von (8a) durch gut konvergente Reihen wie 


az 1 )= 162 D nn) Eil-(8n + - > yo -F(z10)) 
N a \ 

| 

| 


n—=0 


ß (9). 
Seele) ” u 
Hrald. Te : 
Hierbei bedeutet g5(n) die Anzahl der ER nr von nZ0in der Gestalt 
ü 1 
Kill 8: 2 mit k,2>0 
dabei ist 
8 g4,(n)=r,(dn +3). 
Ferner bedeutet 
T T 
Kr Er 
F(x)=— |e | — dt 
Ve) Ya) Ye 
0 0 
das Gaußsche Fehlerintegral und 
un —t 
Bi (— ») = fe di 
den Integrallogarithmus. 
Für diese erhält man bei & > 0 leicht die Abschätzungen 
ee? 
<1 Hals EEE EERITETETTE IE 
und 
Be 


Explizite wurde in (9) die erste Reihe bis n=2, die zweite und dritte bis n = 10 be- 
rechnet. 
Für die Reihenreste ergibt sich: 


0< 32, 1m) Eil— (8n +3)] 


e-8n+3) 


<I,n@n +3) a wegen (8) 
= — (8243) > 
($n +3)- (wa Pit: 
en‘ een DB 
6.027 ex? 1 
SE ar TA : e EEE OR SE ‚ 
rege Ba 12 10 (9a) 


Hier gelangt man also mit der Abschätzung (8) auf eine bequem summierbare geo- 
metrische Reihe. 
Ähnlich ist für die dritte Reihe, wiederum wegen (8) 


ze 


Saul % "or Dur ee eg. 5) (de): 


n=il uaals!" 4 


n=11l 
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Bei der zweiten Reihe ergibt sich: 


If 26 © n 
m Ta (N T 6n 2 _n- 
Zr ahnen 
In=11 A ee / Yn 
_— nn? 
12 2 ut 9 2 EEE 
= —gR »3= 2 12 _3% Br n 
Eh 11 n=1l 
12 -17 
2 ee: 
= 11 IT 1 e = ),36 10 . 


Unter Berücksichtigung der bei den vorangehenden Reihengliedern aufgelaufenen 
Fehlerabschätzungen ergibt sich: 
2 \\| (8) = 3,238 624 927 + 1,25-10-”, d.h. 
2| 


8) — —- 3,238625 . . (11); 


vgl. auch 5), Tabelle. 

Man sieht, daß die Fehler der Restgliedabschätzungen völlig zurücktreten gegenüber 
dem maximalen Gesamtfehler, der zu etwa 80% auf den ersten beiden Gliedern der ersten 
Reihe beruht. 

Bei der zweiten Reihe war die Abschätzung r,(n) = O(n) nicht voll ausgenutzt worden. 


10| 
Bei Berechnung von Ayaı = — ®9Z 14) (1) gelangt man jedoch auf ein Glied der Größen- 
ordnung Era — rn D) ‚ für die man wegen (10a) mit der Abschätzung r,;(n) = O(n) Ö 


gerade auf eine geometrische Reihe und damit ebenfalls auf eine handliche obere Schranke 
gelangt. 


Konforme Abbildung eines Rechtecks mit Halbkreiskerbe 
Von Helmut Epheser in Hannover und Friedemann Stallmann in Gießen 


Die vorliegende Arbeit geht auf nicht mehr veröffentlichte Ansätze von Herrn vonKop- 
penfels zurück, die von den Verfassern genauer untersucht und bis zur numerischen Be- 
herrschung durchgeführt wurden. Ein Teil dieser Untersuchungen betr. die konforme Abbil- 
dung eines Parallelstreifens mit Halbkreiskerbe wurde bereits veröffentlicht (Arch. Math. III 
(1952) S. 276—81). Ein entsprechendes Verfahren für das Rechteck mit Halbkreiskerbe soll 
hier angegeben werden. 

Wir betrachten das elliptische Integral 3. Gattung: 


z 
n(#)= const - je Bi ei 
& — a ya — &2)(1 — RE) 
es 6, k, reell, 0<k<]1) 
Durch dieses Integral wird das Gebiet J(2)>0, |#|>o, e<1 auf ein mit einer Kerbe ver- 6 
sehenes Rechteck abgebildet, sofern nur die singulären Stellen x — + Ya, & = + Yeim Inneren 
des Kreises |x| <oliegen. Die Größe und Form der so entstehenden Kerbe werden durch die 
vier Parameter a, c, k, e bestimmt, und durch geeignete Wahl dieser Parameter kann man er- 
reichen, daß sich die Kerbe einem vorgegebenen Kreisbogen gut anschmiegt. Das Gebiet 
J(2)>0, |e|>o läßt sich dann durch elementäre Funktionen auf den Einheitskreis oder ein 
entsprechendes Normalgebiet abbilden. 
Zur Bestimmung der Parameter schreiben wir das Integral (1) als Thetareihe. Wir setzen: 


- 


(5 Flhg, a=sing g=e ®.)) 


Es wird dannn=v-+ yw, wo w die ER ENtMIOENIBE in vund g bedeutet: 


w =lo sin De in 2nv 
u = u gr N { 
oder (2) 
—_ _ gJo ger Gin 2no’ 
u ee en ” sin nv 


!) F(k, p) sowie K, K’ sind das unvollständige bzw. die vollständigen Integrale 1. Gattung in der Le g en- 
dreschen Normalform. 
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y und a bzw. a’ sollen dabei reell und >0 sein. Sie hängen in einfacher Weise mit den bisher 
betrachteten Größen a und caus (1) zusammen. 


Die leben J(x2)>0 wird durch v(x) auf ein Rechteck der Breite x und der 
e n 
Höhe h = SEK abgebildet. Als Bild des Kreises |z| =o erhält dieses Rechteck eine Kerbe, 


deren Fußpunkte bei +v, und deren Scheitel bei iv, liegen möge. Es mußdann a < v, bzw. a’ <v, 
sein. Die Bilder von v,, %v, in der w-Ebene seien w(v,)— — w, und w(i%,)—= iw,. (Bild 1.) In der 
n-Ebene erhalten wir dann ein Rechteck mit einer Kerbe, deren Fußpunkte bi+ 1 —+(v, — yw,) 
und deren Scheitel bei ?73—=t(%-+ y ws) liegen. Hat die vorgegebene Kerbe den Halbmesser R 
und sollen Fußpunkte und Scheitelan den richtigen Stellen liegen, so müssen die Gleichungen 


Ne il. SR a a RER ARERE 6) 
erfüllt sein. Daraus bestimmt sich y zu: 


NR 
2 W, + Ws n 


Das Rechteck in der n-Ebene ist so normiert, daß es die feste Breite mhat. Seine Höhe ergibt 
sich zu: 


h= ee nn me (BD). 


SS 
a 
= 
oO 
S 
| 
li 
a 


Bildi 


Die Gleichungen (3), (4), (5)können dazu dienen, bei vorgegebenem R und h und bei einem — 
innerhalb der angegebenen Grenzen — zunächst beliebig zu wählenden Wert von a” die 
Größen 0, y und k zu bestimmen. Wollen wir jetzt noch den günstigsten Wert von a“ ermit- 
teln, so müssen wir noch untersuchen, welchen Einfluß dieser Parameter — unter der Voraus- 
setzung, daß Gl. (3), (4), (5) erfüllt sind — auf die Kerbenform hat. Wir können hierzu folgen- 
des sagen: Strebt x > v, so strebt w, > oo, der Einfluß von yw wird daher nur in der Umgebung 
der Fußpunkte wirksam, sonst ist die Kerbe dieselbe wie in der v-Ebene, im ganzen erhält man 
also eine Kerbe mit einer Ausbauchung. Strebt umgekehrt & — v,, so erhält man eine Kerbe, 
die gegenüber derjenigen in der v-Ebene nur in der unmittelbaren Umgebung des Scheitels ver- 
ändert ist, im ganzen also eine Kerbe mit einer Einschnürung. Zwischen diesen beiden Grenz- 
fällen muß dann die gewünschte Kerbenform liegen. Zur praktischen Bestimmung der Para- 
meter sei noch erwähnt, daß die Kerbenform durch das Periodenverhältnis nur wenig beein- 
flußt wird, so daß wir zur Bestimmung von a“ den Grenzfall k =0. d.h. q —=0 zugrunde legen 
können. Wir brauchen dann nur elementare Funktionen im Komplexen heranzuziehen. Die 
so gewonnenen Werte von a2) setzen wir dann in (3), (4), (5) ein und bestimmen die übrigen 
Parameter, wobei wir nur noch mit reellen Größen zu rechnen haben. 

Die Genauigkeit dieses Verfahrens ist — im Verhältnis zum Rechenaufwand — recht 
groß, und es dürfte dadurch den üblichen Näherungsverfahren überlegen sein. Man kann hier- 
mit natürlich neben der halbkreisförmigen auch andere Kerbenformen annähern und so auch 
andere Gebiete der Behandlung zugänglich machen. Erwähnt sei etwa der Hohlzylinder mit 
Keilnuten. 


2) Bei der hier vorliegenden Abbildungs-Aufgabe tritt nur & d. h. die erste Form der Thetareihe für nauf. Bei 
anderen Kerbenformen kann jedoch auch die andere Form der Thetareihe eine Rolle spielen, weshalb sie hier er- 
wähnt wurde. 
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Die mechanische Deutung einer geometrischen Differenzgleichung 
Von Wolfgang Hahn in Braunschweig 


1. Unter einer homogenen linearen geometrischen Differenzengleichung für eine Funk- 
tion f (x) verstehen wir eine Funktionalgleichung der Gestalt 


+) + me fge) + ma) Fee) =0 (1) 


oder der äquivalenten Form 
P,(a) 9" f(@) + Pa-ıla fa) ++ Pl)d fe) + Pa)fe) =0  Q), 
wobei |g| < list und die p;(&) bzw. P,(&) Polynome sind; der Operator ® f(x) ist durch 


df(e) = Eee) fa)=BlHflR)).. 
(—1)% 
definiert !). Die Funktionalgleichungen (1) bzw. (2) sind Verallgemeinerungen der linearen 
Differentialgleichungen, mit denen sie zahlreiche Eigenschaften gemeinsam haben, und lassen 
sich durch Grenzübergänge in diese überführen. Im folgenden wird ein einfaches Problem der 
Mechanik betrachtet, das auf eine Gleichung der Form (1) führt; die Lösungen sind Spezial- 
fälle der in der Literatur mehrfach behandelten Heine schen Reihen. 

2. Ein masseloses, im homogenen Schwerefeld schwingendes Seil sei mit n Massepunkten 
P,=1,2,::-,n) besetzt. Es sei k, der Abstand P; P;}ı, m; die Masse von P;, u;(t) die 
horizontale Auslenkung des t-ten Punktes. Bezeichnet man noch die Seilspannung in P; mit 
S;, so gelten nach bekannten Regeln (im Fall kleiner Schwingungen) die Bewegungsgleichungen 

2), 
m; n —— nn (u; — w-ı) — Tu — u) (=1,23,..,,1 Wem =0) . (4) 


Der Separationsansatz u;(t) = U(t) y; führt auf die Gleichungen 


S 
en ee a a re Ta li) 
S;_ S; . 
Bm; =, Hu) ur W) =, nl 


Darin ist # ein Separationsparameter, die Maximalauslenkungen y; sind Polynome in u, 
und zwar muß man weinen Wert erteilen, der die Auslenkung des Punktes, an dem wir das Seil 
aufgehängt denken, zu Null macht ?). Wir wollen voraussetzen, daß die Massen m; und die Ab- 
stände %k, in geometrischer Reihe abnehmen. Es sei also 


m; —gG Mn gl Vi 0 1 Fe a Zu 


Zur Vereinfachung setzen wir m, =g! (g die Erdbeschleunigung), k, =k/(1 — g). Die 
Punkte sollen zunächst von unten nach oben gezählt werden, so daß P,„,, in Ruhe ist. Die 


Den 
Spannung ist dann gleich (m, + ma + + m.) er . Die Gleichungen (5) und (6) 
führen, wenn man noch ed 


ee 
Tu) = IE m Ya)» = 221, 


r=1 


1 


Il 


einführt, auf 


1— {1 ae 1 — gi-1\2 h 
Tu) = (n— Fe u) Ta) Te) Be) Re (8). 


Das Polynom i-ten Grades in u, das dieser Rekursionsformel genügt, läßt sich unter 
Benutzung der Abkürzung 


1-9, =1l—-2)(1—ge)  (1—-ga) (r=12..), (l-a)=1 (9) 
in geschlossener Form aufschreiben; es ist nämlich, wie man leicht verifiziert, 


_ it) Fi 
KW A ag... 2 do) 


!) Zur Literatur vgl. die Angaben bei W.H abn, Über die höheren Heine schen Reihen und eine einheit- 
liche Theorie der sogenannten speziellen Funktionen. Math. Nachr. 3, 257—294 (1950); W.Hahn, Über un- 
eigentliche Lösungen linea,er geometrischer Differenzengleichungen. Math. Ann. 125, 67—81 (1952). 

?) Vgl. 0. Bottema, Die Schwingungen eines zusammengesetzten Pendels. Jber. Deutsche Math.- 
Vereing. 42, 42—60 (1932). 
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Die Rekursionsformel (8) charakterisiert die 7;(u) als eine Kette von Orthogonalpoly- 
nomen. Wir folgern daraus, daß T7;(u) für alle © reelle und einfache Nullstellen hat. Die 
Veränderliche u ist so zu bestimmen, daß 7,(w) =0 wird. Den n Wurzeln dieser Gleichung 
entsprechen n Hauptschwingungen des Seiles, aus denen sich die tatsächliche Schwingung 
gemäß den Anfangsbedingungen zusammensetzt. 

3. Wir wollen nun annehmen, daß das Seil mit unendlich vielen Massepunkten der 
Gestalt (7) besetzt ist. Seine Länge ist dann 

k 
1-9)? se 


Der ruhende Punkt sei P,. Dann wird ; = (mı + mi42 +) = gi/(l — g), so 
daß (6) in 
i—1 L b 
KITTS TU TUN) Un) 00.0... (2) 


übergeht. Diese Gleichung deuten wir folgendermaßen: sie verknüpft Werte einer Funktion an 
Argumentstellen 2 =g(j =1,2, -- -), ist also eine geometrische Differenzen- 
gleichung. Schreibt man V(x) für y;_,, also x für g’=!, so folgt nach einfacher Umformung 


Nez (RER 2)Vle) FVee)=0... 2.0.3). 
Wir setzen zur Lösung 
i je 
Bere A WE erde nee.» (14) 
j=—o 
und erhalten für die unbestimmten Koeffizienten die Rekursionsformel 
Bi ne graz en Nenn, (ld). 


Im Fall} =0 bricht die Rekursionsformel mit 7 =0 ab. Dann ist (14) eine Potenzreihe, 
und zwar ist, wenn noch A, —=1 gesetzt wird, mit Benutzung der Abkürzung (9) 
ae) 


oo 


Ka) = I Yo okay Bu. (16). 


r=0 

oV(«) 
OA =0 
einen Logarithmus und scheidet für unser Problem aus, desgleichen die Lösungen mit 
beliebigem A. Die Reihe ist ebenso wie (10) vom Typ einer Heineschen Reihe). 
Der Parameter u wird die durch Randbedingung %,—=0, d.h. V (=) -1ı=V(1)=v(u) = Ofest- 
gelegt. Da (16) als Grenzfall von (10) für <— oo angesehen werden kann, folgt aus den Bemer- 
kungen über die Nullstellen von 7,(u), daß V(1) als Funktion von u unendlich viele reelle ein- 
fache Nullstellen 1, Us, ‘hat. Zum Nachweis, daß die entsprechenden Eigenfunktionen 
v(u ®), v(uz%), ein Orthogonalsystem bilden, führen wir die Umkehrung der in (3) defi- 
nierten Operation ein. Ist 9 F(x) =f(z), so ist — Konvergenz der Reihe vorausgesetzt — 


Fa) =F0)—- zg-)IEfd a. 
Wir schreiben daher in Analogie zum Integral 
ee ee (18). 


Mit dieser Bezeichnung gilt die ‚„‚Orthogonalitätsrelation‘ 


Die ‚zweite‘ Lösung von (13), die sich in der Gestalt darstellen läßt, enthält 


1 0 mn 
S vum Y) vn N) AG Weg — 1 (19) 
. jA vun g)vln) man i 
Aus (13) folgt nämlich für die Bildung 
D(2) =v(um 1%) ven ©) — Vltm ©) vet q%) 
die „Lagrangesche Identität“ 
D(g =) — D(x) k ” 
®D EI P— I —— ee m rn ® An % v(ung%) u |. (20). 
= = lin un) Om g2) 


3) Vgl. die Angaben in Anm. 1. 
20* 
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und daraus (19) durch g-Summation zwischen 1 und 1 gemäß (18), wenn man noch v(u„)= 
v(u„)—=0 und v(0)—1 berücksichtigt. Auf die durch (19) nahegelegte Frage der Entwicklung 
einer beliebigen Funktion nach den Eigenfunktionen (16), d.h. nach der Lösung unseres Pro- 
blems bei beliebigen Anfangswerten wollen wir hier nicht eingehen. 

4. Wenn man die Massepunkte so anordnet, daß der Häufungspunkt P, oben, also in 


Ruhe ist, so führt die Rechnung nicht zu der Funktion (16), sondern zu der nach den Aus- 
drücken (9) fortschreitenden Reihe 


I (1 = %), 
w(%, u) = 1% ku), |k m az 
= I Nager Ikal< (21) 
Die Eigenwerte sind die (zwischen 0 und eins gelegenen) Nullstellen von w(0, u); auch hier 
genügen die Eigenfunktionen einer Orthogonalitätsbeziehung vom Typ (19). 

Wählt man in (7) für die Folgen m, und %, verschiedene Quotienten, so geht deren Ver- 
hältnis als zusätzlicher Parameter in die Rechnungen ein. Die Lösungsfunktionen sind eben- 
falls Heine sche Reihen wie (16) und (22). 

5. Wenn man von der diskreten zu einer kontinuierlichen Massenverteilung übergehen 
will, so muß man g gegen eins gehen lassen; vorher ist wegen (11) %k durch k(1 — g)? zu ersetzen. 
Berücksichtigt man die Beziehung lim (1 — g)/(1 — 4) =Yr, so entsteht mit g— 1 aus (16) 


die Reihe 
genen 
— 7er 


Das ist im wesentlichen die Reihe der nullten Besselschen Funktionen, durch die ja 
bekanntlich die Bewegungsgleichungen des schwingenden Seils gelöst werden®). 


Über die Lichtverteilung im Bild eines Linienstückchens bei 
Aberrationsfunktionen 2. Grades in den Pupillenvariablen 


Von Rudolf Krieger in Gießen 


Ein optisches System aus zentrierten Kugelflächen, etwa ein photographisches Objektiv, 
vermittelt von einem Objekt eine mit Aberrationen behaftete Abbildung. Das betrachtete 
Objekt wird als selbstleuchtend angenommen. Die durch das Objektiv vermittelte Abbildung 

wird hier nach der geometrischen Optik be- 
handelt; d.h. alle Interferenz- und Beugungs- 
X erscheinungen werden vernachlässigt. 
Die Strahlen des von einem Objekt- 
‚ punkt P ausgehenden wirksamen Bündels 
werden durch die 3 Größen beschrieben 
(Bild 1): 
l: Abstand des Objektpunktes von der 
Achse, 
Bild 1 (£,n): Maß für die Richtung des Strahls 
im Bildraum. 
Man kann £, n näherungsweise als kartesische Koordinaten einer in der Austrittspupille 
(A.P.) gelegenen Ebene deuten. Die Randstrahlen des Bündels bilden in der A.P. im all- 
gemeinen ein Zweieck, wenn die Strahlenbegrenzung durch die Vignettierung an den Linsen- 


& TE 
u: — 1] annähern. Die bei 


ERBA 


> 


r= 


fassungen gegeben ist. Diesen Rand kann man als Ellipse: jr 
der Näherung verlorengegangenen Strahlen tragen in sehr vielen Fällen wenig zur Helligkeit 
imBild bei und können daher vernachlässigt werden. — Die Aberration vom idealen Bildort hat 
die sagittale Komponente x und die meridionale Komponente %; sie ist im allgemeinen bei den 
verschiedenen Strahlen verschieden; also: 


s=Xl;ö,n)\, y-Ylön. 

Die ‚‚Aberrationsfunktionen‘ X, Y werden als Potenzreihen in den Pupillenvariablen 
&, n geschrieben, wobei die Entwicklung aber nur bis zu den in &, n quadratischen Gliedern 
getrieben werden soll, denn dann ist die Helligkeitsverteilung im Bild eines Linienstückchens 
durch elementare Funktionen darzustellen. Die Aberrationsfunktionen 2. Grades in &, n 
lauten: 

s=SE+2REn, y=Tn+3Km+RE. 
4) Vgl. Anm. 2. 
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Legt man eine Aberrationsfunktion 3. Grades mit beliebigen Koeffizienten zugrunde, so 
muß man zur Darstellung der Helligkeitsverteilung elliptische Integrale 1. Gattung heranziehen. 

Die Lichtverteilung im Bild des Objekts wird durch die Aberrationsfunktion bestimmt, 
also durch die Bildfehler: R (Rinnenfehler), $ (sagittaler Schnittweitenfehler), 7 (tangentieller 
Schnittweitenfehler), X (Koma). 

Als Objekt wird ein Linienstückchen betrachtet. Es wird vorausgesetzt, daß dieses so klein 
ist, daß die Koeffizienten A(l), S(l), T(l), K(l) für die verschiedenen ] des Linienstückchens 
praktisch konstant sind, und daß das Linienstückchen entweder auf einer den Achsenpunkt 
enthaltenden, in der Objektebene gelegenen Gerade liegt (‚‚sagittales‘ Linienstückchen) oder 
daß es eine solche Gerade rechtwinklig schneidet (‚‚meridionales“ Linienstückchen). 

Im Fall des sagittalen Linienstückchens ist die relative Helligkeit H nur eine Funktion 
von zallein: 4=H (x), für die man folgenden Ausdruck erhält: 


a-10-- ((M)2-[(&)m 


wobei das Integral längs dem Teil der Kurve =S&E +2 REn =const. zu nehmen ist, der 


2 2 
in der Ellipse En —=]1 liegt. Beim meridionalen Linienstückchen erhält man für die 


relative Helligkeit 7 = H (y)einen analogen Ausdruck. Die folgenden Betrachtungen beziehen 
sich auf ein sagittales Linienstückchen; sie gelten aber analog für ein meridionales Linien- 
stückchen. 


Die Helligkeit H ist also ein bestimmtes Integral. Die Bestimmung der Integrations- 
grenzen macht Schwierigkeiten, wenn die Integrationsgrenzen die Schnittpunkte der Kurve 
&2 
a2 
chungen 4. Grades ergeben. Man umgeht diese Schwierigkeiten, indem man nicht bei vor- 
gegebenem & die Integrationsgrenzen ausrechnet, sondern indem man die Integrationsgrenze 
als Parameter zur Darstellung der Kurve für die Helligkeitsverteilung vorgibt. Die Kurve 
für die Helligkeitsverteilung erhält man dann wie folgt: 


Es wird eine neue Variable « eingeführt durch: 


u=— | (),®. 


wobei das Integral als Funktion der oberen Grenze betrachtet wird; die untere Grenze kann 

willkürlich gewählt werden. Das unbestimmte Integral «ist danach eine Funktion von x und 
E:u—=U (2). 

Die Integrationsgrenze € wird nun als Parameter aus dem Intervall: — az s<-+xa 

&2 


a? 


2 
%=const. mit dem Rand au, —=1sind, also sich als Wurzeln von allgemeinen Glei- 


2 
vorgegeben, aus + =1 wird n ausgerechnet: n=n(£) und in die Aberrations- 


gleichung: &=X(En) eingesetzt: —=x(£). Setzt man z=x(f) in u=U(z,) ein: 
u—u(E), so erhält man in Parameterdarstellung die Kurve: 2 —=x(E), uw=u(8). 

Aus der Kurve v=u(E), »—= x(E£) kann man die Helligkeit wie folgt ablesen: 

Hat die Kurve x =const. mit dem Rand 2 Schnittpunkte, so gehören zu diesem &-Wert 
2 u-Werte ; deren Differenz ist dann die gesuchte Helligkeit 4. Hat x =const. mit dem Rand 
4 Schnittpunkte, so liefern die zu diesem z gehörenden 4 u-Werte paarweise 2 Strecken, deren 
Summe dann die Helligkeit Hin ist. Entsprechendes gilt für 6 und mehr Schnittpunkte. 

Die Helligkeit H wird also aus einem oder mehreren Paaren krummer Kurven abgelesen; 
ein solches Paar krummer Kurven kann auch eine zusammenhängende Kurve bilden, die auch 
geschlossen sein kann. Nach dieser Methode kann man für alle möglichen Fälle der zugrunde 
gelegten Aberrationsfunktionen 2. Grades die zugehörigen Helligkeitsverteilungen berechnen. 
Ein Überblick über alle überhaupt möglichen Helligkeitsverteilungen erhält man, indem man 
für etwa 20 typische Fälle die Helligkeitsverteilungen zeichnet. 

Eine ausführliche Darstellung erscheint in der Optica Acta, Paris, 
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Allgemeine Betrachtungen über die Randwertbedingungen in der 
Elektrodynamik bewester Körper 
Von @. Lampariello in Rom 
Die Feldgleichungen 
rt E=— BD div®—0 
rt 9 =dD+I divd=o 
bleiben auch für bewegte Körper bestehen, wie H.Minko ws ki zuerst bewiesen hat. 

In ponderablen Körpern aber haben die Größen &, D bzw. ®, 9 keine unmittelbar an- 
schauliche Bedeutung; nur im Vakuum kann man sagen, daß diese Größen die Kraft auf einen 
im Inertialsystem K ruhenden elektrischen bzw. magnetischen Einheitspol sind. 

Die Größe o ist die Ladungsdichte, und $ist der allgemeine Leitungs- + Konvektions- 
strom 

N — Ry = oVv. 
Bei den konkreten physikalischen Problemen muß man auch mit den Bedingungen auf den 
Trennungsflächen der bewegten Körper, die ein verschiedenes elektrisches oder magnetisches 
Verhalten aufweisen oder auf der Grenzfläche der bewegten Körper im Vakuum rechnen. 

Die Untersuchung der Randwertbedingungen führt man am besten durch, indem man 
von der Integralform der Feldgleichungen 


[Bude =— [C*ds [C„de=[H*ds 


ausgeht, in welchen 


Ce =E+VAB HH —=H—-vAD 
B8=-B—-rotvA®B D=-Dd+o9— ro AD 
0=9+%, 


ist. 

Diese Integralform hatte schon H. A. Lorentz vor der Entdeckung der Relativitäts- 
theorie aus der Elektronentheorie hergeleitet. 

Wenn durch die Trennungsfläche der bewegten Körper, die aneinanderstoßen, die Ge- 
schwindigkeit v stetig ist, haben wir die folgenden Randwertbedingungen: 

1) Die Stetigkeit der Normalkomponente von 9; 

2) der Sprung von ®, ist gleich der Oberflächenladung 0,; 

3) die Stetigkeit der tangentiellen Komponenten von &*; 

4) die Stetigkeit der tangentiellen Komponenten von 9*. 

Wenn sich ein Körper im Vakuum bewegt, ist seine Grenzfläche eine Unstetigkeitsfläche 
von v: man wendet dann die Bedingungen 3) und 4) an, indem man annimmt, daß die Ge- 
schwindigkeit auf den beiden Seiten der Grenzfläche dieselbe sei. 

Wenn die folgenden Bedingungen erfüllt sind, a) das Feld in dem Inertialsystem X 
stationär, b) die Geschwindigkeit v auf o tangentiell, oder mit anderen Worten, die Ober- 
fläche des Körpers in sich verschoben ist, hat A. Sommerfeld in Bd. III seiner be- 
rühmten Vorlesungen über theoretische Physik bewiesen, daß die Theoreme 3) und 4) durch 
folgende Theoreme: 

3) Die Stetigkeit der tangentiellen Komponenten von €, 

4’) Die Stetigkeit der tangentiellen Komponenten von 9, ersetzt werden müssen. 

Über das Verhalten von & und $ beim Durchschreiten der Grenzfläche kann man 
überhaupt nichts sagen, wenn die Sommerfeldsche kinematische Bedingung b) nicht 
erfüllt ist, aber folgende Bemerkung kann von Nutzen sein. 

Denken wir uns den Fall, daß sich der Körper im Vakuum bewege und die Oberfläche 
nicht in sich verschoben sei, und nehmen wir an, daß die magnetische Feldstärke ® auf den 
beiden Seiten der Oberfläche gleich sei; dann gilt die Bedingung 3’) auch in dem Fall, daß 
die physikalische Bedingung a) von Sommerfeld nicht erfüllt ist. Dagegen hat die Be- 
dingung 4’) überhaupt keine Gültigkeit. Es wäre natürlich interessant, die Randbedingungen 
systematisch zu erläutern, um die Fortpflanzung der elektromagnetischen Wellen durch be- 
wegte Körper zu erörtern; ich möchte jedoch Ihre Aufmerksamkeit auf ein erstes partikuläres 
Problem lenken: 

Es gebe ein homogenes magnetisches Feld mit der Feldstärke ®, in einem Inertial- 
system K. Wir denken uns einen kugelförmigen isotropen Halbleiter mit dem Radius a, der sich 
mit gleichförmiger Geschwindigkeit v senkrecht gegen die Kraftlinien im Vakuum bewegt. 

Man berechne dann das elektrische Feld und im besonderen die Oberflächenladung. 
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Ein zweites partikuläres Problem ergibt sich, wenn sich die Kugel gleichförmig um die 
Achse dreht, welche dieselbe Richtung wie die magnetische Feldstärke B, hat. 

In dem ersten Problem wird die Sommerfeld sche kinematische Bedingung b) nicht 
erfüllt, sondern es gilt mit hinreichender Genauigkeit meine Bemerkung. 

In dem zweiten Problem wird die Sommerfeldsche kinematische Bedingung 
rerfüllt, doch ergibt sich eine Schwierigkeit, auf die schon Sommerfeld in seinen Vor- 
lesungen hingewiesen hat. 

Sie folgt aus dem Umstand, daß die Minkowskische Theorie für Rotations- 
probleme nicht zuständig ist. Wenn man zum ersten Problem zurückkehrt, so ist die Stetigkeit 
der tangentiellen Komponenten von © gleichbedeutend mit der Stetigkeit des elektrischen 
Potentials y durch o, wenn dieses existiert. 

Zu einem Problem dieser Art gelangen wir, wenn wir das elektrische Feld berechnen 
wollen, das durch die Bewegung der Erde um die Sonne entsteht, sobald man annimmt, daß 
es ein magnetisches Sonnenield gäbe. 

Zweck meiner zweiten Mitteilung ist, die allgemeinen Betrachtungen, die diesem kon- 
kreten Problem vorausgehen, anzuwenden. 

Wir sehen von der Drehung der Erde ab und ziehen nur ihre Translation in Betracht. 


Das elektrische Erdfeld und das magnetische Sonnenfeld; 
ein Versuch ihre Beziehung auf Grund der relativistischen 
Elektrodynamik zu erklären 


Von @. Lampariello in Rom 


$1. Dieexperimentellen Erfahrungen 


Schon vor zwei Jahrhunderten (1752) hat der französische Physiker Lemonnierein 
elektrisches Erdfeld nichtatmosphärischen Ursprungs entdeckt. 

Ich weiß nicht, ob Versuche gemacht wurden, die dieses Phänomen erklären. 

Es ist aber auch bekannt, daß die Sonne ein Magnetfeld besitzt, wie aus der Aufspaltung 
der Spektrallinien (Zeemaneffekt) im Spektrum des Sonnenlichtes geschlossen wurde !) 
und vor kurzer Zeit durch die Untersuchungen von G. Thiessen in Hamburg erforscht 
wurde 2). 

Ich bin der Ansicht, daß das elektrische Erdfeld durch die elektromagnetische Induktion 
des magnetischen Sonnenfeldes bei der Bewegung der Erde um die Sonne erzeugt wird. 

Ich nehme mit Thiessen an, daß die Richtung des magnetischen Sonnenfeldes zur 
Ebene der Ekliptik senkrecht und das Feld homogen ist. 

Ich bezeichne die magnetische Feldstärke mit B;,. 


82. Dieallgemeinen PrinzipienderrelativistischenElektro- 
dynamik und Stellung des Problems 

Denken wir uns zwei Bezugssysteme k und k’, das erste bezieht sich auf die Sonne, das 
zweite auf die Erde. Ich nehme an, daß % und X’ Inertialsysteme seien. Betrachten wir das 
magnetische Sonnenfeld in bezug auf % als statisch. 

In bezug auf k’ nehmen wir ohne weiteres an, daß im Innern der Erde die Ladungsdichte 
o’ und die Stromdichte I’ gleich Null seien. 
Diese Hypothese stimmt mit dem Umstand überein, daß das elektrische Erdfeld in 

bezug auf %’ statisch ist: denn es handelt sich bei X um einen Leiter im Ruhezustand, 

das Verschwinden der Größe o’ und I’ ist die Folge von einer Erweiterung eines durch Experi- 
mente schon lange bekannten Phänomens. 

Hingegen ist anzunehmen, daß die Oberflächenladung o, nicht gleich Null sei; unser 
Problem besteht also in der Berechnung von g.. 

0, I,0,, I,sind Raumladungsdichte, Stromdichte, Oberflächenladung, Oberflächenstrom 
bezogen auf k. 

Ist v die Geschwindigkeit der Erde und 


==, a ar au vr au . Ci), 
so gelten die Formeln ; 
In =a(I—o»v), BEE. N rn» Aaılat), 
Iaıı =4a(l, — 0v), Kri=l, BEE IH DIENT Det (2.2)» 
1) Hale, Astrophys. Journ. 47 (1918), 235. 
R GT hi ess 5 2 Zeitschr. f. Astrophys. 26 (1949), 130; Observatory 69 (1949), 228; Measurements of the 
Sun’s General Magnetic Field, Nature, vol, 169 (1952), p. 147. 
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n— ey; (2.3 
07 =—0 0 8 ® NEE . ); 
age Layg (2.4) 
0, =0|0, — ca DB DR a 5 BE Ei ann s 
und ihre Umkehrungen 5 , 
I,=a(l’+e'v)ı, =D . nl I 
Is =a(ls + 08 ?)]| » Pe] b7, nd, ERBEN ERS BEER BR. (2.2), 
’ 1 4 ' 
o=ale +3 vP) NT ALS vr EM, VEIT Fr (2.3), 
0 ‚ 
0 =a (+5 vL) Be N a 3, 
Nachdem 
De 0, I 0" LE U er) 
ist, folgt aus (2/3), (2/4) daß auch 
dus I=0 13) 
ist. Außerdem ist der Leitungsstrom 
epagyuNn SITES NEE RENREE RER. ER 
Null: 


Hingegen entsteht in bezug auf %k ein Konvektionsstrom auf der Grenzfläche, dessen 
Dichte durch 
== 00% 1.0 00 Kine SE 
gegeben ist. 
Dieser Strom erzeugt, auf Grund des Ro wlandschen Prinzips, ein magnetisches Feld 
B,, das zum bereits bestehenden magnetischen Feld 8, hinzukommt. 
Die erste Gleichung vonMaxwell-Minkowski ist also 


rt. = Be En 
BB, Bl ee vn LE 0 E e, 


wo B durch 


gegeben ist. 

Aber das Magnetfeld, welches von den Ro wlan d strömen auf der Oberfläche der Erde 
herrührt, ist offenbar sehr klein gegen das ursprüngliche induzierende Sonnenfeld und kann 
daher vernachlässigt werden. 

Daraus folgt 

rot & Dip,» ei ee ne ae kn lahs 


es existiert also ein elektrisches Potential y, das wir berechnen wollen. 
Zu diesem Zweck mache ich noch auf folgende relativistischen Formeln aufmerksam 


ol, =o(E+VAB)) NIEREN EST ZEN FE (8.1), 
I, =00(&C+vA DB), arte Ber er (8.2), 


wo co die elektrische Leitfähigkeit bezeichnet. 
Auf Grund von Formel (5), muß 


= + A BER: Fsnulanusee ine mama ie) 
oder 

E=-— vr A Bm AB ulhaitelibeiukiunn ae 
sein und daher, ungefähr 

&E=— v A Bi atahaah Seine a 


Nennen wir y; das elektrische Potential im Innern der Erde. Aus Formel (10) schließen 
wir, daß die elektrische Feldstärke & homogen ist. 

Also ist y; eine lineare Funktion der rechtwinkligen kartesischen Koordinate im Sonnen- 
system k. 

Dieses System Ox yz wird folgendermaßen definiert: Die Achse Oy hat die Richtung der 
Geschwindigkeit der Erde, die Achse Oz die der Feldstärke ®, des magnetischen Sonnenfeldes. 

Aus Formel (10) ergibt sich, daß die innere elektrische Feldstärke die Richtung der 
negativen x hat, 


Ö' 
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Also ist das innere Potential y; durch 


Deoomeı.:.. cken (11) 
gegeben. 
Nun berechnen wir die elektrische Erregung D. Das relativistische Transformations- 
gesetz für Dist durch die Formeln 


2, = (»+5 018), Y=ald+ zuNS) 


gegeben. 
Andererseits hat man 
a or ee) 
und außerdem 
re N ri % Byr: . (14), 


wo e die dielektrische Konstante des isotropen Leiters ist. Daraus a daß in kA’ D’ im 
Innern des Leiters überall gleich Null ist, also in bezug auf k im Innern der Erde die Formel 


1 
gilt. 
Nachdem 
1 
Le er 
besteht, folgt auch 
NS ee ea el). 
Wenn man vom Beitrag der Feldstärke B, absieht, kann man folgende Formel ableiten 
D=— HU N Bd, =HE ee Er eo ERBE. Er (17). 


Und nun gehen wir an die Berechnung des äußeren elektrischen Potentials %,. 

Die Funktion y, muß im Äußeren der Erde und im Unendlichen harmonisch und regulär 
sein. 

Was die Randbedingung auf der Oberfläche des Leiters betrifft wird man behaupten 
können, daß %, gleich y; auf dem Rand ist. 

C ist der Mittelpunkt der Erde mit Geschwindigkeit vlängs der Achse Oy,und so hat man 
Y = vt, wenn %, dessen Abszisse ist. 

Wir führen ein System von polaren Koordinaten ein, dessen Pol C ist. Die Koordinate © 
wird von der Halbgeraden, die von C ausgeht und die Richtung der Achse Ox hat, berechnet. 
Das Innenpotential %; ist durch 


Dee Ban cos Vilnius (11) 
gegeben, und wenn a der Radius der Erde ist, stimmen die Werte von y, auf o 
Ge a ee in ALLE) 


mit denen von y, auf o überein. 
Es kommt also darauf an, y, zu bestimmen, so daß 


Ayı =0 SE a ee ee | 
De com ge. sehn. (19) 
%, Tegulärim Unendlichen ist . - - -» . . . . (20). 


Dieses Problem ist klassisch und wird durch einen Dipol im Mittelpunkt der Erde mit 
einem Moment 
mM—=— vB,a? s 
der die Richtung der Achse Oz hat, gelöst. 
Daher ist das Potential y, durch 


1 ‚a 
Wars 37 =— v B,a® — (22) 
oder durch 
3 
ya vB, ., cos een ‚(22/) 


gegeben. 
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Nun schreiten wir zur Berechnung der Oberflächenladung. Im Innern und Äußeren 


hat man 
D N Ve a a ee re alias 
Daher 
d=—H grad Y;; D, — — &, grad AD nn: Se 16 VERUNERE RER SRRBER (23). 
Die Oberflächenladung e, ist durch 
(Disc + (Don 0m 2 ee ee 
oder durch 
day; N En € ’ 
— ee rndne Ma Silweraeee 
86 & B dn, ( ) 
gegeben. 
Durch die Formeln (11), (227) hat man, wenn man 
Eu ER 
dn; Or dn, Or 
in Betracht zieht, 
dp; My; 98 
FR 3} —=—vB,cos 9 (25), 
dpa My > a? 5 
== 9 ie 26 
In er vB; cos © (26) 
Also, auf o 
dp; | 
I ——— —H—rg 5 ee ee 27 . 
Fr en By 3v»B,cos 9 (27) 
Und als Endergebnis hat man 
D-30860 0: a a ee ne (28) 
und auch 
&=3 VW cos 9 35 = Hn008 0...  uecefsp.nce (28)). 
Das Maximum von g, ist also 
v 
Oamav = 3 5 Ho anst ratelen ar dalhn e lah nel. 


Um die Oberflächenladung 0, in bezug auf %’ zu berechnen, muß man die Formel (2.4) 
berücksichtigen, in der man 


I, — 0% v 
an Stelle von /, den Wert 0,v setzen muß. 
Denn da 
&— zL=o- Po —=a(l—P) Br) DB, Di, AROTB 
ist, hat man nach (2,4) 
ol, B%)0, ih Bio, ae ei 
Also kann man annähernd annehmen, daß 
De ee en 2 


ist. 
Der Wert von H,ist Schwankungen von 0 bis 50 Oersted unterworfen, so daß das Maxi- 
mum von 0, durch 
omas = 9. 10-1 Coulomp/ Me ee a ME 
gegeben ist. 
Die erhaltene Formel für Ogyar 


v U h 
Coma = 3 a Bo a 


kann umgekehrt dazu dienen, die Erregung des magnetischen Sonnenfeldes zu berechnen, wenn 
die Dichte 0, auf experimentellem Wege bekannt würde. 
Man kann die Lösung des Problems der Rotation durch eine analoge Methode erhalten. 
Um das Innenpotential zu bestimmen, muß man den klassischen Ausdruck der Geschwin- 
digkeit v, der durch die Formel 


gegeben ist, an Stelle von v setzen. 


Ö 
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Daraus ergibt sich 
Era EN ne RE rn 
Diese Größe wird vom Potential 
y; =} B,o r?sin? © 
abgeleitet. 
Dann muß man ein erstes Randwertproblem lösen, um das Außenpotential y, zu berech- 
nen. Die Randwerte des Potentials sind durch die Funktion 
(9) =} B,wa? sin? © 


gegeben. Es ist seit langem bekannt, daß die Lösung durch das berühmte Poissonsche 
Integral ausgedrückt werden kann. 


Eine Methode zur Berechnung des charakteristischen Exponenten 
einer Hillschen Differentialgleichung 


Von Friedrich Wilhelm Schäfke in Mainz 


Es wird ein Verfahren angegeben, das es gestattet, den charakteristischen Exponenten 
einer Hillschen Differentialgleichung y'’(t) + (A + f(t)) y(t) =0 zu jedem Werte des Para- 
meters Ain numerisch brauchbarer Weise direkt zu berechnen, wenn nur etwa die Eigenwerte 
»=,(n =0,1,2,--:) der Periodizität bekannt sind. Zur Vermeidung von Fallunterschei- 
dungen, die für den Gedanken der Methode unerheblich sind, wird die Mathieu sche Diffe- 
rentialgleichung 

GE ED DYNO .E ur. 20 0. ll) 


mit einem reellen h? #0 betrachtet. 
Nach dem Thorem von Flo quet gibt es immer (mindestens) eine nichttriviale Lösung 
mit 
ya) = er y(). 
» wird als charakteristischer Exponent bezeichnet; mit vsind alle Zahlen 2 k + » mit ganzem k 
mögliche charakteristische Exponenten. Der Zusammenhang mit A wird durch eine Gleichung 


COSZE VEmET(A) Tee ee Eau von un. (2) 
vermittelt (h? wird als konstant betrachtet). f(A) ist eine ganze analytische Funktion von } und 
als solche höchstens von der Wachstumsordnung 1/2 }). 


Mit f(A) ist auch f(A)— cos zz», bei beliebigem », höchstens von der Wachstumsordnung 1/2. 
Es kann daher die Weierstraßsche Produktdarstellung verwandt werden. Sie liefert 


cos ay—cosnwm +8 A—Ant2r (A?) el 
),— An +2r (A?) 
Dabei sind A, +2r(h?)(r =0, +1, + 2, - -)sämtliche Wurzeln von 

Ja) -cosrm=0, 


ihrer Vielfachheit entsprechend gezählt; sie sind (bei reellem A?#0) für reelles », reell und 
verschieden. Außer diesem System der A, +2r (h?) muß zu einem A, 


COS TY — COSTYy 


v=—o 


A) =eosrmy#cosn 


bekannt sein. Dann liefert (3) zu jedem Aden zugehörigen Wert von cos zz v. Numerisch ist diese 
Formel unbrauchbar; die Konvergenz ist zu schlecht, da die Faktoren nur die Größenordnung 
1 +0(r”2) haben. 

Eine wesentliche Verbesserung erreicht man bei Berücksichtigung des asymptotischen 
Verhaltens 


Ay 2r(A?) —— (v9 + 2 r)2 + Ol): ce I: Ge ar Er er (4), 
Man kann dann leicht zeigen: 
COS TV — COST +® A — A,+2r(h?) ui Au +2r(h?) — (un + Zr)? E 
ie be = J |!- In B). 
Era co run wat  MmFrr" „Zn A— (w+2r) 


(5) bietet zwei Vorzüge gegenüber (3). Einmal ist die Konvergenz verbessert; die Faktoren sind 
1 +0 (r). Zum anderen wird kein zusätzliches Paar »,, A, benötigt. 


1)F.W.Schäfke, Math. Nachr. 3, 20—39 (1949). 
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Ist noch ein solches Paar bekannt, so kann die Konvergenz weiter verbessert werden: 


(A, «> 4) (Ant 2 r(h?) I) (% +2 r)?) (6) 
(A, — (+ 2r)?) ((A — Ay+2r(h?)) 2 


COS TYL—COSNY  COSMYA, — COS 7% in 
coS TV — COS TY cos ii COS TV, r=—» 
Die Faktoren sind 1 +0 (r*). 

Eine zweite Möglichkeit zur Verbesserung der Konvergenz besteht im Übergang zu Par- 


tialbruchformeln mit Hilfe logarithmischer Differentiation von (5). 
Für », =0 und », =1 erhält man insbesondere 


1 + 


Pr => ef S 1 d, L IT t ya 
sur Ti) A—a| ap 2° 
(N), 
zen 5 4 5 1 1 a ey 
cart TATEN | Fayı = 
woraus man 
Be d OHR 
cosnv Se ST 7 °05 7» Sy DL, (8) 


berechnet. Hier zieht man die ganz- und halbperiodischen Eigenwerte heran, die ja im Falle 
der betrachteten Mat hie uschen Differentialgleichung weitgehend tabuliert sind, und erhält 
dafür neben cos x» noch cos av. Die Reihenglieder sind O (r®). 

Wegen der Beweise sei auf meine Arbeit Math. Nachr. 4, 175—183 (1950) und das beim 
Springer-Verlag in Vorbereitung befindliche Werk J. Meixner-F.W.Schäfke,Ma- 
thieusche Funktionen und Sphäroidfunktionen mit Anwendungen auf physikalische und 
technische Probleme, verwiesen. 


Algebraisierung der endlichen Hilbert- Transformation 
Von H. Söhngen in Darmstadt 
In den Anwendungen tritt öfters die endliche Hilbert - Transformation 


rt i1<acı) 


auf, wobei das Integral an der singulären Stelle im Sinne des Cauch yschen Hauptwertes 
zu nehmen ist. Es erscheint daher nützlich nach einer Transformation T zu suchen, die die 
endliche H-Transformation in eine Multiplikation übersetzt, ähnlich, wie die Laplace- 
"Transormation die Operation der Differentiation in eine Multiplikation verwandelt. 
Betrachtet man die Transformation 
1 


za} = re 


so wird dabeieine Funktion F(&) der Klasse L’(p > 1) auf eine in dem Streifen |Rs| <1 — 1/p 
reguläre analytische Funktion f(s) abgebildet. Diese Transformation kann umgekehrt werden 
und man erhält (aus der Verwandtschaft mit der zweiseitigen Laplace-Transformation, 
e —=(1 + «)/(1 — x)) unter geeigneten Voraussetzungen über F(x) die Umkehrformel 


ıT 
1 Fa 1 1 za n 
rer a lee: f(s) ds a a er Yarte ( j* 
Weiter gilt 


{HFC} = eig me ey All FE 


sin Ts 


Als Anwendungsbeispiel soll die Bodenpressung P einer Staumauer bestimmt werden. 
Hierfür hat Tölke das Integralgleichungssystem 


VY()=— eH{H(a)}, H(a)=ck +.9{V lo} -» - -. (4) 
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angegeben, wobei &e>2 durch diePoisson sche -Konstante des Mauerwerks gegeben ist und 
k ebenfalls eine Konstante ist. Gesucht sind Hund V. Setzt man F(x) = H (x) + iV (x), so 
geht das System (4) in die Integralgleichung 


EEE ee 2 5) 
und diese nach Anwendung der T-Transformation in 
re N F(0) 
Aalzs en TTS ie |eotg ei er sin zzs 


über. Mit &e =cotg zg wird dann 
ae llz 2 sk) 
cosr(s + ig) 
und die Umkehrformel (2) liefert als allgemeine Lösung 
des Integralgleichungssystems (4) 


1+ 20% : 
F(x) = = | ) 0) — 27k ink], 
wobei f (0) beliebig vorgegeben werden kann. 
Die Integralgleichung (5) ist ein Spezialfall von 
F(a)+ 2H{F@)} -G (a), 
wobei @(x) einer Klasse LP(p > 1) angehört und A beliebig komplex ist. Gesucht ist die all- 


gemeine Lösung F(x) aus der Klasse Z!t°(e>0). Die oben dargelegte Methode liefert dann, 
falls A nicht in dem Intervall von —: bis +t liegt, 


Bild 1 


+1 
Ze: 21 0me ge) 
un nee) - ar I er. :® 
1 1 
+— 2 100) cos au: 1 a, 
wobei 
ey > | 
rer mis - TR <0 


und f (0) beliebig ist. 
Ganz entsprechend lassen sich natürlich auch Integralgleichungen behandeln, in denen 
der Operator 9 mehrfach auftritt. 


Ein ausführliches Referat des Vortrages: 
E.A. Deuker, Über die Verteilungsfunktionen von Vektorsummen 
wird in einem der nächsten Hefte dieser Zeitschrift erscheinen. 


Beferate der Vorträge: 


P.Mönnig, Über die Lösung der Hamilton - Jacobischen Differential- 
gleichung durch Trennung der Variablen, 


K. Stange, Über die Gesetze der Kornverteilung bei Zerkleinerungsvorgängen 
sind nicht eingegangen. 


C. Dynamik und Stabilität 


Beschränkte Funktionen, Frequenzcharakteristiken elektrischer 
Netzwerke und alsebraische Stabilitätskriterien 


Von H. F. Effertz in Aachen 


Ein elektrisches Netz wird bekanntlich durch das folgende Differentialgleichungssystem 
mit konstanten Koeffizienten beschrieben: 


Z 30 Dri,= Du Beim). Sen: nei (l) 


s=1 7=0 
D bedeutet den Differentialoperator nach der Zeit. Nimmt man die Spannungen in der Form 
u, — U,e* an, so gelangt man mit Hilfe des Ansatzes ti, = I,(2) e* bei Auflösbarkeit zum 
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Gleichungssystem 


L1Q)= fr) U, REST 


Die Frequenzcharakteristiken f,, (2) sind gebrochen rationale Funktionen des komplexen 
Frequenzparameters 2. Der Zusammenhang der Koeffizienten der Frequenzcharakteristiken 
mit dem topologischen Aufbau des Netzwerkes und seinen Schaltelementen kann unter Be- 
nutzung der topologischen Verknüpfungszahlen in Determinantenform angegeben werden. 
Setztman U,= U, = :::—= U„=0, so gehen die beiden ersten Zeilen von (2) in die Vier- 
polgleichungen über. Die Betriebseigenschaften des Vierpols werden durch die Frequenz- 
charakteristiken vollständig bestimmt. Das technische Problem ist aber umgekehrt: Die 
Betriebseigenschaften sind vorgegeben und die zugehörige Schaltung wird gesucht. Ausgangs- 
punkt der sich hier anschließenden Äquivalenz-, Interpolations- und Approximationsprobleme 
sind die funktionentheoretischen Eigenschaften der Frequenzcharakteristiken. Dazu hat 
W.CGauer den Begriff der positiven Funktion eingeführt. Danach heißt eine Funktion f (2) 
positiv, wenn sie in der rechten Halbebene regulär, auf der reellen Achse reell ist und in der 
rechten Halbebene positiven Realteil besitzt. Den positiven Funktionen sind wegen ihrer 
technischen Bedeutung zahlreiche Untersuchungen, insbesondere von W. Cauer, gewidmet 
worden. Aus den klassischen Arbeiten über beschränkte analytische Funktionen von C. 
Carath&odory und J. Schur können weitergehende Resultate gewonnen werden. 
Zum Beispiel erhält man zwei charakteristische Darstellungen für rationale Funktionen 9 (z), 
die in der rechten Halbebene regulär, für imaginäre z (bis auf Pole) imaginär sind und in der 
rechten Halbebene positiven Realteil besitzen. Es gilt die Partialbruchdarstellung 


k 
j (, — 1)z+(1+ 8) 3 
p(2) ee ii . (3), 
wobei breell, r„>0Oist und alle e, voneinander verschieden sind und |e,| = 1gilt. Verlangen 
wir, daß 9(z) für reelle z reell sei, so geht (3) in die bekannte Partialbruchdarstellung von 
R.M. Foster!) für Frequenzcharakteristiken von Netzen ohne Oh msche Widerstände, die 
Reaktanzfunktionen, über. Den Reaktanzfunktionen entsprechen geometrisch und analytisch 
leicht zu charakterisierende Klassen der vonC.Carath&odoryund J.Schur betrachte- 
ten Potenzreihen. Die zweite Darstellung hat die Form 


fa ef 2) EEE A 
Se) + ef 2) 
wobei f(z) Hurwitzpolynom k-ten Grades, f(z) das Polynom mit konjugiert komplexen Koeffi- 
zienten und |e| = list. Verlangen wir hier ebenfalls, daß (z) für reelle z reell sei, so erhalten 
wir das Reaktanztheorem von W.Cauer?). Mit Hilfe von (4) gewinnt man leicht folgendes 
Theorem: 

Gilt %B>0, el, 50 ist, PR) 2020, 
stets und nur dann Hur wit zgleichung, wenn 


P@)= alP() — eP(— Q]+ BLP@)+eP(- 2)] } 
Hurwitzgleichungist. Man kann dieses Theorem auch aus einem Satz von J. Sch ur?) her- 
leiten. Zur Gradreduzierung von Hurwitz polynomen haben J. Schur ®)und H. Bückner %) 
zwei verschiedene Verfahren angegeben. Die beiden Methoden lassen eine elektrotechnische 
Deutung zu, die ihren Zusammenhang aufklärt. Ist 


I) = (2? + 0,2272 4. + 09n) + (a 22° 1 + 0,2273 +» + Qgn-ı2) 
Hurwitzpolynom (2n > 2), so ist nach dem Reaktanztheorem 
2n Se 
F(e) = Io?" este 


p(2) = 


a, zen iL 4 dan 
Frequenzcharakteristik der Kettenbruchschaltung: 
Zo Zz Zzn-2 


Bild 1 


1) R.M.Foster, Bell Syst. Techn. J. 3 (1924), S. 259. 
2) W.Cauer, Elektr. Nachr. Techn. 16 (1939), S. 96. 
3) J.Schur,Z.angew. Matb. Mech. 1 (1921), S. 307. 
#) H.Bückner, Quart. appl. math. 10 (1952), S. 205. 
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Die Größen der Induktivitäten bzw. Kapazitäten Z, ergeben sich nach der Bader- 
schen °) Kettenbruchformel Z,— R, Ry4ı, =0,:--,2n—1). Die R, sind hierbei die 
Routhschen Probefunktionen des Polynoms f(2). Baut man die Schaltung dadurch ab, 
daß man nacheinander Z, Z,***, Zen-2 Null setzt, so erhält man entsprechend 
die Frequenzcharakteristiken Fs, Fas, Fss '°*", Fien-ns. Bildet man die Summen von 
Zähler und Nenner dieser Funktionen, so erhält man die reduzierten Hur wit zpolynome nach 
Schur fis fs» **"» fen-ys. Baut man die Schaltung — jetzt hinten beginnend — ab, 
indem man nacheinander Zen), Zan-a, '**, Z, unendlich groß werden läßt, so erhält man in 
analoger Weise über die Frequenzcharakteristiken Fr, Fa5,°**, Fen-yz die reduzierten 
Hurwitzpolynome nach Bückner fi», fer, ***, fien-yz. Bei ungeradem Grad der 
Ausgangsgleichung ändert sich das Verfahren nur unwesentlich. 

Die Variabilitätsbereiche der Koeffizienten Hur witzscher Polynome und ein darauf 
gegründetes Verfahren zum Aufbau stabiler Frequenzfunktionen wurde mitgeteilt. Aus den 
Werten der Hur witz determinanten können nicht nur Aussagen über die am meisten zu 
befürchtende Eigenfrequenz ®), sondern auch Aussagen über die Stabilitätsgüte gewonnen 
werden. Dazu wurde eine allgemeine Darstellung der quadratischen Regelfläche als Funktion 
der Anfangswerte und der Hur wit z determinanten angegeben. 


Einige Stabilitätskriterien 
Von Friedrich W. Schäfke in Mainz 


Für eine Hillsche Differentialgleichung zweiter Ordnung in Normalform werden zwei 
Scharen von Kriterien angegeben, die die Stabilität sichern, falls bei gegebener unterer (oberer) 
Schranke für die Koeffizientenfunktion deren Mittelwert nicht zu groß (klein) ist. 

Es wird die Differentialgleichung 


RES IPN LIE ur  ER D| 
betrachtet, in der f(x) für — © <x < o stetig, m-periodisch und nicht konstant ist. 
Wir zeigen: 


KııterrtumT? sei 


IR) zen Zar ln IP TREO IHRER RE UR (a) 
und 
iur ER: 
— | fa) da < a8 4 IE ren Je na (b); 


0 


Dann sind die charakteristischen Exponenten von (1) verschieden und stabil?). 
Beweis: Wir führen den Beweis für x <0. Er verläuft für « =0 analog und liefert 
dann das bekannte Kriterium von Liapounoff. 
Wir nehmen an, unsere Behauptung sei nicht erfüllt. Dann gibt es eine reelle Lösung 
y(z) = 0, deren Nullstellen sich z-periodisch wiederholen. Wegen (a) muß dann gelten: 


vn) =0, ya) =0 
für 
SEE RUE N << H<mHrSsm tr; 
auf Grund des Sturmschen Vergleichssatzes ist die Länge jedes der 2n + 2 Teilintervalle 


nicht größer als 37% 


Betrachten’wir das erste Teilintervall und wählen 
n(x) — [sin (2 — %,) 


mit 
n"+0n=0 


5) W.Bader, Arch. Elektrotechn. 34 (1940), S. 293. 

6) L.Cremer,Z.angew. Math. Mech. 25/27 (1947), S. 161. 

1) Für & = 0 ist rechts der Grenzwert 4/r°? zu wählen. j - h 

2) D.h.: Es gibt zwei nichttriviale Lösungen y,, ymit yi(l+ x) = e+tirr yi(t), wobei» reell und nicht 
ganz ist. 
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als Vergleichsfunktion, so erhalten wir in bekannter Schlußweise 
&o 
/ G & 
"9 Yna=/U-)ynde, 
und eine einfache Abschätzung ergibt 
&o 
actg &(&, — ©) < | (fe) — 0°) de; 
To 


dabei kann =nur stehen, falls dort f(x) = 0°. 
Analoges gilt für die übrigen Teilintervalle. Daher folgt: 


m In+1 
1 If 
ai fiade— az | U) a) da 
2 1 
> 2 fg al) + eig) + + eg ee 


im Widerspruch zu (b). 
Bemerkung 1: Die Beweisführung läßt erkennen, daß es zu jedem e> O0 stetige, 
rr-periodische, nicht konstante Funktionen f(x) mit (a) und 


I | fa) da = 


0 


21 2(n +1) Am 


gibt, derart, daß (1) eine Lösung y(x) mit z-periodisch wiederkehrenden O-Stellen, also nicht 
verschiedene und stabile charakteristische Exponenten hat. Unsere Schar von Kriterien kann 
nicht durch Vergrößerung der Schranke für den Mittelwert verbessert werden. 


Bemerkung 2: Die rechte Seite von (b) ist größer als «(n-+1). Dies liefert ein 
schwächeres, jedoch praktisch oft handlicheres Kriterium. 


Bemerkung 3:Füra =1,2,3,4,5sind wesentlich schwächere Kriterien dieser Art 
von R. S. Gusarova?°) angegeben worden. 
Wir beweisen weiter : 
Kriterium 2: Sei 
fe) 20%, M<esh +17 a0, 12270250) 
und 


1 7 
2 | Ha) ran nen sand) 
() 
Dann sind die charakteristischen Exponenten von (1) verschieden und stabil. 
Beweis: Wir betrachten zuerst den Fall n —0, wo die untere Schranke von (d)O und 
von & unabhängig ist. Wäre hier unsere Behauptung falsch, so müßte es wegen (d) eine Lösung 
y(z)==0 von (1) geben, deren O-Stellen z-periodisch wiederkehren?); dem widerspricht jedoch 
(c) auf Grund des Sturmschen Vergleichssatzes. 


Seinun n> 0, wieder die Behauptung nicht erfüllt und y(z) = 0 die reelle Lösung von 
(1) mit z-periodisch wiederkehrenden O-Stellen. Sind dann 


< %< D< ... 
aufeinanderfolgende O-Stellen, so gilt wegen (c) 


BER OD 
2: — %_ a 
D v1 & 


und =nur für f(x) = a®im betreffenden Teilintervall; 


3) Priklad. Math. Mech. Moskau 1949, 241—246; Priklad. Mat. Mech. Moskau 1950, 313—314. 
4) Anderenfalls würde es eine reelle Lösung y(z) == 0 mit z-periodischem ee geben und man hätte daher 
{ y\ıX 


[In de 
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daher ist 
nı<orT 
und mithin 
4. = m. tr 
mit einem 
k<n. 


Betrachten wir das erste Teilintervall und vergleichen mit 
I 
1,(&) = sin (zer) W=0,1), 


so wird, wie beim obigen Beweise, 


+, + Er 
| (@—-f)ynd= ym—ny , 
To 0 
und mithin 
Ber R at 
( — f)de > I 
Y 
To 
wobei = nur steht, falls dort a® = f. Analog folgt 
2, Wr 
i ABER 
= ie 2 — — 
| (ar er). Y\ 
rin 
Nun gilt 
er ET ZT” 234 
‚ EL N 
= a y e I ie — | 2) da + | Fe) da. 
der er 24 5 | 204 37 


Daher ergeben die drei vorstehenden Zeilen 


fe — f(x) da 2.0? (a N —&) 


& 


mit =nurin dem Falle, daß «? = f(x) in [%,, %,]. Entsprechendes gilt für die übrigen Intervalle; 
also hat man 
. n 77. 
je je = )j2=l-2) 
Be uk 2: e- Bi 1—- —)2 21 —). 
z (x? — f(a)) da = (x — f(x)) de > a —)2 % . 
To 

Das widerspricht aber der Voraussetzung (d). 

Bemerkung4:Esgilt Ähnliches wie in Bemerkung 1. 

Bemerkung 5: Kriterium 2 stellt für n—=0 eine Verbesserung eines Kriteriums von 
Gusarova)dar, der statt (c) stärker fordert f(x) + |min f(@)| <1, 
ferner des Kriteriums von S. Wallach), der statt (d) f(x) z O annimmt. 

Bemerkung 6: Aus beiden angegebenen Kriterien folgt das allgemeine Kriterium 
von S. Wallach’), das Stabilität für 

n?<f(a2) s(n+1)? (n —=0, 12,..+.) 
behauptet. Es ist dann nämlich entweder 
7 
1 
I | ya) da < nn +1) 
ö 

und Kriterium 1 anwendbar oder 


1 [so de n(n-+1) 


und Kriterium 2 anwendbar. 


5) Die erste Arbeit in Anm. 3. 
6) Proc. Nat. Acad. Sei. USA 34, 203—204 (1948). 
?) Vgl. Anm. 6. 
21 
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Vereinfachung der Grundsleichungen der elastischen Stabilität 
mit Anwendung auf Stäbe, Platten und Schalen 


Von H. Neuber in Dresden 
1. Übersicht 


In einer früheren Arbeit [1] hat der Verfasser nachgewiesen, daß bei Aufstellung der 
Grundgleichungen der elastischen Stabilität unter Anwendung der exakten Rechenverfahren 
des Ricci- Kalküls innerhalb der ‚‚linearen‘‘ Stabilitätstheorie insofern eine Unsicherheit 
besteht, als die Veränderung des Vorspannungstensors während der bei Instabilität eintretenden 
Deformation genaue Aussagen über Größen höherer Ordnung sowohl bei der Formulierung des 
Elastizitätsgesetzes, als auch der bis zum indifferenten Gleichgewicht erfolgten Vordeformation 
erforderlich machen würde [2]. Die diesbezüglichen Größen sind jedoch so klein, daß sie sich 
den Meßmethoden entziehen. Andererseits liegen diese Effekte meist auch innerhalb der Her- 
stellungsgenauigkeit der jeweils zu untersuchenden Bauteile. Aus diesem Grunde konnte vom 
ingenieurmäßigen Standpunkt aus von vornherein erwartet werden, daß sie keinen wesent- 
lichen Einfluß für die kritische Last zur Folge haben. Diese Vermutung läßt sich für das Bei- 
spiel der Knickung prismatischer Stäbe bei reiner Druckbeanspruchung bestätigen. Hieraus 
ergibt sich für die Theorie der elastischen Stabilität folgender Sachverhalt: Die Veränderung 
des Vorspannungstensors während der bei der Instabilität möglichen Deformation ist bei einer 
Reihe technisch wichtiger Stabilitätsprobleme ohne wesentlichen Einfluß auf die Größe der 
kritischen Last. Die mit der Veränderung des Vorspannungstensors in Zusammenhang stehen- 
den kleinen Größen beteiligen sich jedoch an der Formulierung der für indifferentes Gleichge- 
wicht geltenden linearen Differentialgleichungen und sind für den Schwierigkeitsgrad der zur 
Integration erforderlichen Rechnungsgänge von ausschlaggebender Bedeutung. Daraus ergibt 
sich die lohnende Aufgabe, innerhalb der vom Standpunkt der Werkstoffmechanik aus zu er- 
wartenden Möglichkeiten der Formulierung eine Darstellung zu finden, welche für eine mög- 
lichst einfache Integration am vorteilhaftesten ist. In der erwähnten Arbeit des Verfassers 
wurden bereits zwei derartige Ansätze diskutiert. In einer weiteren Arbeit des Verfassers über 
das Stabilitätsverhalten von Sandwich platten [3] wird ein weiterer Ansatz verwendet, 
der gewissen naheliegenden Vorstellungen entgegenkommt. 

In vorliegender Arbeit wird ein neuer Ansatz für den Vorspannungstensor angegeben, 
welcher die Integration wesentlich vereinfacht. wobei sich die mit den Komponenten des Vor- 
spannungstensors behafteten Glieder innerhalb der linearen Stabilitätstheorie gegenseitig Kom- 
pensieren, so daßsich die Gleichgewichtsbedingungen auf das Verschwin- 
den der Divergenz des durch die indifferente Deformation geweck- 
ten (bzw. mit derselben durch das Elastizitätsgesetz kausal verknüpften) Zusatzspan- 
nungstensors reduzieren. Sie gleichen daher formal den Differentialglei- 
chungen der gewöhnlichen Elastizitätstheorie bzw. den Gleichge- 
wichtsbedingungen des undeformierten Körpers. Damit stehen aber zur 
Lösung die Verfahren der gewöhnlichen Elastizitätstheorie zur Verfügung, insbesondere 
im allgemeinen Fall des dreidimensionalen Problems der vom Verfasser eingeführte 
„Dreifunktionenansatz“ [4]. Wie nachstehend gezeigt wird, tritt der Vorspan- 
nungstensor dann lediglich noch bei der Aufstellung der Randbedingungen in Erschei- 
nung. Die Brauchbarkeit der gewonnenen Methode wird zunächst am Beispiel der 
Knickung druckbeanspruchter prismatischer Stäbe nachgewiesen (Bestätigung der Euler- 
last). Ferner wird für den Fall der Plattenknickung gezeigt, daß sich die allgemeine 
Differentialgleichung dieser Problemgruppe mit Hilfe der neuen Formulierung unmittel- 
bar aus der gewöhnlichen Plattentheorie ergibt, wobei sich Übereinstimmung zeigt mit bis- 
herigen Formulierungen, insbesondere den Arbeiten von K. Marguerre [5] (Anwendung 
der energetischen Methode) und Fr. A. Willers [6] (Formulierung als Variationsproblem 
mit zahlreichen Anwendungen). Ebenso lassen sich die verwickelten Differentialgleichungen 
der Schalenknickung in verhältnismäßig einfacher Weise aus der tensoriellen Schalentheorie 
[7] des Verfassers herleiten, worauf an anderer Stelle noch näher eingegangen werden soll. 


2. Darstellung der Grundgleichungen der elastischen Stabilität mit Hilfe des Ricci- 
Kalküls und Aufstellung eines vereinfachten Ansatzes für den Vorspannungstensor 

Mit den kartesischen Koordinaten = %; y = 3; 2 =7, allgemein %,, %, &, %, mit 
k, l,m, n=1, 2. 3) und den kartesischen Komponenten V,, V, der V,, des elastischen Ver- 
schiebungsvektors gilt für die Koordinaten nach der Deformation 


ge = ze 4 yrl „We AB ABI BE EI 
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Werden zur Abkürzung dieselben Bezeichnungen verwendet, die vom Verfassef in den unter 
Anm. [1]und [7]angeführten Arbeiten benutzt wurden, so gilt für das Gleichgewicht des Span- 
nungstensors nach der Deformation 


a Ru AR nun. (2). 


Der Spannungstensor nach der Deformation K?“ besteht aus dem Vorspannungstensor, der 
infolge der kleinen Formänderungen des Werkstoffes, wie sie bei Formulierung des indifferenten 
Gleichgewichtes (entsprechend der „linearen“ Stabilitätstheorie) in Betracht gezogen werden 
müssen, vom ursprünglichen Vorspannungstensor 7?“ in den etwas veränderten Tensor Ta 
übergegangen ist, sowie einem kleinen Zusatztensor, der mit diesen kleinen Formänderungs- 
größen durch das Ho o ke sche Elastizitätsgesetz kausal verknüpft ist. Wird der Zusatztensor 
mit 1?# bezeichnet, so gilt: 


Be ra); 


Für die exakte Herleitung von 7?“ wäre eine Erweiterung des Formänderungsgeselzes auf 
Größen höherer Ordnung erforderlich, die sich den Methoden der Meßtechnik infolge des Hin- 
zukommens verschiedener Fehlerquellen entziehen. Deshalb wurden in der unter [1] ge- 
nannten Arbeit des Verfassers bereits zwei heuristische Darstellungen angegeben. Eine weitere 
Darstellung, welche gewissen Annahmen über die Unabhängigkeit von der starren Drehung 
des Werkstoffteilchens entspricht, findet sich in der unter [3] angegebenen Arbeit des Ver- 
fassers. Die nähere Untersuchung dieser Ansätze zeigt am Beispiel des gewöhnlichen Knick- 
stabes, daß die kritische Last stets mit der Euler last übereinstimmt, d. h. von der speziellen 
Art der Ansätze unabhängig ist. Der zugehörige Nachweis macht jedoch bei den erwähnten 
Ansätzen die Durchführung der Integration der erweiterten Gleichgewichtsbedingungen erfor- 
derlich, wobei der Vorspannungstensor als Parameter auftritt und dadurch die Rechnung 
wesentlich erschwert. 

Da das Endergebnis offenbar innerhalb gewisser Grenzen von der Art des Ansatzes über 
den veränderten Vorspannungstensor unabhängig ist, soll hier ein neuer Ansatz aufgestellt wer- 
den, der für die Integration besondere Vorteile mit sich bringt. Zunächst seiein Restspannungs- 
tensor »*# eingeführt, welcher der Veränderung des Vorspannungstensors Rechnung trägt, so- 
weit diese nicht unmittelbar mit dem Elastizitätsgesetz in Zusammenhang steht. 

Dann gilt allgemein: 

Tree ARE Er uR ee Ev Em Tape (A) 
wobei r*# klein gegenüber 7** ist, und zwar von derselben Ordnung wie i**. Aus Gl.(3) folgt 
nunmehr: 

VIE Sl Fr 
Setzt man diesen Ausdruck in die Gleichgewichtsaussage (2) ein, so folgt mit Rücksicht auf die 
Tatsache, daß der Vorspannungstensor im undeformierten Körper im Gleichgewicht war, d.h. 
die Bedingung 


YES rn Al Di }) 1 Be Ve (6) 
erfüllt, unter Beschränkung auf Glieder erster Ordnung: 
Te ya Te 0. 2.2.2. lM. 


Wie eine kurze Rechnung bestätigt, bedeutet diese Gleichung das Verschwinden der Diver- 


genz eines Hilfstensors #: 
tn ae (8) 


wobei 

AR RaVa, Tar L Vu, Tr — VTıe, 2.22... 
wird. 
Mit Bezug auf die anfangs dargelegten Gedankengänge, kann über den Tensor r*# inner- 
halb gewisser Grenzen frei verfügt werden, ohne die Endergebnisse zu beeinflussen. Derjenige 
Ansatz verdient hierbei besonderen Vorzug, der eine Vereinfachung der Rechenmethode herbei- 
führt. Der einfachste diesbezügliche Ansatz ist offenbar der, für den 


Be en... nee. . (10) 
d.h. 
Aa — VA, Ter _ ya, Tr {vr Teu,.....2000.. (ll) 
wird, denn in diesem Falle liefert Gl. (8): 
ur LE Br Er EEE Er E77) 


91% 


-_ 
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also eine Bedingung, welche mit der Gleichgewichtsbetrachtung des undeformierten 
Körpersim Einklang steht. Da der Zusatztensor durch das H o o k e sche Elastizitätsgesetz 
mit den Ableitungen des Verschiebungsvektors linear verknüpft ist: 


| > 


gilt für das so entstehende Integrationsproblem die gewöhnliche Elastizitätstheorie und damit 
der vom Verfasser aufgestellte Dreifunktionenansatzs [4]. Die jeweiligen Randbedingungen 
sind allerdings für den Tensor K?* zu formulieren, bzw. da der Vorspannungstensor vor der 
Deformation 7*# bereits die Randbedingung erfüllt, für die Tensorsumme Ar rin, 


2: 
Ws ENT) 7 ulA ı Vr 
1 - (m Ha gie 


3. Nachweis der Methode am Beispiel der Knickung des druckbeanspruchten 
prismatischen Stabes 

Um die Brauchbarkeit des Ansatzes nachzuweisen, sei als Beispiel der druckbeanspruchte 
prismatische Stab betrachtet. Wird angenommen, daß es sich um eine dünne Platte handelt, 
so kann ebener Formänderungszustand vorausgesetzt werden, und der Dreifunktionenansatz 
reduziert sich infolge der Abhängigkeit von nur zwei Koordinaten auf zwei Funktionen. Wird 
gemäß Bild 1 die x-Koordinate in Richtung der Plattenmittelfläche (zugleich Druckrichtung) 
angenommen, so gilt für den Dreifunktionensatz [4] das Gleichungssystem (mit ®, = 0; 
8,0): 


Lg File F,, eat ( ii; = AF 


lay = — Piss +0 D,,« 
A a 1- >) ar r2080,, ne 
26V, =—-F: 
26V, =—F,+20%, 
v0 
Hierbei ist: 
F=6 20, m AD 0, AD ee DT ee . (15) 


Entsprechend der zu erwartenden Knickdeformation sei gesetzt: 


&=0 Sin” sin ” 


D,— 007 Co)” sin” 1 BLEI A MORE ARSTER 


I l 
Ne o(ein"! +0" 60” Yin | 


"Nach Einsetzen in das Gleichungssystem (14) ergibt sich: 


2 0% I H2—o)e] Mer: sin 

- or ı a8 Co) 185 0 ein“ cos" 

er O1 + ac] Ein? —o Yo si sin er 
DAY ERO a Sin”? I 08) cos" 
207, = (7 u + (2% —1) ce] CT —o "Vein 2) sin®? 


Durch den harmonischen Charakter der Ausgangsfunktionen sind gemäß dem Drei- 
funktionenansatz [4] sämtliche Elastizitätsgleichungen (Gleichgewichts- und Kompatibilitäts- 
bedingungen) von vornherein erfüllt, so daß nur noch die Randbedingungen in Betracht zu 
ziehen sind. Diese lauten: 


trat; bay + ray = Se 
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Entsprechend (Gl. (11) ergibt sich hieraus: 
vl, A 0 Se =(19) 


Hierbei wurde beachtet, daß entsprechend den eingeleiteten Druckkräften in der z- 
Richtung (Bild 1) für die Komponenten des Vorspannungstensors die Beziehungen gelten: 


Tadel Dee Del. =0 anti n(20). 


Die in den Gl. (19 enthaltenen Aussagen für den gewählten Ansatz ergeben sich durch 
Einsetzen der entsprechenden Ausdrücke aus dem Gleichungssystem (17). Hieraus folgen 
unter Umgehung der Zwischenrechnung die beiden nachstehenden Bedingungsgleichungen: 

— Einy+l(aSiny—yYCojy)e=0 | 


ME A(OT 
g—-1)&ofy+f{a—1+g1—20)]Coy+(g— 1)ySinye=0| Si 


Es handelt sich dabei um zwei Gleichungen für die Konstante c, die nur dann eine nicht- 
triviale Lösung besitzen, wenn die Determinante verschwindet. Hieraus folgt: 

Ra 17, FI SyPRaE (12 y=0% .n...2...%..010% 3 (22). 

Die Auflösung nach g liefert: 

2y — Sin(2y) 
) } 3% 
9=;5 een; (23). 
2y + (a — 1) ©in(2y) 

Wird noch beachtet, daß der Quotient 5 G = 
der Vordehnung des Werkstoffes entspricht und in- 
folgedessen außerordentlich klein gegenüber 1 ist, so 
folgt, daß y sehr klein sein wird. Es ist daher nach 
diesem Parameter zu entwickeln und man erhält bei 
Berücksichtigung des ersten, von Null verschiedenen 
Gliedes die Übereinstimmung mit der Eulerschen 
Knickspannung der Platte: 

272 E 
I; T=ET & 


1 
12(1 ) 12 


Damit ist die Brauchbarkeit des Verfahrens 
für den gedrückten prismatischen Stab nachgewie- 
sen. Es sei nochmals darauf hingewiesen, daß sämt- 
liche Differentiälgleichungen aus der gewöhnlichen 
Elastizitätstheorie entnommen werden können, bzw. 
nach bekannten Methoden, insbesondere mittels des 
Dreifunktionensatzes, integriert werden können. Vor 
allem tritt der mit dem Vorspannungstensor in Zu- 
sammenhang stehende Parameter erst bei For- 
mulierung der Randbedingungen auf. 


il, 


Bild 1 


4. Anwendung auf die Theorie der Plattenknickung 
Die sinngemäße Anwendung auf die Plattentheorie führt dazu, daß die Differential- 
gleichungen der gewöhnlichen Plattentheorie auch für die Plattenknickung Gültigkeit behalten. 
Allerdings muß das an der Plattenoberfläche angreifende Kräftesystem hierbei vollständig 
(d.h. sowohl die Normal- wie auch die Tangentialkomponenten der Oberflächenkräfte) berück- 
sichtigt werden. In tensorieller Form ergibt sich das zugehörige Gleichungssystem aus der 
allgemeinen Schalentheorie des Verfassers [7], wenn der Krümmungstensor b*? der Schalen- 

mittelfläche vernachlässigt wird. Dann gilt die Gleichgewichtsbedingung: 


NW. tmk HP =0.:. 2... (25). 
Der Tensor der Biegemomente steht hierbei in folgendem Zusammenhang mit den 
Ableitungen der Verschiebung senkrecht zur Platte: 


Gh? axß ; 
mer - HE wer— _ Wr) ee ae > (26), 
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Nach Einsetzen dieses Ausdruckes in Gl. (25) ergibt sich folgende Differentialgleichung 
für die Plattendurchbiegung: 
m Gh° 
To 
Für die Komponenten des Spannungstensors an der Plattenoberfläche gilt unter Voraus- 


setzung konstanter Plattendicke und bei Vernachlässigung von Gliedern höherer Ordnung [7] 
Gl. (85)!) 


wis + m, +9 == (re re (27). 


Ba a ne ne 


Im Falle der Knickung der Platte unter einem Vorspannungstensor mit zweidimensiona- 
len Komponenten, d.h. mit 


1320, TEE (ou: 320 ee 
lauten die Randbedingungen für die Plattenoberfläche: 
De 0, 133 728 (cn Be re 
Unter Beachtung der Gl. (11) ergibt sich hieraus: 
18 — TeW, =0 für y-+t re) 
Bei Beachtung von Gl. (28) folgt: 
m? — h TAB Wie). all. Tears Re 1223 


Nach Einsetzen in Gl. (27) ergibt sich, wenn gleichzeitig beachtet wird, daß die Platten- 
oberfläche unbelastet sein soll (K® —p —=0): 


1 
al) 


AAW 7 


Ta, „EE0 Tea IN DEE ST 


Hierbei ist zur Abkürzung 


(a ee 


gesetzt. Gl. (33) steht in Übereinstimmung mit den bisherigen Formulierungen des Problems 
der Plattenknickung, insbesondere den Untersuchungen von K.Marguerre mit Hilfe der 
energetischen Methode [5] und den zahlreichen Arbeiten von Fr. A. Willers, wobei die 
Plattenknickung als Variationsproblem aufgefaßt ist und auf diese Weise die Anwendung auf 
kompliziertere Fälle der technischen Praxis ermöglicht wurde [6]. 

Die dargelegten Gedankengänge führen auch zur Formulierung einer vereinfachten 


Stabilitätstheorie der elastischen Schalen, worüber in einer späteren Arbeit berichtet werden 
soll. 
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D. Elastizität und Plastizität 


Über eine Verbesserung des plastischen Druckversuchs durch 
aufgezwungene Gleitreibuns!) 


Von E. Mönch in München 


Anläßlich der Untersuchung der Materialeigenschaften von Zelluloid für den 
plastischen spannungsoptischen Versuch hat der Vortragende ein neuartiges Verfahren des 
Druckversuchs erprobt. Beim Druckversuch besteht die bekannte Schwierigkeit, daß an den 
gedrückten Endflächen des Probekörpers die Reibung dessen seitliche Ausdehnung behindert. 
Die Folge davon ist, daß der beabsichtigte gleichförmige einachsige Spannungs- und Verfor- 
mungszustand von den Druckflächen her gestört wird. Die freie Oberfläche baucht sich nach 
außen aus. Zu den angreifenden äußeren Spannungen in axialer Richtung kommen die 
seitlichen Komponenten der Reibung hinzu, so daß die resultierenden äußeren Spannungen 
nicht, wie es sein sollte, parallel, sondern schräg nach innen gerichtet sind (Bild la). Nament- 


a b 


Bild 1. Lastangriff beim Druckversuch (schematisch). a) ohne, b) mit Druckflächenverschiebung 


lich bei plastischen Versuchen, wo die Probenhöhe wegen der Ausweichgefahr klein gehalten 
werden muß, ist es, auch bei Schmierung der Druckflächen, unmöglich auch nur annähernd 
Gleichförmigkeit zu erzielen, wenn nicht besondere Maßnahmen ergriffen werden. Ein Ver- 
fahren des Druckversuchs, das die besagten Schwierigkeiten weitgehend behebt, ist das von 
Siebelund Pomp 1927 angegebene sog. Kegelstauchverfahren, das sich jedoch für den 
vorliegenden Zweck aus verschiedenen Gründen nicht empfahl. Das neue Verfahren des Vor- 
tragenden besteht darin, daß man die Druckplatten zwangläufig mit der Stauchung auf dem 
Probekörper gleiten läßt, auf Ober- und Unterseite in entgegengesetzter Richtung (Bild 1b). 
Dann haben die Reibungskräfte auf der ganzen Druckfläche gleiche Richtung, im Gegensatz 
zu Fall @«. Die resultierenden äußeren Spannungen sind daher, wenngleich ebenfalls um den 
Reibungswinkel o gegen die Druckflächennormale geneigt, sämtlich jedoch parallel gerichtet. 
Gibt man nun schließlich den seitlichen Begrenzungsflächen des Probekörpers ebenfalls die 
Neigung 0, so wird er in allen Teilen auf reinen einachsigen Druck in der durch o gegebenen 
Richtung beansprucht. Senkrecht zu ihr wirken keinerlei Spannungen. Die Querdehnung 
kann also nicht behindert werden, wenn es gelingt wirklich gleichmäßige Reibung längs der 
Druckflächen zu erzielen. Es können daher dann auch ziemlich breite Probekörper verwendet 
werden. 

In der zur Durchführung des Verfahrens gebauten Versuchsanordnung blieb die obere 
Druckplatte in Ruhe, während die untere während des Stauchvorgangs zugleich seitlich bewegt 
wurde. Eine zwangläufig mit der Bewegung der unteren Druckplatte arbeitende Führung sorgte 
dafür, daß der Probekörper sich mit der halben Geschwindigkeit bewegte. Diese Führung ist 
notwendig. Ohne sie würde der Probekörper auf einer der Druckflächen haften’ bleiben. Die 
Druckflächen wurden mit Paraffin geschmiert. 

Das während des Stauchversuchs beobachtete spannungsoptische Isochromatenbild, 
das zur Messung des optischen Effektes diente, gibt gleichzeitig Aufschluß über den Grad der 
erzielten Gleichmäßigkeit des Spannungszustandes. Bild 2 zeigt die Isochromaten, die bei der 
Stauchung von Zellulojdproben, mit und ohne Druckflächenverschiebung, unter sonst gleichen 
Versuchsbedingungen beobachtet wurden. Die bessere Gleichmäßigkeit beim Versuch mit 
Druckflächenverschiebung ist offensichtlich. Hier weicht der optische Effekt im wesentlichen 


1) Eine ausführliche Veröffentlichung wird in der Zeitschrift für angewandte Physik erscheinen, 
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nur um eine Isochromatenordnung von der höchsten Ordnung 27% ab, während beim gewöhn- 
lichen Stauchversuch (a) der Unterschied 7 Ordnungen beträgt. 


Bild 2. Isochromatenbilder von Zelluloidproben bei einer Stauchung von 30 % 
a) ohne, b) mit Druckflächenverschiebung 


Die mit dem neuen Verfahren gemessene Spannungs-Dehnungskurve und ihr Vergleich 
mit derjenigen für Zug ergab das für die plastische Spannungsoptik wichtige Ergebnis, daß bei 
Zelluloid diese Kurven den entsprechenden für Stahl weitgehend ähnlich sind. 


Bei Versuchen mit Aluminium fielen die Spannungs-Dehnungskurven für Druck und Zug 
im vergleichbaren Bereich fast völlig zusammen, in Übereinstimmung mit den Messungen, die 
Siebelund Pomp mit dem Kegelstauchverfahren durchführten 


Der senkrecht zu seiner Ebene belastete, elastisch gebettete 
Kreisringträger 


Von Karl Federhofer in Graz 


Die Formänderung eines senkrecht zu seiner Ebene belasteten Kreisringträgers ist gekenn- 
zeichnet durch die Verschiebung v(p) des Querschnittschwerpunktes in Richtung der Nor- 
malen zur Ringebene und durch die Verdrehung ß (g) des Ringquerschnittes, wobei der Winkel 
p die jeweilige” Lage des Ringquerschnittes festlegt, Der Träger sei durch eine beliebig mit 
veränderliche Belastung 9 (p) je Längeneinheit des Ringes und durch Drehmomente © (9) 
belastet und sei in einem Medium gebettet, das jeder Verschiebung v einen elastischen Wider- 
stand c,v(Winklersche Annahme) und jeder Verdrehung ß ein Moment c, ß entgegensetzt 
(elastische Drehbettung). 


Die Differentialgleichungen für v und ß gewinnt man nach der Energiemethode, indem 
das Gesamtpotential y zu einem Extremum gemacht wird. Dieses setzt sich zusammen aus der 
Formänderungsarbeit der Biegung und Torsion (einschließlich der Wölbkrafttorsion), aus dem 
Arbeitsbeitrage der elastischen Widerstände und aus der Arbeit der gegebenen Belastung q 
und ©. Solche als Fundamente dienende Kreisringträger haben rechteckige Querschnitte, 
bei denen der Einfluß der Querschnittsverwölbung, wie genauere Rechnungen t) zeigten, sehr 
klein ausfällt; beim Rechteckquerschnitt alt auch die Berücksichtigung der Wölbmomente 
des Querschnittes. 


Die Forderung öy —=0 führt dann zu den simultanen Gleichungen 


j Re: 3 
year — ypyD Lay (1+») BD = BT, | 


EJ,' 
| B 9a: De 1 AREAL) 
Army mHmEd —- (1+)B=— —— | 
EJ, 
$ ZEN v GJn 
worin mit « als Halbmesser der Ringachse: y= — , w= BJ, 2 (Verhältnis der EN Hin zur 
4 
Biegungssteifigkeit), c,, €, die dimensionslosen Betlungsziffern © Guss ET, m 57. bedeuten. 


Durch Beseitigung von ß folgt als Grundgleichung für die Berechnung des elastisch 
gebetteten und drehelastisch gelagerten Kreisträgers bei Belastung durch g(p) und @ (p): 


ra Ay HALS) yE-Al+S)y=TQ) :..2.2 


')K.Federhofer, Sitz.-Ber. Akad. Wiss. Wien, Abt. IIa, Bd. 157 (1949), S. 307. 
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hierin bedeutet /'(p) die Belastungsfunktion 
Ma qga(ll+c O'’a(1-+A) 1 


IS Se m — 

(9) BJ, GT, EI. und 4 uf 
Mit der Lösung dieser Differentialgleichung ist nach den Gln. (1) auch ß bestimmt und zwar 
durch 


Be Eu 5 ange 9a(1+9) 
(1 in 4) (1 t %)P Y Ar (2 Fi 7) Y t 1% EJ, EJ, 
Der praktisch wichtige Sonderfall des nur durch einander gleiche und voneinander gleichweit 
entfernte lotrechte Lasten belasteten Kreisträgers, für den g=0, 0 =(, somit T= 0, ist 
mit der Annahme c, —=0 ohne Herleitung der Grundgleichung (2) bereits kürzlich von E. 
Volterra?)ausführlich unter Beigabe numerischer Ergebnisse behandelt worden. 
Bezüglich der Einzelheiten der vorliegenden Arbeit und einer im Hinblick auf die 
Unsicherheit der Angaben für die Bettungsziffern ausreichenden Näherungslösung der Gln. (1) 
sei auf meinen Beitrag zum v. MisesAnniversary?) Volume verwiesen. 


Bruchflächentheorie für das ebene anisotrope Plastizitätsproblem ') 
Von H. Neuber in Dresden 


1. Übersicht 

In einer früheren Arbeit ?) zeigte der Verfasser, daß sich die allgemeine Lösung des 
ebenen Plastizitätsproblems bei beliebigem anisotropen Fließgesetz auf graphischem Wege mit 
Hilfe eines von dem Spannungskomponenten gebildeten Koordinatensystems finden läßt. Die 
ebene Plastizitätsbedingung er- 
scheint hierbei als Fläche, deren j Grdy-P 
45°-Böschungslinien mit den 
Charakteristiken zusammen- 
fallen. Die fürdie Bruchrichtung 
im isotropen Medium von 
Mohr aufgestellten Überle- 
gungen, welche auf die Mohr- 
sche Grenzkurve führen, rei- 
chen im anisotropen Medium 
nicht mehr aus, da ein zusätz- N_ 


Bild 1. Beispiel einer anisotropen Plastizitätsfläche in axonometrischer Darstellung 


licher Richtungseffekt berücksichtigt werden muß. Zur Ermittlung der Bruchrichtung im 
anisotropen Medium wird nunmehr eine Erweiterung der Mohrschen Gedankengänge 


2)E. Volterra, Journ. Appl. Mech. 19 (1952), S. 1—4. 

3) Academic Preß, Inc. 

1) Der Vortrag war ursprünglich für den internationalen Kongreß für Mechanik in Istanbul angemel- 
det, an dem der Verfasser nicht teilnehmen konnte. 

2)H. Neuber, Z. angew. Math. Mech. 28 (1948), S. 253—257. 
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vorgenommen, wobei an die Stelle der Mohrschen Grenzkurve eine Grenzfläche tritt. 


Im Endergebnis wird der Satz bewiesen, daß die Charakteristiken mit den Bruchflächen 
zusammenfallen. 


2. Die geometrischen Zusammenhänge als Ergebnis der Integrationstheorie 
der Gleichgewiehtsbedingungen 

Zur Einführung seien die Ergebnisse der Integrationstheorie des Verfassers?)im Hin- 
blick auf ihre geometrischen Folgerungen kurz dargestellt. Es erwies sich als zweckmäßig, 
ein Koordinatensystem zu verwenden, das von der Summe der Normalspannungen o,-+-o, 
der Differenz der Normalspannungen o,— co, und der doppelten Schubspannung 27,, ge- 
bildet wird. In diesem Koordinatensystem erscheint die Plastizitätsbedingung als Fläche 
(‚„‚Plastizitätsfläche‘‘); ein Beispiel zeigt Bild 1, die zugehörige Grundrißprojektion Bild 2. 
Für den Sonderfall der Isotropie ergibt sich eine Drehfläche, während bei Anisotropie die 
Drehsymmetrie verloren geht und nur eine Symmetrieebene erhalten bleibt (als Folge der 
Benutzung des doppelten Winkels der Hauptspannungsrichtung). 


5 4 0x-04-Q9 005(2 9) 
ee 5 
a 
EEE = —r 


\ N | Er BZ 


27-g sin(2%) 


N 


N 
Bild 2. Grundrißprojektion der anisotropen Plastizitätsfläche 
Sind o, und o, die Hauptspannungen und 
ter und er 


ihre Summe bzw. Differenz, so folgt aus der erwähnten Arbeit des Verfassers für die Spannungen 
längs der Charakteristiken «, v mit den Bezeichnungen von Bild 3 und 4: 


00, 00, Be a) 00, 00, E 5 
- = = 0 = — = == ® 
ou eu ou o W ev y 2) 
oder allgemein 
(do,)=i(do,) —.(drsy)e= 0,22 2 re IRRE!) 


bzw. 
Id(o, + 0,)]? = [d(o, — 0,)]? #A (dr...) See Ei 
Beim Übergung auf die Hauptspannungen ist | 
a t+Yy=P, %—y=gcs(2p), 2my=gsinl2p) . -.....0). 
Die Koordinaten 7, g, 29 bilden Zylinderkoordinaten. 
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Aus Gl. (4) ergibt sich 


ka ya A a tn er (6). 
Demnach sind die Charakteristiken die 45°-Böschungslinien der Plastizitätsfläche. 
Für den Winkel 2 «ihrer auf die Grundrißebene projizierten Tangente mit der (o, — 0,)- 


Achse folgt 
do, Tyy 
Bean Ydd, 


bzw. 


vgl. hierzu Bild 3. Ist ferner 25 der Winkel der 
auf die Grundrißebene projizierten Tangente mit 
der Meridianebene, welche ihrerseits den Winkel 29 
mit der (o, — 0,)-Achse bildet, so gilt 


ee 


U 
Wekonst) 


Bild 3. Bezeichnungen beim Übergang auf beliebige Koordinaten, sowie Bild 4. Ermittlung des Tangentenwinkels 
auf die Hauptlinien 


Diese Winkel finden sich sämtlich im Körper selbst wieder, und zwar ist (entspr. Bild 3) a der 
Winkel zwischen der Charakteristik und der &-Achse, ß der Winkel zwischen der Charakteristik 
und der Hauptspannungsrichtung. 

Beim Übergang von der z, y-Ebene des Körpers auf die Grundrißebene der Plastizitäts- 
fläche gilt also die „‚Abbildungsvorschrift‘‘, daß die Winkel 9, a, ß verdoppelt werden; ferner 
ist beim Ablesen der Winkel «x und feine Vertauschung der beiden Charakteristiken (u —= konst. 
mit v = konst.) zu beachten, d.h. der Tangentenwinkel an eine Linie % = konst. in der z, Y- 
Ebene ist gleich dem halben Tangentenwinkel an die Grundrißprojektion der 45°-Böschungs- 
linie v — konst. in der (o, — 0,), (2 7,„)-Ebene. 

® 
3. Die Mohrsche Theorie für den isotropen Stoff 

Die Vorstellung der Bruchlinien, die ja auch durch das Experiment bestätigt ist’), geht 
auf grundlegende Untersuchungen von Mohr) zurück. Hiernach gibt es für eine Bean- 
spruchung, die an der Plastizitätsgrenze liegt, zwei ganz bestimmte Schnittflächen in welchen 
die Schnittspannungen bei geringfügiger Überschreitung der Belastung die Zerstörung des 
Stoffzusammenhanges herbeiführen. Da man früher meist die Schubspannung als Haupt- 
ursache des Bruches ansah (Schubspannungshypothese), sprach Mohr von ‚Gleitung‘ in 
diesen Schnittflächen, nannte jene selbst ‚,Gleitflächen‘“, obwohl diese Bezeichnungsart eigent- 
lich mit einer Beschränkung des ebenen Plastizitätsproblems auf den Gültigkeitsbereich der 
Schubspannungshypothese verbunden ist. In Wirklichkeit sindin jenen Schnittflächen außer 
Gleitungen (d.h. Tangentialverschiebungen) auch Normalverschiebungen der Schnittufer 
möglich, so daß es sich um ‚„‚Bruchflächen‘“ handelt. 

Wird zunächst eine beliebige Schnittfläche betrachtet, welche mit der Hauptspannungs- 
richtung den Winkel ß bilden möge und an welcher die Normalspannung o und die Schub- 
spannung 7 angreifen (Bild 5), so folgt aus der sinngemäßen Anwendung der Gleichungen (5), 


3) A. Nadai, Der bildsame Zustand der Werkstoffe. Berlin 1927, 
4) O. Mohr, Z. VDI 44 (1900), S. 1528—1530. 
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indem an die Stelle der Schnittebene % = konst. die hier betrachtete Schnittebene tritt, d.h. 
indem — ß statt , o statt o, und 7 statt T,, gesetzt wird: 
20=p—geos(2Pß), 27=-.gsin (2 BP) Fi Arena 

Wird ferner zunächst nur der isotrope Stoff betrachtet, so entsprechen jedem 
Wertesystem p, q, welches der isotropen Plastizitätsbedingung genügt (bzw. auf der Plasti- 
zitätsfläche des betreffenden isotropen Stoffes gegeben ist), zugehörige Werte von (20) und 
(2 r), welche in einem rechtwinkligen, von (20) und (2 r) gebildeten Koordinatensystem auf 
einem Kreise liegen, dem sogenannten Mohrschen Spannungskreis des 
isotropen Körpers (Bild 6), dessen Mittelpunkt sich auf der (2 0)-Achse im Abstande 9 be- 
findet und dessen Radius gleich q ist. 

Alle Spannungskreise, welche im Einklang mit der isotropen Plastizitätsbedingung 
konstruiert werden können, haben dann nach Mohr die Eigenschaft, daß sie eine Grenz- 
kurve (die Mohrsche Grenzkurvedessisotroap-plastischen Körpens; 
Bild 6) berühren, was aus folgender Betrachtung hervorgeht: 

Die der Grenzkurve entsprechenden Werte der Schnittspannungen sind jene Grenzwerte, 
welche gerade an der Grenze des plastischen Verhaltens des Stoffes liegen ; sie seien mit o, und 
7, bezeichnet. Die zugehörige Schnittfläche ist die ‚„‚„Bruchfläche‘“, der zugehörige Winkel 
gegenüber der Hauptspannungsrichtung 1 sei mit ß, bezeichnet (Bild 5). Da o, und 7, den 


Gleitebene 


{ 


(Hauptspannungsrichlung) 


Bild 5. Bezeichnungen zur Bruchflächentheorie Bild6. Die Mohrsche Grenzkurve des isotropen Körpers 


Grenzwerten entsprechen, so darf für keine andere Schnittfläche, d.h. an keiner anderen 
Stelle des jeweiligen Spannungskreises die Grenzkurve erreicht oder überschritten werden 
(der untere Zweig der Grenzkurve verläuft spiegelbildlich zur (2 0)-Achse). An der Stelle 
o=0,Tr=r,trittdahereine Berührung mitder Grenzkurveein, d.h. die Grenz- 
<kurveist Hüllkurve aller Spannungskwesse „dıezin rede manee 
trachteten isotrop-plastischen Körper möglichsind (Bild 6). 


4. Allgemeine Bruchflächentheorie für den anisotrop-plastischen Stoff 


Für den anisotropen Körper ist diese Darstellung nicht mehr brauch- 
bar, da bei festem p der Radius gesmitr der FTaupispan nee 
richtung veränderlich ist. Deshalb muß eine sinngemäße Erweiterung und Ver- 
allgemeinerung der Mohrschen Vorstellung vorgenommen werden. Aus dem Einfluß des 
Schnittwinkels auf die Länge des Radius geht hervor, daß die zu einem festen p-Werte ge- 
hörenden Spannungswerte bei fester Bruchfläche nicht auf einem Kreise, sondern auf einer 
ellipsenähnlichen geschlossenen Kurve liegen, welche punktweise aus der betreffenden Kurve 
p — konst. der Plastizitätsfläche oder ihrer Grundrißfigur (Bild 2) konstruiert werden muß. 
Für eine andere Bruchrichtung erscheinen die ellipsenähnlichen Spannungskurven gedreht, so 
daß sich eine andere Grenzkurve ergibt. Hieraus folgt, daß man nicht mehr mit einer 
Grenzkurve auskommt, sondern unendlich viele Grenzkurven benötigt, d.h. an die 
Stelle Grenzkurve tritt eine Grenzfläche. 

Zu den beiden Koordinaten o und 7 des Grenzspannungs- oder Bruchrich- 
tungsdiagrammes von Mohr muß mithin folgerichtig als dritte Koordinate der 
Winkel & hinzugenommen werden, welcher de Lage der Schnittebene charakteri- 
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siert, an welcher die Spannungskomponenten o und r angreifen. Werden Zylinderkoordinaten 
gewählt, mit (2 0) in Achsrichtung (27) in Radialrichtung und (2) als Breitenwinkel, so er- 
gibt sich ein anschauliches Bild der anisotropen Grenzfläche; eine solche ist in Bild 7 veran- 
schaulicht. 

Mathematisch besteht nunmehr folgender Sachverhalt: Die drei für die Spannungen in 
einer Schnittfläche maßgeblichen Größen o, rt, x sind Funktionen der Größen 9, 9,9 und ß, 
welche die Rolle von Parametern spielen, wobei p, qg und p noch durch die Plastizitätsbedin- 
gung, bzw. durch die Plastizitätsfläche funktional verknüpft sind. In der Tat bestehen außer 
der Plastizitätsbedingung noch 3 Beziehungen zwischen diesen Größen, und zwar mit 
Bezug auf Gl. (10) und Bild 5: 


20=p—gcos2ß, 27r=—qsin2ß, Dean ee. (11). 


Bei gegebenen Werten 9, 9, ß 
und Beachtung der durch die 
Plastizitätsfläche festgelegten 
funktionalen Verknüpfung von 
?, qund @ können o, r und a 
ermittelt werden. Dieser Sach- 
verhalt besteht für jede belie- . 
bige Schnittfläche. Die entspre- 
chenden Spannungskurven bil- 
den jeweils eine Fläche, die als 
„Spannungsfläche‘“ bezeichnet 
werden möge. In Bild 7sind die 
Spuren von zwei solchen Span- 
nungsflächen (p—=4 und p=0) 
in den Ebenen 2x —=0 und 

7 i 
2 eg; axonometrisch veran- 
schaulicht. Indem nunmehr ver- 
langt wird, daß alle möglichen 
Wertetripel o, 7, innerhalb der 
Grenzfläche liegen sollen, so 
folgt, daß alle Spannungsflächen 
diese Grenzfläche gerade be- 
rühren; d.h. die Grenz- 
fläche ist Hüllfläche aller 
durch mögliche Wertetripel o, 
t, & gegebenen Flächen. 

Das Wertetripel 0 — 0: 
T= Tr, = &g ZU welchen die Bild 7. Beispiel einer anisotropen Grenzfläche 
Parameterwerte 9, q, 9, ß ge- 
hören, stellt daher den gemeinsamen Berührpunkt zweier unendlich benachbarter Spannungs- 
flächen dar und kann daher zugleich auch durch die Parameterwerte p + dp, q + dq, 9 + dp, 
ß + dß dargestellt weren; d.h. die Ableitungen der Gl. (11) nach diesen Parametern ver- 


schwinden: 
d(2 0) = dp — dq cos (2B) + 2qdßsin 2P)=0, 
d(2 7) = — dqsin 2 BP) — 2 gdß cos (2P)—0 Be (12) 
d(20)=d(2p)+d4@P)=0 | 


Aus der dritten Gleichung folgt zunächst 


Ra ea a (©) 
mithin ergibt sich aus der zweiten Gl. (12) 
2 gdp 
tg (2 ß) = —— (14). 
AI er 
Durch Einsetzen in die erste der Gleichungen (12) ergibt sich nach Umformung 
dq — dp cos Ne sungletanh, wen (15). 


Andererseits folgt hieraus mit Bezug auf Gl. (14) 
2 q dp = dp sin (2) Da (16). 
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Durch Quadieren und Addieren ergibt sich 

(dp)? = l@0)?. (2 2,doPr ve EEE 
Der Vergleich mit den für die Charakteristiken gewonnenen Ergebnissen führt demnach zu der 
Festellung: 

Die, «Bruch hin jen. des Haklg ei meinen ne be nenn an TER Pe 
plastischen Zustandes fallen mitrden CGharakterisitikenzüer 
Differentialgleichungen des Gleichgewichtes’zusammen’ und 
sind zugleich die 45°-Böschungslinien der Plastizitätsfläche., 

Es ist bemerkenswert, daß die Bruchflächentheorie als solche mit dem beschrittenen Wege 
der Parametervariation aufgestellt werden konnte, ohne auf die Gleichgewichtsdifferential- 
gleichungen Bezug zu nehmen. Allerdings führen erst jene zu der ergänzenden Feststellung, 
daß die Bruchlinien mit den Charakteristiken zusammenfallen, welche für die Integration von 
ausschlaggebender Bedeutung sind, sowie zu dem wichtigen Gesetz der Vertauschung und Ver- 
dopplung der Tangentenwinkel bei der Abbildung der x, y-Ebene, auf die’(o, — 0,), (2 T,,)- 
Ebene. Es sei noch bemerkt, daß man zwar zu den Gl. (15—18) formal auch ohne Berück- 
sichtigung der Anisotropie gelangen kann ?). Dies darf aber nicht zu dem Trugschluß führen, 
daß die Anisotropie bei Bestimmung der Bruchrichtung außer Acht gelassen werden könnte. 
Eine solche Schlußweise ist grundsätzlich falsch, weil hierbei der durch den funktionalen Zu- 
sammenhang der Hauptspannungen mit der Hauptspannungsrichtung hinzukommende 
Richtungseffekt völlig vernachlässigt werden würde. Auf die Abweichung der durch die 
quasi-isotrope Auffassung ermittelten Bruchrichtung von ihrer wirklichen lL.age wird im nächsten 
Abschnitt an Hand eines Beispieles noch näher eingegangen. Andererseits sei an dieser Stelle 
noch ausdrücklich darauf hingewiesen, daß die schon lange bekannten Beziehungen der 
Spannungsdifferentiale beider Mohrschen Grenzkurve sich infolge ihres speziellen Charak- 
ters stets auf den isotropen Stoff bezogen haben und auch niemals die Ergebnisse der Inte- 
grationstheorie für beliebiges isotropes oder anisotropes Fließgesetz vorwegnehmen Konnten. 
Die Aufstellung einer allgemeinen Bruchflächentheoriesundder 
allgemeine Beweis für die formale Übereinstimmungihrer Er- 
gebnisse mit der Integrationstheorie sind erst: mitider vor- 
liegenden Acpertdargelegt. 


5. Abweichung bei Ermittlung der Bruchrichtung nach der quasiisotropen 
Vorstellung gegenüber der Wirklichkeit 
Die erwähnte formale Übereinstimmung einiger Differentialbeziehungen für den iso- 
tropen Stoff mit den allgemeineren Gleichungen des anisotrop-plastischen Mediums darf den 
Praktiker nicht dazu verleiten, bei Bestimmung der Bruchrichtung von der Vorstellung des 


p- + 


Exakte Gleitrichtungen 


p-1 

p*2 . 7 
23 
p=4 
PR | 
N 
= H, P kKonst, exakt / 
0 P = Konst., Mäherung (quası -1solrop) 


Bild 8. Ermittlung der Bruchrichtung; Darstellung der Abweichungen der quasi-isotropen Näherungsmethode von der exakten Methode 


isotrop-plastischen Mediums auszugehen und in bekannter Weise die Mohrsche Grenzkurve 
anzuwenden. Es erscheint insbesondere deshalb wichtig, auf die so entstehenden Fehler hin- 
zuweisen, weil sich in der Bergbaupraxis dieses Verfahren seit längerer Zeit eingebürgert hat. 
Entnimmt man aus einem bestimmten Stoff (z.B. einer Gesteinsart mit ausgesprochen 
anisotropen Festigkeitseigenschaften) prismatische oder zylindrische Proben in einer be- 
stimmten Richtung, so erhält man aus den zugehörigen Festigkeitsversuchen Aussagen 


5) H. Geiringer, Z. angew. Math. Mech. 32 (1952), S. 379; C. Torre, Österr. Ing. Arch. I 
(1946), S. 36 u. 316; IV (1950), S. 174; VI (1953), S. 32. 
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für die Bestimmung einer funktionalen Beziehung zwischen den Hauptspannungen für die 
betreffende Hauptspannungsrichtung. Man erhält auf diese Weise nur die zu dem be- 
treffenden Hauptspannungswinkel gehörende Meridiankurve der Plastizitätsfläche. Be- 
schränkt man sich auf solche Materialproben und entnimmt keine weiteren Proben in anderen 
Richtungen, sondern betrachtet das Medium als quasi-isotrop, so würde die zugehörige Plasti- 
zitätsfläche eine Drehfläche sein. Die Höhenlinien p — konst. wären sodann Kreise, welche 
lediglich auf der speziellen Meridiankurve mit den wirklichen p-Werten zur Deckung kommen. 
Die Konstruktion der Böschungslinien an dieser Ersatzdrehfläche führt offensichtlich zu 
wesentlich anderen Bruchrichtungen, als den Böschungslinien der Plastizitätsfläche entsprechen. 
Die Zusammenhänge sind im Bild 8 veranschaulicht. Man erkennt deutlich den großen Fehler 
bei der Bruchrichtungsermittlung nach der Näherungsmethode gegenüber der exakten Be- 
stimmung. 


Eingegangen am 9. Juni 1953. 


Neue Methoden zur Bestimmung von Fließkurven und Relaxationszeiten 
an Rohr- und Rotationsviskosimetern und Diskussion des Prandtl- 
schen Gesetzes an Hand experimenteller Unterlagen 


Von F. Schultz-Grunow und H. Weymann in Aachen 


Zur Messung von exakten Fließkurven bei nichtnewtonschen Flüssigkeiten sind von den 
üblichen Viskosimetern nur das Rohr- und das Rotationsviskosimeter brauchbar. 

Beim Rohrviskosimeter war bisher wegen der Einlaufeffekte der Druckgradient nicht 
exakt bestimmbar. Es wird über ein neuartiges Rohrviskosimeter berichtet, in dem der Druck 
über die Deformation, die das Rohr erleidet, an beliebig vielen Stellen gemessen werden kann. 
Zu dieser Messung werden elektrische Dehnungsmeßstreifen benutzt, die bis zu Dehnungen 
von 10 eine Genauigkeit von 3% haben. Je nach Wahl der Meßrichtung kann der statische 
Druck oder die Schubspannung bestimmt werden; eine wichtige Tatsache, da in allgemein 
viskosen Flüssigkeiten der achsiale und der radiale Druck nicht notwendig gleich sind. Mit 
diesem Rohrviskosimeter lassen sich somit Druckgradient, Einlaufeffekte, Schub- und Normal- 
spannungen messen. Weiter besteht die Möglichkeit, die Abhängigkeit der Viskosität vom 
Druck zu untersuchen. 

Es wird dann daraufhingewiesen, daß bei der bisherigen Bestimmung von Fließkurven 
unzulässige Näherungen und Mittelbildungen verwendet worden sind. Die exakte Auswertung 
von Rohrversuchen kann mit Hilfe der Formel von Rabinowitsch erfolgen, für das 
Rotationsviskosimeter hat der Vortragende eine ähnliche Formel abgeleitet, 


Bu le.,da 
R r) =2M +) 


die als Rekursionsformel für das Fließgesetz 


=D 
zu verstehen ist. 

Es wird anschließend auf die Prandtlsche Formel D = ( ©in r/A für normalzähe 
Stoffe eingegangen. Die Berechnung der Strömung einer normalzähen Flüssigkeit im Rohr 
war bereits von Prandtlund Vandrey durchgeführt worden. Die entsprechende vom 
Vortragenden durchgeführte exakte Rechnung für das Rotationsviskosimeter für beliebige 
Spaltbreiten wurde sodann kurz skizziert und anschließend Meßergebnisse mitgeteilt, die mit 
Rohr- und Rotationsviskosimetern an Polystyrol-Toluol-Lösungen und einem Acrylsäurebutyl- 
ester gewonnen wurden. Es zeigt sich, daß das Prandtlsche Gesetz für diese Stoffe sehr gut 
bestätigt wird. Anschließend wird die aus der Prandtlschen Theorie folgende Temperatur- 
abhängigkeit der Konstanten A und Ü entwickelt und gezeigt, daß die Beziehung 

20 = Dis ealt 
durch die Versuche ebenfalls sehr gut bestätigt wird. Dann wird der gemessene Verlauf von A 
in Abhängigkeit von der Temperatur bei Polystyrol und dem Acrylsäurebutylester gezeigt, und 
das von der Theorie abweichende Verhalten erklärt. 

Im Anschluß wird über eine Rechnung berichtet, die die Strömung einer Ma xwell- 
schen Flüssigkeit im Rohr einschließlich der Einlaufströmung wiedergibt und auf die Möglich- 
keit hingewiesen, nach dieser Methode die Relaxationszeit Max wellscher Flüssigkeiten zu 
bestimmen. 
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Theorie zähe-fester und firmo-viskoser Stoffe 
Von (C. Torre in Wien 


Zur Untersuchung des Verhaltens fester Stoffe bei langandauernder Beanspruchung wur- 
den diezähe-festen Stoffe durch die Spannungs- Geschwindigkeits- Gleichungen 


(0) zen 
"—9 = EN Met: sw. . 2 3 . . . . 1 
07 242 = +A(V »d), Tıy ae } oe USW (1) 
in Anlehnung an die ER. Ms, der Elastizitätstheorie 
Se. 00 26 136, OD A 
Or 2 6, I en Ex 5) (F Re u) > ya G Er + 5) sTlL ILS-WieZie ie Lie, Soger (2) 


definiert. Der firmo-viskose Stoff wird durch die Summe der Gl. (1) und (2) mathe- 
matisch ausgedrückt, was bei kleinen Verschiebungen (vd = Qu/ot) ergibt 


el). welter. 


Zum Vergleich geben wir noch die Spannungs-Geschwindigkeits-Beziehungen der Hydro- 
dynamikan, die lauten 


m m 0 
Deren > ’ © R vu = Y c h S 2 x x P R 4 h 
Or p > u ax HA(p 2); Tıy U | we =) USW ( ) 
Um Gl. (4) mit Gl. (1) zu vergleichen, bilden wir die mittlere Spannung o,„ nach Gl. (1) 
1 2 
= (taten =[R+z ud ERS REKEN 


Addieren und Subtrahieren wir den Ausdruck Gl. (5) auf der rechten Seite der Gl. (1), dann 
ergeben sich die Normalspannungen in folgender Form 


u, 
= an + 2 ZW), uw). 


Die hydrodynamische Beziehungen Gl. (4) werden mit Gl. (6) erst durch die Annahme von 
Stokes4=— 2 u/3 gleich. Da die hydrodynamischen Berechnungen fast ausnahmslos mit 
dieser Annahme durchgeführt werden, so gelten die Ergebnisse der hydrodynamischen Be- 
rechnungen eher für zähe-feste Stoffe, bei welchen die elastischen Kräfte gegenüber den Rei- 
bungskräften vernachlässigt werden können!). In Gl. (4) kommt noch der Druck (— p) hinzu, 
der den ebenen Verzerrungszustand stets ermöglicht, während bei zähe-festen Stoffen, wie bei 
elastischen, auch ein ebener Spannungszustand möglich ist. Das entspricht dem gegebenen 
Volumen und gegebener Gestalt bei festen Stoffen gegenüber der gegebenen Gestalt und unbe- 
stimmtem Volumen bei tropfbaren Flüssigkeiten. 


DieBewegungsgleichung eines firmo-viskosen Kontinuums lautet 


ou 4 (&- »2) Au -| "+ +)alrw:w, den, 


01? m—2 


Zur Gl. (7) kommen noch die Rand- und Anfangswertbedingungen hinzu, die wir fallweise bei 
gegebenen Beispielen angeben werden. Wir haben mit Gl. (7) als Ausgangsgleichung ver- 
schiedene Fälle behandelt, wie den einachsigen und ebenen Spannungs- und Verzerrungszu- 
stand, die wir an anderer Stelle ausführlich besprechen werden. 


Zwei Probleme der ebenen Elastizitätstheorie 
Von U. Wegner in Heidelberg 
I. Das 1. Problem betrifft die Berechnung der Spannungen in einer Rahmenscheibe 
(Dicke —=1), die durch Zentrifugalkräfte beansprucht ist. In der Praxis interessieren die 
Spannungen an den Stellen 1 und 2 
=), y=H ... kudloksnls au 2 


z-o und g,z-0. 
y=h y=H 


1.12.09, 


1)C. Torre, Österr. Ing.-Archiv 6 (1952), S. 417. 
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Für 1> H — h, also sehr dünne Riegel kann 0, genähert aus der elementaren Balken- 
biegungslehre ermittelt werden. Ist das Verhältnis 1: (H — h) nicht sehr groß, so gibt die 


Balkenbiegungslehre, aus der ae 0,z=0o sich ergibt, falsche Resultate. Gerade bei 
z=R 


” . y=H 
zweifach zusammenhängenden Bereichen und bei der Ermittlung von Randspannungen an 
Stellen, die nicht in der Nähe von Ecken des Bereiches liegen, sind die Variationsmethoden zur 
Berechnung der Spannungen besonders praktisch und genau. Für unser Problem ist 


02, 0? F, 0? F, 
2 i U : — % 2 ni 
0% dy? Sun 0, da2 +U; Toy 0%0y 
mit 
2 
: Dee, 
209 y 
& — Winkelgeschwindigkeit, y =spez. Ge- 
wicht, g = 981 cm/sec?, also 
AAF. =10% 
i a 
% 
TEN 
Ti 
TEEN 
TZSSN 
Bild 1 Bild 2 
ML We. Erden: _rf AN : r 
se ur IT y*ist partikuläres Integral, so daß F,—=F, + F, mit „= 0 wird. Die 


Randbedingungen sind: 


1 Re w?y 
NR y 4g u 


am äußeren Rand 7’, und: 


F,r; — 


ER LLITR B 
54m® a De HAx + B+H 


aminneren Rand J';,, wobei A, Bund CO noch unbestimmte Konstanten sind (bedingt durch den 
zweifach zusammenhängenden Bereich). 
OF, 


on 
folgenden Variationsproblem errechnen, was den Randbedingungen von F,, äquivalent ist: 


an] pdedy+ | (16) £ — g(s) ) ds = Extrem. 
® T4+T. fr 


ü 


—= const. auf T, und 7',. Ist AF,—=o(z, y), so ist Ap = 0 und 9 läßt sich aus dem 


mit: 
fa=Fır+r; > 9(s) = 
Das Randintegral wird für Ax + By + C aber Null, so daß die komplizierte Bestimmung 


der Konstanten A, B C herausfällt. Auf Grund des funktionentheoretischen Verhaltens von 
und der Symmetrie des Problems wird für @ der dreigliedrige Approximationsansatz: 


on Tahlı 


1 2 x y® 
= ot+ca1Ne2? + Ne (4) =o+.(® — y)+% (@ + y) 
gemacht. Dann ist: 


m — 1w?y 


AR=AF + AR= 9 — Ag Y=0o+%,— 2 U 
d2y m+1 
dale: 0, + 0,7 ta, y? Pre 


22 
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IsbzuB. H= 56, re 23h 80, Lesen; solwird 


we y wey 
er 3,2 10 u +0,16 -1102 2,1 
nn 2g a 0) 


Nach der elementaren Balkenbiegungslehre wird: 


2 

wy 

Ou|x—=0 Oulz=0 IT — 0,8 a ln > 
y=h y=H a 


II. Zwei verschieden beheizte gleich große Rechteckscheiben aus dem gleichen Material 
stoßen an einer vollen Seite aneinander. Wie verlaufen die Spannungen an dieser Seite? Wir 
nehmen an, daß das Rechteck ein Halbstreifen der Breite 2 HZ sei. Die Dicke der Scheibe sei h 
26. TW = Temperatur, A = Wärmeleitzahl, o = Wärmeübergangszahl vom Scheiben- 
material zurLuft, t=Zeit. NachMarguerre?) genügt T® bei einer dünnen Scheibe, deren 
Ober- und Unterfläche die Wärme nach dem Newto.nschen Gesetz ableitet, der Gleichung 


r] = 


ATV = 7? (mo — T, + ae T, = Lufttemperatur 


Bei uns sei TO (x, Y; t) Rand : 


TDd(x, HH; )=T, =0; T®(0, y;t)=Acos it A,c=const., u 
T®(o,y;t)=0; und TDd(z, y;0)=4Acos SH: 
Dann ist: 
Er ae r os m? 
Toy = Ans gnann| © da; erg 
Es sei: E 


TV (z, Y; Ü)|Rana; TD(g,+ Ja 0)=0 < TD(—oo;y; DEE; 


TWD(0; y;t)=4Acos U. gi (Randbeheizung) und T®(z, y;0)=0 (z=+D0). 


2H 
Dann ist: 
RAN 
4 ny BIN) 
II a a RA ER EFT 
TOD (x, y;t)= a aH® ee dE. 
0 
Für das Spannungsproblem ist: 
j oa rn e m 2G(m +1) € 
OT aa 2a(e +4) SE mEeT Es 
02 Fi) em) 2G(m +1) 
v) — C m: Te ey ee n 
% 00? a (er u Mm — ) MRNIER: Er 
02. Fw 
WM — — (v).» ——: 1Zi 
9 3003 @G vie & —= Ausdehnungskoeffizient 
v=], 11. gesetzt. 
Dann wird: 
1 
A AF® — er u GxATw; Fo = Fo 4 Fo 
gesetzt, wobei F® ein partikuläres Integral ist; mit 
Ara = — ud )aTn gibt AAFW—0. | 


!) Genauere Approximationen und Diskussion des Spannungsverlaufes s. die demnächst erscheinende Diss. 
von Frl.Fr.Henn, Darmstadt. 


2) K. Marguerre, Ing.-Archiv Bd. 8. 1937, S. 216ff. 
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oF@ OF 
. a 1 . [ 7 . . 
Randbedingungen: F® und FD sowie H und er sind gleich Null auf dem ganzen 


Rand außer dem für &=0. Esist: 
(d n 
F\ ER = Zee ee / y—go.e-odE = fly; ti), 
f) 


oF® ortı ) 


ON zo on 


-| Tayla=0 dE = g(y;t); 
0 


z=0 


da Oz|z—0 und T,,'z„—o unbekannt sind, ist fund g noch unbekannt. Die Variationsprobleme als 
Aquivalent der Randbedingungen für Ff” und F{® Jauten mit gm = AF,:v=LIN) 


+H 
1 od) 
3 | | + army azay — [sin +00 900..J]ay = war. 


®] —n 
und 
1 yz dp(ED 
53 [ | (iD + AFWPdEdy+ [ KR 2) nn —9(y,t) OL) dy = Exstr. 
"Br BT L |z=0 
Ansatz: 
THT TAT 
TRY rY 


+. 
PD gED — 3 B,(t)e ?* cos 3 


ID zoD—= N A,(t er cos 
gm zol ) IH 


He 


Unbekannt bleiben die 2» Fourierkoeffizienten von f(y;t) und g(y;t). Weitere 
Bedingungen 
ud — u für =0 


vd — yD für 20, 
wobei ist: 


u ef 
um — 3G (1 Al AF®(E, y; t) de — Er a | TW&yo)dtw-—yY, 


1 Dr rw i 
al En AFW (a, Lt) dE — dY SER J To n,C;)dE+uy2. 
0 


Es ist 
oFw or oFn or 
E =. und ——— 
0Y |z=0 0Y \z=0 0% 0% 
Somit gibt es An Bedingungen für AF® — g® und AF(D — ID („ — 0 wegen der Symmetrie; 


Ug, %, sind ebenfalls an den obigen Integrationsgrenzen = 0) und damit können die 2» Fou- 
rier koeffizienten eindeutig festgelegt werden. Ansatz wurde für n —= 2 durchgerechnet °). 


z=0 


3) Ausführliche Berechnung im ‚Bauingenieur‘. 


Referate der Vorträge: 
H. Jung, Ein Beitrag zur Plastizitätstheorie, 


H. Wolf, Zur Bestimmung von Abplattungswerten im Erdinnern aus einer 
vorgegebenen Dichteverteilung durch numerische Integration der Clairautschen 
Differentialgleichung 


sind nicht eingegangen. 


& 


22* 
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E. Hydrodynamik 


Lösung eines Randwertproblems der Magneto-Hydrodynamik 
Von R. Nardini in Bologna 
Es wird eine ideale, inkompressible, homogene, leitende Flüssigkeit im Halbraum 2>0 be- 
trachtet, und es werden ihre Bewegungen studiert, die unter der Wirkung eines magnetischen Fel- 
des stattfinden und nur von der Zeit iund der Koordinate zabhängen. Für z2>0 werden der 
Anfangswert der magnetischen Feldstärke 4 als konstant und parallel zur z-Achse und der 
Anfangswert der Geschwindigkeit d als Funktion der Zeit normal zur z-Achse angenommen. 
Auf der Ebene z =0ist H normal zur z-Achse als Funktion der Zeit vorgeschrieben und 
die Komponente von v gleich Null angenommen. 
Bewiesen wird, daß die z-Komponente von H immer mit dem Anfangswert H, überein- 
stimmt und Z,und H,, d.h. x- und y-Komponenten von H, der Gleichung 
0? H,(2, t) de) „02H.(z, 0) uH; 
- rn — + 6b? - ee 
ot? ot 02? 02? uy 0) 
genügen, wo u, y und ö bzw. magnetische Permeabilität, Leitfähigkeit und Dichte sind. Die 
Anfangs- und Randbedingungen liefern 


oH;(z, t) 


(i=2, Y; EI 


lim Z,(2,i)=0 lim —=— =9,() (<z<+ ©) 

0 +0 08 =, 9%), 
lim 7,(2, 1) =@;(t) (07) 

z>+0 


wo ®;(z) und G;(t) bekannt sind, während die Konvergenzbedingung im Unendlichen lautet 


lim H,z,1)=0 (<ı). 
z>+» 


Durch die Anwendung der Laplace-Transformation erhält man die Lösung des obigen 
Problems. Für B;(z) = 0 gilt z. B. 
Er br z 
t s 
H,(z, =) *e a: > FIN Bi ) (—%, Y); 


wo 
B,=y B=y 


t t 
es 2 (t — r)m1 | E ie a Ile A, 
Bam = IC i 2 (m — 1)! Rs Bamıı = U \a“ (m — 1)! Re 


(y und x sind die bekannten Funktionen der Wärmeleitung) sind. 
Die Komponenten von v und der Druck können mittels 7, und H, ausgedrückt werden. 
Man kann endlich hinreichende Bedingungen angeben, so daß die Lösung eindeutig ist. 


Berechnung von Strömungsfunktionen bei Konturen mit Ecken 
Von @untram v. Gorup in Darmstadt 


Die kürzlich in dieser Zeitschrift [1] veröffentlichte Methode zur Berechnung der Strö- 
mungsfunktionen im Außenraum einer zeitlich veränderlichen oder bewegten Kontur wird hier 
für den Fall weitergeführt, daß die Kontur eine Ecke besitzt. 

Dazu muß die Abbildungsfunktion 2 (r) (Abbildung auf den Einheitskreis der 7-Ebene) 
Verzweigungen besitzen, kann also nicht mehr, wie in [1] vorausgesetzt, meromorph sein. 

Hier wird nun gezeigt, wie in dem Sonderfall einer bestimmten Kontur mit rechtem 
Winkel für den Fall der starren Drehung der Kontur die erforderliche additive Zerlegung von 
z z* (z* geht aus z durch Ersatz von r durch 1/7 und i durch —i hervor) in einen im Außen- 
gebiet und einen im Innengebiet des Einheitskreises regulären Anteil gelingt. 

Wir gehen aus von der Abbildungsfunktion: 


(HE 


1i—-— — 
Ta 
Sie liefert eine tropfenförmige Kontur mit rechtem Winkel. 
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Wir bilden 22*: 


3/1 3 3 2 F 
22*=(T+ 1)? | u= ı)' -(r+ 2-4 ll -(Yz+ | : 
T TR Yr 
Diese ungerade Funktion in Yr multiplizieren wir mit der ebenfalls in Yr ungeraden 
Funktion 


ll Fe 1 
— (arctgyr-tarct = 
7 Ir 


die andererseits in einem kompletten r-Blatt, auf das wir uns beschränken dürfen, identisch 
gleich Eins ist. 


Mit diesem Faktor läßt sich die verlangte additive Zerlegung dann durchführen: 


* (74 : ) (are t Yr + arctg - ) 
= FA} = _- — (02 arc I ST 
2 Yr a Yr 


3 j# 
-(fr+,) arctg fr + ır 
T 


regulär im Innengebiet 


r 1213 ar 1 
++) eyiasit ne 


regulär im Außengebiet, 
5 
mit — ?% — multipliziert, 
7 
gleich der gesuchten Strömungsfunktion. 


Eine Verallgemeinerung läßt sich durch Hinzunahme eines meromorphen im Außen- 
gebiet regulären Faktors f (r) zur Abbildungsfunktion anbringen. — Das führt zu: 
A 
= (Ve + a) fe 


yr 
und ist ebenfalls aufspaltbar. 


Durch passende Wahl des Faktors f(r) läßt sich jede Kontur mit rechtem Winkel 
beliebig genau annähern. 


Sehrifttum 


[1] G. v. Gorup, Eine neue Methode zur Berechnung der Strömungsfunktionen bei zeitlich veränderlicher 
Kontur. Z. angew. Math. Mech. 32 (1952), S. 371. 


Von einem bewegten Störpunkt ausgehende Wellen in einem 
dispersionsfähigen Medium (Wasserwellen) 


Von E. Becker in Göttingen 


Eine punktförmige Störquelle, die in einem schwingungsfähigen Raum mit konstanter 
Geschwindigkeit wandert, erzeugt ein in bezug auf die Störquelle zeitlich stationäres Wellen- 
feld. (Beispiel: Bewegung eines Störpunktes über eine Wasseroberfläche.) Für die Kurven kon- 
stanter Phase (z. B. Wellenberge) in diesem Wellenfeld gelten dabei folgende Beziehungen: In 
einem Aufpunkt mit den Koordinaten x und y (x-Achse in Fortschreitungsrichtung der Stö- 
rung) ist die örtliche Wellenlänge A gegeben durch V(A) —=ccos ©, wobei V(7) die nach Maßgabe 
des für die Wellenfortpflanzung in dem betr. Medium gültigen Dispersionsgesetzes zu A gehörige 
Phasengeschwindigkeit bedeutet, c die Geschwindigkeit der Störung und © den Neigungswinkel 
zwischen Fortschreitungsrichtung der Störung und örtlicher Wellenfront. Außerdem muß die 
zu A gehörige Gruppengeschwindigkeit U(A) einen solchen Betrag haben, daß der Energieflluß 
in dem Feld radial von der Störquelle weggerichtet ist. Aus diesen beiden Bedingungen ergibt 
sich folgende Differentialgleichung für die Kurven konstanter Phase: 


y __U/VsinO cosO (1) 
a a en. 


d 
mit O=arctg 4 . 
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Nach Gl. (1) gehen alle Kurven konstanter Phase durch eine Ähnlichkeitstransformation 
auseinander hervor. 

Es kann nun vorkommen, daß das Wellenfeld auf einen keil- bzw. kegelförmigen Raum- 
teil (je nachdem sich die Störung in einem zwei- bzw. dreidimensionalen Medium bewegt) be- 
schränkt bleibt: Machsches Phänomen. Voraussetzung ist ein Extremum von y/&. Durch 
Differentiation von (1) nach © findet man als Bedingung für ein solches Extremum: 


ec? ya te arme Fe Be ea (2), 
wobei % wie folgt definiert ist: 
U/k 
— 1 a Ban 
vr Nun 8), 


Du 
mit k= FE — Wellenzahl. 


Die Funktion yist eine Transformierte des Dispersionsgesetzes des Mediums. Trägt man 
yin Abhängigkeit von V? in einem Diagramm auf, so erhält man ein einfaches Kriterium für 
Mach sches Phänomen, wenn man in diesem Diagramm im Abstand c? von der Abszisse eine 
Parallele zu dieser zieht. Schneidet diese Gerade die y-Kurve, so tritt das Machsche Phä- 
nomen auf. Es läßt sich zeigen, daß auf den durch diese Schnittpunkte gegebenen Extremal- 
geraden im Wellenfeld die Kurven konstanter Phase Spitzen bilden und daß die Amplituden 
dort über alle Grenzen wachsen. Dies folgt aus der nach der Kelvinschen Methode der 
stationären Phase abgeleiteten Formel für die Wellenamplituden n: 


TE aa | 
_4A[lk yı De j 2m: U/k R Rn; 
1 m e En | 22) (© y) E dU/Ar e.\w tische Ir (4). 


A ist eine Dimensionskonstante, R der Abstand des Aufpunktes vom Ort der Störung, Ö der 
Winkelabstand von der Fortschreitungsrichtung, @ die Phase im Aufpunkt, n —=2 bzw. 3 die 
Zahl der räumlichen Dimensionen. Man erkennt außerdem, daß die Phase in einer Spitze 
(ce? =y) um »/2 springt, da der Nenner unter der Wurzel dort sein Vorzeichen ändert. 


Ein Referat des Vortrages: 


D. Riabouchinsky, Einige hydrodynamische Anwendungen verallge- 
meinerter komplexer Größen 


ist nicht eingegangen. 


F. Strömungslehre und Gasdynamik 


: Überschallströmungen um schlanke Drehkörper 
Von ©. Heinz in Lörrach/Baden 


Bei der Berechnung der Potentialströmung um schlanke, gegen die Strömungsrichtung 
schwach angestellte Drehkörper führt die Randbedingung der tangentialen Umströmung auf 
gewisse lineare Integralgleichungen für die Belegungsfunktionen f} bzw. fs, durch die das Poten- 
tial einer Quell- bzw. Dipolverteilung dargestellt werden kann. Im Falle einer Überschall- 


strömung sind diese Integralgleichungen vom Volterraschen Typ und haben folgende 
Form: . 


Ao 
„Sfr — Ay Sol 7) Sf Ad=s, y=kla)  r=l3RBe ei: 
Dabei bedeutet: 


x, y: Zylinderkoordinaten mit Symmetrieachse des Körpers als &-Achse; 2 —=0, y =0 
ist die Spitze des Drehkörpers. 


A =YM®— 1=cot a;a—=Mach scher Winkel der ungestörten Strömung. 
y =k(x) ist die Gleichung der Meridiankurve des Drehkörpers. 


45 - N, 


s, sind durch 4 = k(x) bestimmte Funktionen und somit bekannt. 


N 
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Für die Druckverteilung sind weiterhin, nachdem fi und fe aus (1) ermittelt sind, die Integrale 
2 
= ta Ayo NER year)... ......@ 


zu berechnen. 
Um Integralausdrücke der Form 


As 
L,()= [Ka — AyCof 7) Cor 2a 


zu diskutieren, führen wir 2=x— Ay, e = Ay/z als neue Variable ein und finden mit 


Z=e(lofi— 1) 


1 
1 
L, l — em. FE} = 5 € ——. 2 
()=e Sie dr: R,(s,ö)dE &= ed), 


wobei der Kern 


K, (0) = I 


Gelingt es nun, die Kerne K,(e,&) im ganzen Integrationsintervall O<£<1 durch ein 
Polynom in & zu approximieren, 


ist. 


Kyle, 8) 3 Burelo)le - a BUS TAN en, 


so wirdin der Genauigkeit dieser Approximation 


1 
L.() =&7 Bere N 2 — 26)E° ? de 
oder mit T =z£ a 
L, (23. & Burele) z-all(2)e20 3] 
Für die Faltungsprodukte der rechten Seite findet man 


\(z) 03 — Pet) md (2), m(z)=Uz) * 2%, 


7($) 
ML 1 
so daß mit rn „re Yu,roSchließlich 
ee Dr AOL 16 Be 
0o=0 


wird. Hierdurch ist der Integralausdruck auf einen linearen Differentialoperator zurückge- 
führt. 
Die Approximation von K,(e,£) durch ein Polynom kann im Falle e>% durch die 


Ta ylorentwicklung von (2 Ee+ 3 bewirkt werden. Da aber gerade für schlanke Dreh- 
körper & klein sein soll, müssen wir K,(e, &)im Bereich 0 <Z < 1 durch orthogonale Polynome 
D,(£) (hypergeometrische Polynome) approximieren: 


R,(e, 2) 6100) Bel) ran... (5). 


Diese Reihe brechen wir beio =xab und ordnen den so entstehenden endlichen Ausdruck nach 
Potenzen von £ um, wodurch wir (3) gewinnen. 

Die Wahl von x richtet sich nach der geforderten Genauigkeit und nach &; denn für e =0 
kann (4) nicht mehr konvergieren, für kleine e sind demgemäß für eine geforderte Genauigkeit 
mehr Glieder der Entwicklung (5) erforderlich. Eine genauere Analyse zeigt jedoch, daß auch 
für sehr schlanke Drehkörper (e< 1) die Approximation sehr gut ist. 

Mit Hilfe von (4) sind die Integralgleichungen (1) auf Differentialgleichungen zurückge- 
führt, die man (etwa nach dem Verfahren von RungeundKutt a) numerisch integrieren 
kann. Auch die Berechnung der Druckverteilung nach (2) bietet dann keine Schwierigkeit. 

Die Reduktion der Integralgleichungen (1) euf Differentialgleichungen bietet einige Vor- 
teile gegenüber der unmittelbaren numerischen Lösung der Integralgleichungen, die in irgend- 
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einer Form darauf hinausläuft, die Integrale in (1) durch endliche Summen zu ersetzen und so 
ein System von linearen Gleichungen herzustellen: 


1. Die Auflösung eines Systems von n linearen Gleichungen mit n Unbekannten ist sehr 
ungenau, wenn n nicht gerade sehr klein ist. Das ist auch dann der Fall, wenn, wie in unserem 
Beispiel, die Gleichungen gestaffelt sind. Die Schwierigkeiten bei einer numerischen Lösung 
einer Integralgleichung lassen sich mit den Schwierigkeiten bei der Differentiation einer will- 
kürlichen Funktion vergleichen, während die Integration einer gewöhnlichen Differentialglei- 
chung mit der Integration einer willkürlichen Funktion verglichen werden kann. 

2. Der Rechenaufwand bei der numerischen Lösung einer Integralgleichung wächst qua- 
dratisch mit der Anzahl der Schritte an, während er bei der Integration einer Differentialglei- 
chung nur linear mit der Anzahl der Schritte anwächst. 

3. Die rechten Seiten s, der Integralgleichungen (1) haben denselben Unstetigkeitscharak- 
ter wie die aus den Differentialgleichungen zu ermittelnden Funktionen, während die Unstetig- 
keitsordnungen der f, dagegen um #$ verschoben sind. Dies erleichtert wesentlich die Behand- 
lung von Drehkörpern, deren Konturen Knicke aufweisen. 

Die Werte der y„,‚.(e), die als Koeffizienten unserer Differentialoperatoren auftreten, 
sind berechnet und werden in einer ausführlichen Veröffentlichung mitgeteilt werden. Das Ver- 
fahren kann auf nichtstationäre Vorgänge (langsam pendelnde Drehkörper) ausgedehnt werden. 


Über eine Grenzlösung für schwinsende Flügel kleiner Streckung *) 
Von H. Merbt und M. Landahl in Stockholm 
Für Flügel sehr kleiner Streckung kann man in der linearisierten Potentialgleichung 


2M 1 
ZN Das 5 Pi ll. Yas i r (1). 


(1 a M%) Dr w Dyy . D;; — 


Änderungen des Potentials in &#-Richtung vernachlässigen (vgl. Jones’ Stationäre 
Theorie [2]). Es wird dabei angenommen, daß die Störbewegung wesentlich zweidimensional 
ist in Schnitten senkrecht zur Anströmrichtung. Die Änderungen des Flügelquerschnittes in 


Flügelrichtung gehen dann ein bei der Berechnung des Drucksprunges und der Kräfte am 
Flügel. 


Mit dieser Annahme und harmonischem Ansatz 


RE A es 
geht Gl. (1) über in die Wellengleichung 


2 
YatBat+tz Bliiechnmepi®. 222.220) 


Gl. (3).ist zu lösen unter folgenden Randbedingungen 
a) 0\P/öz soll am Flügel vorgegebene Werte annehmen entsprechend der vorgegebenen 
Schwingungsform. 


b) in der Flügelebene außerhalb ces Flügels ist Y—=0 wegen der Antisymmetrie des 
Potentials. 

)V,YP,P,>0 für YP+z2—o. 

d) die Lösung genügt der Sommerfeldschen Ausstrahlungsbedingung. 

Transformiert man Gl. (3) und die Randbedingungen auf elliptische Koordinaten, so 
läßt sich das sich ergebende System lösen mit Hilfe eines Separationsansatzes. Der Lösungs- 
weg entspricht dm Timmans[3]. 

Das Ergebnis läßt sich schreiben in Form einer unendlichen Reihe inMathieuschen 
Funktionen. Die Reihenkoeffizienten folgen aus Randbedingung a). Die Lösung gilt sowohl 


im Unterschallgebiet als auch im Überschallgebiet. Für kleine Schwingungsfrequenzen geht 
sie in die bekannte stationäre Grenzlösung über. 


Abschätzungen durch Einsetzen der Lösung in die ursprüngliche Gl. (1) ergeben Bedin- 


gungen über die zugelassenen Bereiche von aspect ratio A, reduzierte Frequenz ne 


oo 


*) erscheint in der Serie Technical Notes der Technischen Hochschule Stockholm [1]. 
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und Machzahl der Anströmung M,„ (list die ganze Flügeltiefe). Für Schlagschwingung folgt 
beispielsweise 
(1 — M3) 4?< 16 


VE 
>20. 


I 


Rechnungen wurden durchgeführt für Dreiecksflügel. 


Schrifttum 
[1]H.MerbtundM.Landahl, Aerodynamic forces on oscillating low aspeet ratio wings in compressible flow. 
KTH-AERO TN 30 (1953). 
[2] R.T. Jones, Properties of low-aspect-ratio wings at speeds below and above the speed of sound. NACA 
TN 1032 (1946). 
[3] R. Tim man, Beschouwingen over de luchtkrachten op trillende vliegtuigvleugels. Diss. Delft. (1946). 


A PROPOS du POINT d’IMPACT du CYLINDRE CIRCULAIRE 
et de PEPAISSEUR d’IMPULSION 


Adalbert Oudart, Paris 


1-I s’agit!) de la determination de l’&paisseur d’impulsion ö,(0) (Impulsverlustdicke) 
au point d’impact (Staupunkt) pour un fluide incompressible en &coulement plan dans le cas 
ou les equations de Prandtlsont appliquees jusqu’au bord d’attaque. 


2- On etablit que pour le cylindre eirculaire erayon R,ona: 


(0) = d, 2 (d, — coefficient nume6rique fee... 4,0%) 


x 


alors que pour la plaque normale au vent (Hiemenz, Falkner et Skan): 


2 = Bm au = 3,0) Yrlimite / 2. en) 


Us x—0 


avec ö,(1) —= 0,292; x = abcisse a partir du point d’impact; u,= vitesse de l’&coulement 
potentiel. 

On constate que pour R=w, (A) fournit ö&,—=, alors que (B) fournit ö, fini. Cette 
anomalie n’est qu’une apparence, car il faut prendre: uv,—=TU, y=», de sorte qu’il faut 
faire augmenter «, comme R. 


© 
On a d’ailleurs: u, (2) =2 us sin 17 2 R® de sorte que l’on peut ecrire (A) sous la 
forme: 


,(0)— 5, Y2Yr imite/ 2 Bee. |. ...KAN), 
x—0 Ue 
On en conclut qu’iln’ya pas d’anomalie et que: 
6, = 5,(1)/Y2 = 0,292/Y2. 
3- Pour un profil quelconque A bord d’attaque arrondi,ona dans lamöthoded’Eckert 
(Profils d Falkner et Skan — loc. cit. (2) p 63): 
ö,(0) — 0,292 Y v limite Vz en (Ch 
x—0 Ue 
Dans la möthode de Karl Pohlhausen (Profils P4, loc. eit. (1) p 196), on a 
20.077 A’ 0u: Me 
ö,(0) — 0,278 Yalmite / 5 N rare re (Di) 
x>0 Ue 
1) References: 


(1) H.Schliehting, Grenzschichttheorie (1951). 
(2) E.Brun, Cours du C.E.S.M. (1947). 
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I] ya une bonne concordance. 

La methode rigoureuse (loc. eit. (1)pp 125—132) conduit d’ailleurs aisement & la for- 
mule (C). 

4- Pour un profil aA bord d’attaque tranchant et point d’impact en ce bord, on justifie 
par similitude et passage & la limite que: 


5, (0) =, (0) =0. 


Über die maximale Ablenkung im gegabelten Verdichtungssstoß 
Von W. Kofink in Karlsruhe 


Im Weise diagramm für Gabelstöße in Luft (vgl. z. B. Bild 2 der Abh.?)) erkennt man 
im 2. Lösungsgebiet (rechts der ö — ö,-Geraden) zwei verschiedene Arten von ® =const.- 
Kurven, die mankurz als ‚einfach gekrümmte‘‘ für | < 2,038und ‚‚geschlängelte‘ ‘für o, > 2,038 
charakterisieren kann. Wir erörtern an Hand des Bildes 1, wie diese Gestaltsänderung der 
Kurven mit wachsendem @, zu- 
standekommt. Die Maxima in 
ö-Richtung der ‚einfach gekrümm- 
ten“ (imBild1 nicht eingezeichneten) 
©, = const.-Kurven liegen für 


1,749 < 0, < 2,038 


auf der Kurve RS; der zugehörige 
Abszissenwinkel ö ist der maximale 
Ablenkungswinkel im Gabelstoß 


Omar), der unterhalb o, = 2,038 
kleinerist als der maximale Ab- 
lenkungswinkel im einfachen schie- 


fen Stoß Omas (@1). Diese Öhaz(}) 
sind eine Folge der in $2f der Abh.?) 
beschriebenen 2. Ursache: (Om/dx), 
—=0. Mit wachsendem o, wird jedoch 

Äurve, worauf j 2 Fr 
(mk)m2+(3-Ym-2l=0 die Differenz Ömaz(®1) — Ömaz(®,) 
immer kleiner und verschwindet in 
S für o, = 2,038. Für größere w,, 
nämlich 2,038 < @, < Omas = 23,437, 
gibt es dagegen für jedes &, zwei 
Zunkis mit nn Ordina- 
ten 6% ’(w,) und 64.(@,) bei gleicher 
Wr Abszisse ö= Öyan(@,), wo in den zu- 
/ gehörigen 2 Gabelstoßkonfiguratio- 
nen der maximale Ablenkungswinkel 
des einfachen schiefen Stoßes er - 
reicht wird. ö®(w,) bzw. 6(w,) 
liegen auf den Kurventeilen SU bzw. 
ST. Diese ö-Maxima folgen aus der 
ersten in Il $2f beschriebenen Ur- 


: : U. säche: ‚II Gl. (7). Auf. der” Fort- 
Bild 1. Wenn man sich die Kurven RSZ und UST dieses Bildes indas Weisedia- 7 B 
gramm für Gabelstöße in Luft eingezeichnet denkt, so schneiden die Kurven kon- setzung SZ der Kurve RS liegen nun 


stanter Anströmgeschwindigkeit &, = const. jenes Diagramms unsere Kurven 1 1 1 -Richtuns der w, = 
parallel zur ö,-Achse. Die an den Kurven des Bildes besonders markierten Minima in ö-Riec 8 1 


Punkte bezeichnen die Schnittpunkte mit den entsprechenden ®, = const.- const.-Kurven für ®, > 2,038. Das 
Kurven des Weise diagramms. „Auseinanderspreizen‘‘ der Maxima 

zusammen mit dem Auftreten der Mi- 

nima für ©, > 2,038 verleiht den (im Bild 1 nicht eingezeichneten) geschlängelten &, = const.- 
Kurven für wachsende &, immer stärkere Windungen. Schließlich rückt das ö-Minimum für 
©, = @paz in den Punkt Z, welcher mit der „‚Ecke“ 28 in II Bild 9 identisch ist, und die beiden 
Maxima rücken nach U bzw. T d.h. nach den Punkten 34 bzw. 22 in II Bild 9. Das Durch- 


1)W.Kofink, Ann.d. Phys. 9, 200, 1951. Hier mit I zitiert. 

2)W.KofinkundTh.Vollmer,ZaMM 33, 73, 1953. Hier mit II zitiert. 

Druckfehlerberichtigung inII: Auf $. 88, Zeile 19 von oben sollte im Koeffizienten von 1? 
20 %k + k? statt 20 k? + k steben. 
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laufen von Ohaz(0,) auf den Kurventeilen RS und SU bei positiven ö, macht deutlich, daß die 
‚„‚Fußpunkte‘ (s. II $2a) aller &, =const.-Kurven im Weise diagramm für 1,749 < 
< 2,437 auf der Strecke RU bei kleineren ö-Winkeln liegen als die Abszissen jener Maxima. 
Alle , —=const.-Kurven dieses w,-Intervalls laufen also unter einem spitzen Winkel aus der 
Abszisse aus. Die Kurven für &, =2,0 und 2,2 in I Bild 2 müssen daher vor ihrem FEinlaufen 
in die ö-Achse ein wenig zurückbiegen und diesem Sachverhalt entsprechend korrigiert werden. 


Punkte auf RSZ befriedigen simultan die GIn. II (5) und (17), Punkte auf UST die Gln. II 
(5) und (7). Der Übergangspunkt S befriedigt als Schnittpunkt beider Kurven simultan alle 
3 GIn., woraus seine Daten für Luft (x —= 1,405) I; = 10,4891, m; = 8,5177, &; =1,2321; 
(@1)s = 2,0379, (@,,)s = 1,9773, (on, .)s = 1,3596; ds — 36°11,6’, (ö,)s — 11°16,6’ folgen. 
Nach Srücken für ©, > (®,)s nicht nur die ö-Maxima und Minima, sondern auch die Wende- 
punkte der geschlängelten Kurven hinein: dö ist in S auf der ® =(w,)s-Kurve proportional 
döa)*. 
Ein Vergleich mit dem 1. Lösungsbereich ergibt, daß den Stoßkonfigurationen der Kurve 
RSZ dort diejenigen senkrechter Hauptstöße entsprechen; in beiden Fällen ist (9m/dx); —0. 
Schließlich sei daran erinnert, daß im 1. Lösungsbereich für jedes &, > 1,2448 der maximale 
Ablenkungswinkel des einfachen schiefen Stoßes auch im Gabelstoß erreicht werden kann, 
nämlich auf Kurve (, in Bild 2 der Abh. II. Unserer Kurve RSZ dagegen entspricht dort die 
Kurve (,. 


Die Strömung um rotierende Körper bei axialer Strömung 
Von E. Truckenbrodt in Braunschweig 


Führt ein axial angeströmter Körper eine Drehbewegung um seine Achse aus, so entsteht 
ein die Drehbewegung hemmendes Drehmoment, das abhängig ist vom Verhältnis der Um- 
fangsgeschwindigkeit zur Anströmungsgeschwindigkeit. Die Ursache hierfür ist in den Vor- 
gängen der wandnahen Reibungsschicht zu suchen, die infolge der Drehbewegung des Körpers 
mitumläuft. Diese Reibungsschicht hat dreidimensionalen Charakter und läßt sich mit den 
Methoden der Grenzschichttheorie behandeln. Der Rechnung am einfachsten zugänglich sind 
ä Pl k-0:Faden 

Z ; k-1: Kugel 


f 
as 

== =mn 
a a 


turbulent 
41 1 Ei 


Die Drehmomentenbeiwerte von axial angeströmten rotierenden 
R.otationsellipsoiden. 


die Fälle der angeblasenen rotierenden Scheibe und des angeblasenen rotierenden Drehkörpers. 
Während man für die angeblasene rotierende Scheibe den ganzen Bereich von der angeblasenen 
stillstehenden Scheibe bis zur nicht angeblasenen rotierenden Scheibe erfassen kann, gilt dieses 
für den angeblasenen rotierenden Drehkörper nur für mäßige Verhältnisse von Umfangsge- 
schwindigkeit zu Anströmungsgeschwindigkeit. Für den Jaminaren Strömungszustand 
wurden solche Rechnungen für die angeblasene rotierende Scheibe von Fi$Sc h lichting 
und E. Truckenbrodt [1] und für den angeblasenen rotierenden Drehkörper von H. 
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Schliehting [2] durchgeführt. Für den turbulenten Strömungszustand wurden 
die beiden Fälle von E.Truckenbrodt [3], [4] behandelt. Die hierbei gewonnenen For- 
meln ergeben bei sinnvoller Übertragung auf den laminaren Fall hierfür recht brauchbare Nähe- 
rungsformeln. 

Mit den Vereinfachungen der Grenzschichttheorie lassen sich aus den Bewegungsglei- 
chungen in meridionaler (bzw. in radialer Richtung bei der Scheibe) und in azimutaler Richtung 
(Umfangsrichtung) die entsprechenden Grenzschichtgleichungen anschreiben. Aus diesen er- 
hält man durch Integrationen über den senkrechten Wandabstand und unter Berücksichtigung 
der Kontinuitätsgleichung die Impulsgleichungen für die meridionale und azimutale Richtung. 
Während für die Jaminare Strömung ein unmittelbarer Zusammenhang zwischen den Wand- 
schubspannungskomponenten und den Geschwindigkeitsgradienten an der Wand besteht, ist 
dieses für die turbulente Strömung nicht so einfach. Wie in [4] näher gezeigt wird, findet man 
aus bestimmten Analogiebetrachtungen zum schiebenden unendlich langen Flügel brauchbare 
Näherungsansätze für die Wandschubspannungskomponenten in meridionaler und azimutaler 
Richtung. Aus der Impulsgleichung in azimutaler Richtung (Umfangsrichtung) erhält man dann 
unter Einführen des Ansatzes für die Wandschubspannungskomponente in Umfangsrichtung 
eine einfache Quadraturformel zur Berechnung des Drehmomentenbeiwertes eines angeblasenen 
rotierenden Drehkörpers, in der nur die Verteilung des Körperhalbmessers sowie die potential- 
theoretische Geschwindigkeitsverleilung längs eines Meridianschnittes vorkommen. Der lami- 
nare bzw. turbulente Strömungszustand geht ein in Form des Reibungsbeiwertes der längs- 
angeströmten ebenen Platte sowie eines bestimmten Exponenten. Das Drehmoment ist bei 
festgehaltener Anströmungsgeschwindigkeit proportional der Umfangsgeschwindigkeit. Das 
Bild zeigt d'e Abhängigkeit des Drehmomentenbeiwertes von der Reynoldsschen Zahl für 
verschiedene Rotationsellipsoide vom unendlich langen Faden (k=0) bis zur Kugel (k=1). 


Sehrifttum 


[1]H.Schlichting, E. Truckenbrodt: Die Strömung an einer angeblasenen rotierenden Scheibe 
Z. angew. Math. Mech. 32 (1952), S. 97, vgl. auch Auszug in Journal Aeron. Sci. 18 (1951), S. 638. 

[2])H.Schliehting: Dielaminare Strömung um einen axial angeströmten rotierenden Drehkörper. Erscheint 
demnächst im Ingenieur-Archiv, vgl. auch Auszug in Wissenschaftliche Zeitschrift der T. H. Dresden, 2 
(1952/53), S. 423. 

[3] E. Truckenbrodt: Die Strömung an einer angeblasenen rotierenden Scheibe bei turbulenter Strömung. 
Erscheint demnächst in Z. angew. Math. Mech. 

[4] E.Truckenbrodt: Ein Quadraturverfahren zur Berechnung der Reibungsschicht an axial angeströmten 
rotierenden Drehkörpern. Erscheint demnächst im Ingenieur-Archiv. 


Über den Strömunsswiderstand von Sieben 
Von K. Wieghardt in Hamburg 


Der Strömungswiderstand von Sieben ist in den letzten zwanzig Jahren öfters gemessen 
worden, da es anscheinend unsicher war, welchen Druckverlust ein gegebenes Sieb verursachen 
würde, wenn es einem schon gemessenen Sieb nicht genau entsprach. Beschränkt man sich auf 
die meistens benutzten Runddrahtsiebe mit quadratischen Maschen, so ist das Sieb geometrisch 
bestimmt durch die Drahtstärke d und die Maschenweite Z, bzw. durch das Öffnungsverhältnis 


offene d\? 
- 1 


gesamte Fläche in 
Die einfachsten Randbedingungen für die Siebströmung werden erzwungen, wenn das Sieb 
senkrecht in einem Rohr angeordnet wird. Dann ist die Geschwindigkeit vor und hinter dem 
Sieb gleich U und nur der statische Druck fällt von dem Wert p, vor dem Sieb ab auf den Wert 
?, dahinter. Wie üblich sei die Druckverlustziffer definiert durch 


k=(m—mj% U:  (e=Dichte), 


Offenbar muß diese Zahl k von ß abhängen und vermutlich auch noch von irgendeiner Re-Zahl, 
etwa von Ud/» (» = kinem. Zähigkeit). 

Zum Ordnen der Meßergebnisse liegt es nun nahe, den Widerstand des ganzen Siebes zu 
vergleichen mit dem des Einzeldrahtes, also eines Zylinders, und OÖ. Flachsbart[l] hat 
dies qualitativ für sehr offene Siebe bereits 1932 getan. Man kann diese Reduktion aber auch 
quantitativ noch weiter durchführen. Zunächst könnte man dazu so vorgehen: der Anteil 
(1 — ß) des gesamten Rohrquerschnittes wird durch Drähte versperrt, deren Widerstands- 
beiwert ce noch von Ud/» abhängt; dann müßte der Widerstand des Siebes k = (1 — ß)c sein, 
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Diese Formel ist aber offensichtlich falsch, da für ß > 0, also für ein sehr dichtes Sieb, der 
Widerstand oder % beliebig groß werden muß. 

Um die Abhängigkeit von ß oder die gegenseitige Beeinflussung der Drahtstücke besser 
zu erfassen, muß man berücksichtigen, daß der Widerstand eines Einzelzylinders fast aus- 
schließlich von den Zuständen im Totwasser herrührt. Bei den hier in Frage kommenden unter- 
kritischen Re-Zahlen ist die Grenzschicht laminar und löst sich etwa bei & — 80° ab (& Zentri- 
winkel vom Staupunkt aus gemessen), und es entsteht dort ein Totwasser mit Unterdruck. 
Angewandt auf das Sieb heißt das, daß die Geschwindigkeit etwa bei x = 90°, also im Sieb 
selbst, wesentlich sein wird für die Ausbildung des Totwassers und damit für den Widerstand. 
Die mittlere Geschwindigkeit dort ist aber aus Kontinuitätsgründen größer als vor oder hinter 
dem Sieb, nämlich gleich U/ß. 

Wenn man nun diese Geschwindigkeit dimensionsanalytisch als charakteristisch ansieht 
anstelle von U, so wird jetzt 


i 1 
hı -m=(l—P)e 20 (U/ß)°, und nach Definition von k, PL — Ps ko oU®, 


Daraus folgt | 
k= = P O4 
P? 


wo jetzt ce noch von (U/ß)d/v abhängen kann, aber sonst nicht mehr von ß. Tatsächlich fallen 
nun alle bisher veröffentlichten Messungen gut auf eine Straße, wenn man 


o=* E k über (U/ß) d/v 
1—Bß 

aufträgt. Und zwar bestimmen gerade die besonders genauen Messungen von Simmons 
undCowdrey[2] sehr gut eine einzige Kurve hierfür, während andere Messungen stärker 
streuen. Obwohl bei dieser Reduktion der Einfluß der Überkreuzung und Verbiegung der 
Drähte im Sieb außer acht gelassen wurde, liegt diese Kurve überraschend nahe bei der des 
Widerstand beiwertes c, für den Einzelzylinder in unendlicher Flüssigkeit über der Re-Zahl. 
Für 60 < (U/ß) d/v < 600 kann man die Kurve annähern durch 


> Ud —1/3 
2 (75) 


Dann liegen ce zwischen 1,5 und 0,7 und das ursprüngliche % zwischen 0,3 und 17 für Siebe mit 
Öffnungsverhältnissen zwischen ß = 0,2 und 0,8. 


Schubauer, Spangenberg und Klebanoff [3] haben ferner festgestellt, 
wann ein Sieb als Turbulenzerzeuger wirkt. Sie beobachteten mit einem Hitzdraht dicht 
hinter dem Sieb, bei welcher kleinsten Re-Zahl der Oszillograph turbulente Schwankungen 
anzeigte. Auch dieses Ergebnis kann zusammengefaßt werden, wenn man die Re-Zahl (U/ß) d/v 
als charakteristisch ansieht. Turbulenz wird dicht hinter dem Sieb dann erzeugt, wenn 
28 für Siebe mit verschiedenen Öffnungsverhältnissen ß. Auch nach den Ho- 

v 


mannschen Strömungsaufnahmen [4] an Einzelzylindern würde man eine Re-Zahl von etwa 
dieser Größe als kleinste Re-Zahl ansehen, bei dereine Kärmän sche Wirbelstraße entsteht. 
Das bisher Berichtete wird demnächst im englischen Aeronautical Quarterly (1953) ver- 
öffentlicht werden. j 
Weiterhin wurden noch Messungen an Parallelgittern aus Rundstäben durchgeführt, 
(Hier ist 8 = 1 — d/L, mit L = Gitterteilung.). Die Messungen geschahen in einem kleinen 
Gebläse und waren nicht sehr genau; immerhin ergab sich, falls nur genügend viele Stäbe über 
dem Rohrquerschnitt angeordnet waren, e »# const = 0,8. Ein Einfluß der Re-Zahl war inner- 
halb der Streuung der Meßergebnisse nicht zu erkennen; die Re-Zahlen entsprachen dem 
Bereich, in dem auch der Widerstandsbeiwert des Einzelzylinders sich nur zwischen 0,9 und 1 
bewegt. 
Um weitere Einzelheiten der Strömung zu erkennen, wurde die Druckverleilung um 
einen Rundstab im Gitter gemessen, s. Bild 1. Durch Integration wurde daraus der Druck- 
widerstand c, eines einzelnen Stabes im Gitter berechnet. (Der Reibungswiderstand ist bei 
diesen Re-Zahlen bereits vernachlässigbar klein.) Wie aus der Tabelle in Bild 1 hervorgeht, 
war tatsächlich ß?c, = const = 0,8, wodurch der obige Wert für den Widerstand des ganzen 


Gitters bestätigt wird. 
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Schwieriger zu übersehen sind die Verhältnisse bei Parallelgittern aus quadratischen 
Stäben oder Blechstreifen. Nach Betzund Petersohn [5] ergibt sich hier eine Instabi- 
lität der Strahlen hinter dem Gitter. Zunächst ist das Totwasser breiter als die Stäbe, und falls 
zwei Strahlen durch eine zufällige Störung 
sich nähern, so saugen sie sich weiter an 
und bilden einen großen Strahl. Dies führt 
zu Unsymmetrie und schließlich zum An- 
saugen der Strömung an einer der Rohr- 
wände, außer wenn man solche Vorgänge 
durch besondere Führungsbleche ver- 
hindert. 

Bei Doppelgittern aus zwei Parallel- 
\ Zylinder gittern aus quadratischen Stäben, die hin- 

N“ (Potential-Theorie) 7 tereinander und senkrecht zueinander be- 
B% A festigt sind, scheint diese Instabilität je- 
doch nicht aufzutreten, vermutlich weil 
hier die Totwassergebiete der einen Stab- 
gruppe über die der quergestellten anderen 
Stäbe belüftet werden. Messungen an sol- 
chen Doppelgittern in Wasser von Baines 
und Peterson [6] sind in der folgenden 
Tabelle umgerechnet. (Die Seitenlänge der 


Druckverlustin einem 
Doppelgitter aus qua- 
dratischen Stäben 


nach [6] 
ß k [ 
0,111 116 1,61 
0,250 19 1,58 
0,391 6,2 1,56 
Bild 1. Druckverteilung er in verschieden dichten SE nn ne 


quadratischen Stäbe war d = 1’’ —= 25,4 mm.) Die reduzierte Druckverlustziffer c ist danach 
fast konstant und liegt in der Nähe des Widerstandsbeiwertes eines Einzelstabes (2,0 für ebene 
Strömung), während die ursprüngliche Ziffer k von 116 auf 0,58 abnimmt. 

Die Messungen von Flachsbart[1] an Doppelgittern aus Blechstreifen zeigen leider 
in der Auftragung c über ß zu große Streuung, da diese Auftragung viel empfindlicher ist als die 
frühere % über ß. 
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[6] W.D.BainesundE.G.Peterson, An Investigation of Flow through Screens. Transactions of the Ame- 
rican Society of Mechanical Engineers Bd. 73 (1951), S. 467. 


Überdie Differenzenverfahren zur Berechnung laminarer Grenzschichten 
Von H. Witting in Freiburg i. Br. 


Zur Lösung der Prandtlschen Differentialgleichungen ebener, stationärer, Jaminarer 
und inkompressibler Grenzschichtströmungen bei vorgegebenen Randbedingungen wurden von 
K. Schröder (Math. Nachr. 4) und H. Görtler (Ing.-Arch. 16) Differenzenverfahren 
entwickelt, die im wesentlichen darauf beruhen, daß man die Differentialgleichungen in den 
Punkten eines rechtwinkligen Gitternetzes durch Differenzengleichungen ersetzt. Diese 
Näherungsgleichungen lassen sich unter den Randbedingungen rekursiv auflösen, und man kann 
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so durch Fortsetzen vorgegebener Einlaufprofile in Strömungsrichtung Näherungslösungen 
gewinnen. 


Die Berechnung eines neuen Geschwindigkeitsprofils erfolgt bei K. Schröder durch 
die Nacheinanderausführung zweier Iterationsprozesse; von diesen kann man zeigen, daß sie 
gegen die Werte konvergieren, die bereits das einfacher arbeitende Verfahren von H. Gört- 
ler unmittelbar liefert. 


L.Prandtl(ZaMM 18) hat durch Umformung der Differentialgleichungen einen strengen 
Ausdruck zur Lösung der Fortsetzungsaufgabe gewonnen. Neben anderen Schwierigkeiten bei 
der praktischen Durchführung ist es besonders die, daß der Fortsetzungsausdruck an der 
Wand des umströmten Körpers unbestimmt von der Form 0/0 ist, und sich somit die Ge- 
schwindigkeitsverteilung in Wandnähe bei ungünstiger Fehlerfortpflanzung errechnet. 


Diese Linie der Unbestimmtheit ist es auch, die die Genauigkeit der genannten Differen- 
zenverfahren sehr beeinträchtigt. Es zeigt sich nämlich, daß sie infolge der finiten Ersetzungen 
in das Innere des Strömungsgebietes verschoben ist; insbesondere geht sie bei Annäherung an 
die Ablösungsstelle in die der Wand benachbarte Gittergerade über: Der Fortsetzungsaus- 
druck zur Berechnung des wandnächsten Wertes wird also bei Annäherung an die Ablösungs- 
stelle unbestimmt. Dieses erklärt den Genauigkeitsverlust der genannten Verfahren in diesem 
Gebiet. 


Diese für die numerische Rechnung ungünstige Verschiebung wird vermieden, wenn man 
die Geschwindigkeitsverteilung in Wandnähe geeignet durch Polynome approximiert. Indem 
wir ferner die finiten Approximationen in diesem Gebiet besonders genau vornehmen, tragen 
wir der Unbestimmtheit des neuen Fortsetzungsausdrucks an der Wand Rechnung. Die neuen 
Ersetzungen bedingen bei jedem Fortsetzungsschritt die Durchführung eines Iterationsprozes- 
ses, der theoretisch sehr durchsichtig und praktisch überaus einfach zu handhaben ist. Trotz 
der genannten Abänderungen bleibt der Grundcharakter des ursprünglichen Verfahrens 
erhalten. Die Erhöhung des Arbeitsaufwandes — um 10 Minuten pro Fortsetzungsschritt, 
während dieser bisher 1 Stunde in Anspruch nahm — ist unwesentlich und wird durch die 
erzielte Genauigkeitserhöhung mehr als wettgemacht. Das Verfahren wurde an mehreren 
Beispielen erprobt ; dabei wurden — auch in der Nähe der Ablösungsstelle — gute Erfahrungen 
gemacht. 


Wie bei den genannten Verfahren so zeigt auch die nach dem verbesserten Differenzen- 
verfahren errechnete Näherungslösung eine charakteristische Streuungserscheinung; diese 
Aufspaltung der Lösung läßt sich als eine durch gewisse finite Ersetzungen bedingte Instabili- 
tät des Näherungsverfahrens erklären und durch eine einfache mechanische Glättungsvorschrift 
beheben. 


Wegen der Einzelheiten vergleiche man zwei Arbeiten, die in Kürze in der Zeitschrift 
f. angew. Math. u. Phys. bzw. im Archiv der Mathematik erscheinen. 


Strömung durch Schlitz- und Lochblenden bei sehr kleinen 
Reynoldsschen Zahlen 


Von W. Wuest, in Göttingen 


Bei der Strömung durch Schlitz- und Lochblenden ist bei sehr kleinen Reynolds 
schen Zahlen (R < 2, wobei R auf die Schlitzbreite bzw. den Lochdurchmesser bezogen wird), 
der Druckabfall proportional der durchfließenden Menge. Erst bei größeren Reynolds- 
schen Zahlen bildet sich eine quadratische Abhängigkeit heraus. Für den Grenzfall R —=0 
wird eine exakte Lösung der Navier-Stokesschen Differentialgleichungen mit auf der 
Ein- und Ausströmseite symmetrischen Stromlinien hergeleitet und eine Formel für den 
Druckabfall aufgestellt. Für endliche Reynoldssche Zahlen existieren nur Näherungs- 
lösungen, die in unmittelbarer Umgebung des Schlitzes nicht mehr gelten. Auf der Ausström- 
seite löst sich die Strömung bei R> x? ab und bildet einen laminaren Strahl, der bei Erreichen 
einer kritischen Reynoldsschen Zahl periodische Randstruktur aufweist, die schließlich 
zum Umschlag in einen turbulenten Strahl führt. Der Hauptanteil des Druckabfalles vollzieht 
sich in unmittelbarer Nähe des Schlitzes, sodaß die genannten Näherungslösungen nur eine 
Abschätzung des Druckverlaufes in Abhängigkeit von der Reynoldsschen Zahl liefern. 
Ähnliche Verhältnisse gelten auch für die Lochblende. Der rechnerisch ermittelte Wert für den 
Druckabfall bei sehr kleiner Reynoldsscher Zahl liegt etwa 25%, unter dem von Johann- 
sen (1930) gemessenen Wert. Der Unterschied ist vermutlich auf die konische Ausführung 
der Lochblende bei den genannten Versuchen zurückzuführen. 
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Buchbesprechungen — Eingegangene Bücher — Nachrichten 


Z. angew. Math. Mech. 
Bd. 33 Nr. 8/9 Aug./Sept. 1953 


Referate der folgenden Vorträge: 


VandenDungen, Les variants integraux et l’&quation des ondes, 
K. Jaeckel, Analytische Behandlung gewisser Überschallströmungen, 


IsKamperderfieTtieie 
Statistique de la Turbulence, 

Daran 
eines Zyklonabscheiders 


sind nicht eingegangen. 


Sur la definition de la moyenne dans la Theorie 


Rechnerische Untersuchung der Strömungsvorgänge am Eintritt 


BUCHBESPRECHUNGEN 


C. B. Biezeno (Prof. a. d. Techn. Hochschule Delft) 
und R. Grammel (Prof. a. d. Techn. Hochschule 


Stuttgart, Technische Dynamik. Zweite 
erw. Aufl. 
Bd.1: Grundlagen und einzelne Ma- 


schinenteile. XI + 699 S. m. 413 Abb. 
u. 2 Anhängen. Preis geb. 66,— DM. 
Dampfturbinen und Brenn- 
kraftmaschinen. VIII +4528. m. 
315 Abb. Preis geb. 44,— DM. 
Berlin-Göttingen-Heidelberg 1953. Springer-Verlag. 
Die erste Auflage dieses ausgezeichneten Werkes, 
das seit Jahren vergriffen war, wurde in dieser Zeit- 
schrift Bd. 20 (1940), S. 123/124 ausführlich be- 
sprochen. Der etwas unhandliche Band dieser ersten 
Auflage ist jetzt in zwei Bände, die auch einzeln 
käuflich sind, geteilt. Der erste behandelt die Grund- 
lagen und einzelne Maschinenteile, der zweite Dampf- 
turbinen und Brennkraftmaschinen. Der Charakter 
des Buches ist unverändert geblieben. Es zeichnet 
sich durch seine klare, der Denkweise des Ingenieurs 
angepaßte Darstellung aus, fordert aber, um zum 
wirklichen Verständnis durchzudringen, intensive 
Mitarbeit. Man muß den Verfassern sehr dankbar 
sein, daß sie die neue Auflage einer gründlichen Durch- 
arbeitung unterzogen haben. Dabei ist der Inhalt der 
alten in allem Wesentlichen erhalten geblieben. Auf 
diese Abschnitte im einzelnen einzugehen erübrigt 
sich, da das Buch bisher in Kreisen wissenschaftlich 
arbeitender Ingenieure allgemein bekannt ist und 
gern gebraucht wird. Kleinere Verbesserungen und 
Änderungen finden sich an vielen Stellen. Auf die 
größeren Zusätze, die durch die Entwicklung des be- 
handelten Gebietes in den letzten zehn Jahren nötig 
geworden sind, eine Entwicklung, die vielfach auf die 
Arbeiten der beiden Verfasser selbst zurückgeht oder 
doch auf Arbeiten, die durch sie angeregt wurden, 
soll hier kurz hingewiesen werden. Erwähnt sei in 


Bde2: 


Kap. I die Ableitung eines nichtlinearen Elastizitäts- 
gesetzes, in Kap. III die Behandlung einiger neueren 
Methoden der Praktischen Analysis wie der Schwarz- 
Templeschen beiderseitigen Einschränkung des 
niedrigsten Eigenwertes, des Verfahrens von Hardy 
Cross, des Relaxationsverfahrens von Southwell, ferner 
die Benutzung von Meßstreifen zur elektronischen 
Dehnungsmessung, in Kap. IV die Berechnung der 
Stützmomente bei elastischer Lagerung nach Koiter, 
in Kap. V die schwache Umstülpung symmetrieartiger 
Ringe, in Kap. VI die unendliche Halbscheibe mit 
Randbelastung, die unendliche Platte mit einem 
kreisrunden randbelasteten Loch und die zweimal 
gelochte, auf Zug beanspruchte unendliche Platte, 
ferner die Berechnung der mittragenden Breite eines 
periodisch durch Versteifungsringe verstärkten Zy- 
linders und die Berechnung dieser Versteifungsringe, 
in Kap. VII das Scherproblem bei Stäben, Wellen 
und Kreisplatten, die gedrückten und tordierten 
Schraubenfedern. Im zweiten Band ist nur der Ab- 
schnitt über die Biegeschwingung von Schaufeln neu 
geschrieben und der Abschnitt über die Schwingungs- 
tilgung durch Fliehkraftpendel hinzugefügt. 

Das Buch ist ein ausgezeichnetes Lehrbuch, das der 
Einführung des Anfängers in die höheren Abschnitte 
der technischen Dynamik dienen kann. Da aber die 
einzelnen Teile so geschrieben sind, daß sie ohne 
wesentliche Rückgriffe auf die anderen durchge- 
arbeitet werden können, wird es auch vielfach als 
Nachschlagewerk benutzt werden. Das wird dadurch 
wesentlich erleichtert, daß beiden Bänden der neuen 
Auflage ein ausführliches Sachverzeichnis für das 
ganze Werk beigegeben ist. Das Wiedererscheinen 
des vielgebrauchten Buches wird von allen auf diesem 
Gebiet: arbeitenden Wissenschaftlern und Ingenieuren 
lebhaft begrüßt werden. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen, ausführliche Besprechung bleibt vorbehalten. 


F. G. Trieomi (Professor an der Univers. Turin), 
Equazioni differenziali. 2 Aufl. 353 S. 
mit 26 Abb. Turin 1953. Edizioni Scientifiche Ein- 
audi. Preis geb. 4000 L. 


A.Timpe (em. o. Professor an der T. U. Berlin), 
EinführungindieFinanz-und Wirt- 
schaftsmathematik. 2. verb. Auflage. VI 

® 


-+ 217 S. m. 70 Abb. Wiesbaden 1953. Verlag für an- 
gewandte Wissenschaften. Preis geb. 17,80 DM. 


W. Meyer zur Capellen, Leitfaden der No- 
mographie. IV -178 S. m. 203 Abb. Berlin/ 
Göttingen/Heidelberg 1953. Springer-Verlag. Preis 
geb. 17,40 DM. 


NACHRICHTEN 


Freiburg i.Br. Der 0. Honorar-Prof. an der Universi- 
tät Freiburg i. Br., früherer ord. Prof. der Mathematik, 
an der Universität Zürich, Dr. Ernst Zermelo, 
verschied am 21. Mai 1953 im 82. Lebensjahr. 


Mainz. Die diesjährige Tagung der Deutschen Ma- 
thematiker-Vereinigung findet in Mainz vom 20. bis 
26. September statt. 
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Am 14. Juli 1953 verstarb in Boston (Mass.) im Alter von 
70 Jahren der ordentliche Professor der Aerodynamik und der 


angewandten Mathematik an der Harvard-Universität 


Dr. Dr. h. c. 


Richard von Mises 


Es gibt wohl kaum ein Gebiet der angewandten Mathematik, 


das dieser vielseitige Wissenschaftler nicht durch wertvolle Ar- 


beiten bereichert hat. Im Jahre 1920 — damals Professor der 
angewandten Mathematik an der Universität Berlin — gründete 


von Mises die „Zeitschrift für angewandte Mathematik und Me- 
chanik‘“, die er bis zu seiner Berufung nach Istanbul 1934 leitete. 
Die Ziele, die von Mises der Zeitschrift in einem einführenden 
Aufsatz des ersten Heftes weitschauend gesteckt hat, verfolgt 
sie bis heute. Zu seinem 50. und seinem 70. Geburtstage wid- 
mete sie ihm Sonderhefte. Schriftleitung und Verlag werden 


Richard von Mises ein ehrendes Andenken bewahren. 


Schriftleitung Akademie-Verlag 
Professor Dr. Willers Gabelin 
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